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Abstract 


Despite all the recent progress in the held, nanophotonics is still a step behind nanoelec¬ 
tronics in transmitting information using nanometric circuits. A lot of effort is being 
put into making very elaborate structures that can guide light and control light-matter 
interactions at the nano-scale. The held of plasmonics has been especially successful in 
this. In this thesis, a diherent approach is taken to control the light-matter interactions 
at the nano-scale. The approach is based on considering light and sample as a whole 
system and exploiting its symmetries. Thanks to this new perspective, new phenomena 
have been unveiled. These new phenomena have been developed theoretically and/or 
experimentally and are scattered across this thesis. 

In chapter 2, the theoretical grounds of this thesis are settled. Even though every physi¬ 
cist is familiar with the concept of symmetry, a formalism to systematically describe 
the symmetries of electromagnetic helds is explained. With this formalism, some well- 
known symmetry considerations can be as easily retrieved as some much less intuitive. 
For example, it can be demonstrated that a linearly polarised Bessel beam is not cylin- 
drically symmetric; whilst a circularly polarised Bessel beam is both cylindrically and 
dual symmetric. Furthermore, the mathematical tools to describe non-paraxial elec¬ 
tromagnetic helds are given. Due to the fact that this work deals with sub-wavelength 
scatterers, the light-matter interaction cannot usually be described within the paraxial 
approximation. As a result, the polarisation and intensity prohle of the light beams 
cannot be modihed independently as they are linked via the Maxwell equations. 
Chapters 3, 4 and 5 deepen in the study of Generalized Lorenz-Mie Theory. Using 
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the formalism developed in chapter 2, various new effects are discovered. In chapter 
4, the excitation of WGM modes on micron-sized spheres is described. Indeed, using 
cylindrically symmetric beams, light can be coupled into spherical resonators without 
the use of evanescent coupling. Furthermore, it is shown that the use of cylindrically 
symmetric modes also allows for the enhancement of the ripple structure in scattering. 
Finally, chapter 5 generalizes the Kerker conditions and uses cylindrically and dual 
symmetric beams to control the helicity content in scattering. It is shown that non¬ 
dual materials such as Ti 02 spheres can behave as dual if the correct excitation beam 
and wavelength is used to illuminate them. 

Chapters 6, 7 and 8 are devoted to experiments. In chapter 6, a description of the ex¬ 
perimental techniques used in chapters 7 and 8 is carried out. In particular, the basics 
of Spatial Light Modulators and Computer Generated Holograms are given. Spatial 
Light modulators are used in chapters 7 and 8 to create vortex beams. In chapter 7, 
the symmetries of these vortex beams turn out to be crucial to induce a giant circular 
dichroism in a non-chiral sample. Furthermore, the far-held transmission of vortex 
beams through a sub-wavelength nano-aperture is shown for the hrst time. Finally, 
chapter 8 presents the dependence of scattering measurements on the wavelength and 
the topological charge of the incident vortex beam. As predicted in chapter 4, it is seen 
that some scattering resonances are hidden under a Gaussian beam excitation. These 
resonances can be unveiled when the illumination is a vortex beam. 

Overall, this work shows a number of new effects (theoretical and/or experimental) 
produced by the excitation of symmetric structures with symmetric light. These new 
discoveries will help to provide new ideas and design paths to fabricate new nanopho¬ 
tonic devices such as nano-antennas or nano-resonators. A study of the symmetries of 
the system should always be kept in mind for any photonic device where the spatial 
degrees of freedom and the polarisation cannot be decoupled. 
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“Growing up in the place I did I never was aware of any other option 
but to question everything." 


Noam Chomsky 
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Introduction 


Invented by the Sumerians in ancient Mesopotamia around 3000 BC, [1-3], writing may 
be considered one of the greatest human inventions of all-times. Since then, it has been 
used to store, transmit and manipulate information. Concurrently, the technology em¬ 
ployed in writing has changed over the years. The development of computers, perhaps 
one of mankind’s most signihcant technological revolutions in the last century, was in¬ 
trinsically linked to writing. Despite the principles of modern computers having been 
described by Alan Turing in 1936 [4], the hrst programmable electronic computer was 
not built until 1946, at Pennsylvania University [5]. In the following 6 years, John von 
Neumann and co-workers laid the ground work upon which modern computers rest to 
this day [6-8] . The hrst generation of commercial electronic computers was made avail¬ 
able in the US in 1951, but the giant leap towards our current technologies occurred in 
1955, when the so-called ‘second generation computers’ reached the market [9]. Since 
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Figure 1.1; Transmission data rates of different technologies and critical dimension of 
their active devices. The past refers to the semiconductor electronics and vacuum tubes 
used before the 1980’s, whilst modern semiconductor electronics is the technology used in 
nano-electronics. 


that time, the evolution of computers has been linked to the advent of semiconductor 
electronics. Simply put, this is the technology associated with the control of how the 
electrons move in semiconductor materials^. The majority of our current information 
devices function because of our ability to control electron currents from one spot to 
another one, and re-interpret this current as chains of O’s and I’s, which are used to 
encode information. Since the 1970’s, our means of re-interpreting electron currents as 
chains of O’s and I’s have remained fundamentally unchanged. However, the techniques 
and devices to move electrons around in a controlled manner have changed drastically. 
Of particular note, the size of the electronic devices has decreased by many orders of 
magnitude. This can be observed in Figure 1.1. ‘The past’ includes the electronic 
circuits (using semiconductor electronics, or vacuum tubes) used up until the 1980’s. 
From then onwards, the semiconductor manufacturing processes achieved sizes below 

^The development of semiconductor electronics was made possible thanks to the advances in quan¬ 
tum mechanics, and in particular in solid-state physics. 
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1/im currently known as nano-electronics. In fact, nano-electronics has advanced at a 
very fast pace, consistently following Moore’s Law [10]. At present, InteL is fabricat¬ 
ing transistors at 22nm, and current projections are that the llnm technology will be 
obtained by 2015. 

The fact that electronics has been so successful over the past 60 years does not imply 
that it is the ultimate technology to carry out digital writing or information processing. 
Photonics has been following a delayed but similar development. The term photonics 
began to be used in the late 1960’s, once it became evident that photons could be used 
as bits of information in the same way as electrons were being used. However, it did 
not reach the common use until the 1980’s, when the IEEE Lasers and Electro-Optics 
Society established an archival journal named Photonics Technology Letters. Some 
of the early crucial developments in the held were the invention of the laser in 1958 
[11, 12]; the diode laser in the 1962 [13]; the hrst demonstration of optical communica¬ 
tion using optical hbres in 1966, done at the University of Ulm, Germany [14, 15]; and 
the invention of the optical hbre ampliher in the early 1980’s [16]. Photonics presents 
some signihcant advantages over electronics. The most important one is speed. Due 
to its much larger carrier frequencies, photonics has a much larger bandwidth (or data 
transmission rate) than electronics. This is depicted in Figure 1.2, which shows that 
the operating speed of photonic devices can be up to hve orders of magnitude higher 
than electronics. Another important feature of photons is their weak interaction with 
the environment. Photons can travel long distances almost without changing their 
state. This implies much lower losses in transmission. Actually, this feature is also 
exploited in quantum information, in particular in the held of quantum cryptography 
[17] to create secure protocols to data exchange. As a consequence, photonic integrated 
circuits have been gaining a lot of attention in recent times [18]. The most ambitious 
goal is to eventually usurp the predominant role of electronic integrated circuits in in¬ 
formation devices. However, the goal seems to be unrealistic at the moment. Photonic 


^http: //www.intel.com/content / www/us/en/silicon-innovations/intel-22nm-technology.html 




4 


Introduction 



—I-1-1-1-1-1— 

lOnm lOOnm Ijim lOum 100|.im 1mm 



Critical dimension active devices (nm) 




Figure 1.2: Transmission data rates of different technologies and critical dimension of 
their active devices. Dielectric photonics refers to fibre optics technologies. It can be seen 
that photonics is much faster than electronics. 


circuits suffer from a lack of storage capacity, i.e. it is difficult to stop and store light. 
Hence, the storage needs to be done in an electronic medium. Moreover, integrated 
electronics is non-linear, whereas non-linear integrated photonics still needs to evolve 
to be competitive in the market. Last but not least, the typical size of photonic com¬ 
ponents is limited by the wavelength of light. Due to the diffraction limit of light, 
dielectric photonic circuit cannot resolve or retrieve information smaller than a frac¬ 
tion of the wavelength (200nm approximately) [19, 20], as opposed to nano-electronics 
where 22nm transistors are already a reality. 

Most of the efforts to overcome the diffraction limitations in photonics have been 
happening in the held of plasmonics. In fact, lots of authors consider that nanopho¬ 
tonics should be based on metallic nano-plasmonics (see Figure 1.3) [21]. Plasmonics 
is the science that studies the interaction between light and free electrons on a metal. 
The hrst theoretical studies in plasmonics were done in the 1950’s [22-25], and the 
hrst experiments with metallic surfaces in the 1970’s [26, 27]. In general, the strategy 
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Figure 1.3: Transmission data rates of different technologies and critical dimension of 
their active devices. Metallic nano-plasmonics is expected to be the technology that will 
bring photonics to the nano-scale. 


followed in plasmonics to bring photonics down to the nano-scale has a sample-based 
perspective. That is, given a hxed external excitation beam (see inset in Figure 1.4(a) 
for a typical example), samples are engineered (shape and materials) so that they pro¬ 
duce certain effects. For example, in Figure 1.4(a), an array of metallic nano-rods is 
fabricated to guide light from one location to another. Similarly, in Figure 1.4(b) a 
circular nano-hole drilled in a gold him deposited on top of a glass layer is inserted 
into a liquid chamber. The excitation of the structure enabled the authors of [28] to 
trap 50nm particles that had been dissolved in the water solution. There is no doubt 
that fabrication and characterization of samples play a huge role in the current plas- 
monic technologies. Interestingly, the spatial prohle of the incident beam also plays 
a very important role. However, most of the work in the held is done with Gaussian 
beams (the beam on the inset of Figure 1.4). In this thesis, it will be shown that 
more elaborate beams of light can be used to retrieve plenty of additional information 
from nano-structures. An illustrative example is the Stimulated Emission Depletion 
(STED) microscopy. STED microscopy was invented by Stefan Hell and co-workers 
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Figure 1.4: a) Plasmonic wave-guide made of metallic nano-rods. An incident Gaussian 
beam shone upon the structnre is carried along the plasmonic nano-rods. The inset is an 
intensity plot of a transverse section of the excitation beam. The image has been taken from 
activeplasmonics.org with permission of Anatoly Zayats. Credits to Dr. Ryan McCarron, 
King’s College, London, b) A metallic film with a nano-hole is excited with a Gaussian beam 
to excite a surface plasmon. The surface plasmon is used to trap a polystyrene particle. The 
image has been taken from [28], with permission of Mathieu L. Juan. 


in 1994 [29]. It is one of the so-called super-resolution microscopy methods, as it can 
resolve defects as tiny as 30nm. A very complete schematic of its working principle is 
shown in Figure 1.5, and its working principle can be found in [29, 30]. Even though 
Figure 1.5 is complicated, the message should be clear. Probing a sub-wavelength 
specimen with a Gaussian beam (green beam in Figure 1.5) results in a blurry image 
(Figure 1.5(c)). Nevertheless, the combined use of a Gaussian and a doughnut-shaped 
beam (red) results in a much neater image (see Figure 1.5(d)). There are different kinds 
of doughnut-shaped beams, but the ones used in STED are called vortex beams [31-33]. 
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Figure 1.5: Schematics of STED microscopy, a) A depletion laser beam (red) is added 
to a typical microscopy set-up. b) The use of a vortex beam considerably reduces the Point 
Spread Function of the system. The Point Spread Function describes the response of an 
imaging system to a point source. The largest the point spread function, the blurrier the 
image will be. c) Standard optical microscopy image of a sub-wavelength specimen, d) STED 
image of the same sub-wavelength specimen. Image from www.zeiss.com/campus. Copyright: 
Mike Davidson, FSU, Tallahassee. 


Vortex beams are defined by their topological charge Z, which is an integer nnmber. 
The topological charge I acconnts for the nnmber of times that the phase of the beam 
wraps aronnd its centre in a 27r circle. Examples of fonr different vortex beams are 
shown in Fignre 1.6. It can be seen that when I = —1, the phase goes from 0 to 2n 
one time in connter-clockwise direction. In contrast, when I = 3, the phase goes three 
times from 0 to 2?! in a clockwise direction. In chapter 6, sections 6.2 and 6.3 have been 
dedicated to explain how to characterize and generate vortex beams in the laboratory. 
Note that the dehnition of vortex beam and topological charge is independent of the 
polarisation. However, for many applications in nano-photonics (STED microscopy is 
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Figure 1.6: Intensity and phase profiles of vortex beams with topological charges I = 
— 1,0,1,3. For the intensity plots, red indicates maximum and dark blue minimum. For the 
phase, white means 0 phase, and black means 27r. 


one of them), vortex beams need to be tightly focnsed. When that is the case, vortex 
beams become much more complex. In fact, in that regime, the dehnition of the phase 
singularity is polarisation-dependent. This is what many authors have described as 
one of the manifestations of the spin-to-orbit conversion [34-38]. Then, characteris¬ 
ing vortex beams in terms of their symmetries becomes enormously useful. It allows 
for a systematic study of the beam without the need of considering different focusing 
regimes^. 


Sections 2.3 and 2.4 are dedicated to a systematic characterisation of beams in terms of 


^For the advanced reader: note the parallelism between the two focusing regimes (paraxial/non- 
paraxial) in optics with the two regimes in quantum mechanics (non-relativistic/relativistic). In the 
paraxial approximation, light can be described by a scalar function and polarization is added as a 
different degree of freedom. Similarly, the scalar wave-function of a particle is obtained with Schodinger 
equation, and the vectorial behaviour given by the spin is added by hand. In contrast, when light is 
focused a lot, Maxwell equations have to be used, and the polarization and spatial degrees of freedom 
are intermingled. Analogously, when a quantum particle is boosted and its speed is relativistic, Dirac 
equations must be used. In Dirac equations, the spin degrees of freedom cannot be modified without 
modifying the spatial properties of the field. 
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their symmetries. Thanks to the understanding and analysis of the different transfor¬ 
mations under which light beams are symmetric, different results have been obtained, 
both theoretically and experimentally. In the succeeding chapters, three symmetries 
have been specially studied: mirror symmetry, rotational symmetry, and duality [39- 
41]. Both samples and beams of light which are (or are not) symmetric under these 
transformations have been used and as a consequence new fundamental phenomena 
have been found^. In chapter 4, the interaction between single spheres and cylindri- 
cally symmetric beams has been theoretically studied. Among others, a technique to 
excite Whispering Gallery Modes (WGMs) has been developed. In chapter 5, the cylin¬ 
drical symmetry of a single sphere and the duality symmetry of a light beam have been 
used to cancel the forward scattering of a sample. Furthermore, it has been proven that 
a sample can effectively be dual in a controlled manner if its dimensions are changed 
in a certain way with respect to an incident cylindrically symmetric beam. In chap¬ 
ter 7, a mirror and cylindrically symmetric sample has been used to induce a giant 
Gircular Dichroism (GD). Again, the use of dual, cylindrically and mirror symmetric 
modes has been a key point to be able to achieve it. Finally, chapter 8 merges all the 
knowledge accumulated in the previous chapters, both theoretical and experimental. 
Using various dual cylindrically symmetric beams, the resonant scattering behaviour of 
a dielectric sphere is unveiled. This resonant behaviour is hidden under the Gaussian 
beam (or plane wave) excitation, and the use of vortex beams is crucial to find it. Due 
to time constraints, this last study is not as thorough as the others presented in the 
previous chapters. Nevertheless, it allowed for an experimental observation of some of 
the theoretical predictions of the thesis. 

Now, even though the details of what has been done in each chapter may be too 
technical for a non-specialised reader, these should only be considered to emphasise 
the point made before: nanophotonics can also be done from another perspective dif¬ 
ferent from the prevailing sample-based view. In particular, the consideration of the 

^Saying that a system is symmetric under a transformation means that when that transformation 
is applied to the system, its physical properties remain unchanged. 
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sample and the exciting light as a whole system, along with symmetry considerations, 
allows for the discovery of new phenomena at the nano-scale. The futnre in this nnex- 
plored direction looks as exciting as promising. In the next chapters, the reader will 
hnd a taste of it. 



“If I could explain it to the average person, I wouldn't have been 
worth the Nobel Prize. ” 


Richard P. Feynman 



Theoretical methods in nano-optics 


2.1 Maxwell Equations for a linear, non-dispersive, 
homogeneous, isotropic and source-free medium 


The material presented in this thesis can be described using the tools of EM field 
theory [42-44], In particular, the subdivisions of EM field theory that this thesis 
especially deals with are electromagnetic optics and vectorial scattering theory [45- 
52], In this chapter, I will introduce all the necessary concepts and tools to understand 
the theoretical developments scattered across the next chapters. For completeness, the 
starting point will be the macroscopic Maxwell Equations for a source-free medium 
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[42, 45]: 

V-T> = 0 Vx£ = 

V-i3 = 0 V xK = 

where £ and "H are the electric and magnetic fields, V is the electric displacement, and 
B is the magnetic indnction. The relations between V and £, as well as between B and 
Ti are given by the properties of the medinm: 


€()£ + V 

(2.2) 

-I- 

(2.3) 


with eo,/io being the electric permittivity and magnetic permeability of vacnnm. V is 
defined as the polarization density, which is the density of molecular dipole moments 
induced by the electric field [42, 45]. Similarly, A4, defined as the magnetization density, 
is the density of molecular magnetic dipole moments induced by the magnetic field. 
I will be particularly interested in linear, non-dispersive, homogeneous^ and isotropic 
media. In these cases, a linear relation can be established between V and £, and 
therefore the following relation holds [45] : 

V = e£ (2.4) 

where e is a scalar quantity defining the electrical properties of the medium called 
electric permittivity. In the same way, a linear relation can be cast between Ai and 
Ti, yielding the following linear relation between B and H 


dB 

dV 


( 2 , 1 ) 


B = iM 


(2.5) 


with /i being the so-called magnetic permeability constant of the medium. Hence, the 
macroscopic Maxwell Equations for a linear, non-dispersive, homogeneous, isotropic 
and source-free medium can be written only as functions of £ and H [45]: 


V£ = 0 

v-n = 0 


V X ^ 

VxH 




( 2 . 6 ) 


may consider also linear, non-dispersive, inhomogeneous media formed by different linear, non- 


dispersive, homogeneous, and isotropic submedia. 
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Here, I will be especially interested in the monochromatic solutions of these equations. 
That is, I will suppose that all the helds have a time dependence, with uj being 
the frequency of light. Then, any held can be expressed as the real part of a product 
of a purely temporal part and a purely spatial part: A{r,t) = Re [A(r)e“*‘^*], where 
A(r) represents a complex-amplitude vector. Then, equations (2.6) can be re-written 
just in terms of their spatially dependent parts: 


V-E = 0 VxE= zca/iH 
V-H = 0 VxH = —icaeE 


(2.7) 


Now, applying the curl operation to equations (2.7) and making use of the vector 
identity V x (V x A) = V(V ■ A) - V^A, one obtains that both E and H fulhl the 
vectorial Helmholtz Equation [42]: 


V^E + k‘^E = 0 
V^H + fc^H = 0 


( 2 . 8 ) 


where k = is known as the wave-number. Nonetheless, equations (2.7, 2.8) 

cannot be uniquely solved unless some boundary conditions (BC) are applied. Firstly, 
the EM helds must have a proper decay at inhnity and be hnite everywhere else. 
Secondly, when equations (2.7) are applied between two diherent media, the following 
BC need to hold: 


(e2E2 - eiEi) ■ fi = cr (/i2H2 - /iiHi) ■ n 

n X (E2 — El) = 0 n X (H2 — Hi) 

where a is the surface charge density at the surface; js is the surface current density; 
n is the unitary normal vector to the interface between 1 and 2; and {Ei,Hi} and 
{E2,H2} account for the electric and magnetic helds at two diherent media. It is 
observed that the tangential components of E and H must be continuous, while the 
normal components of D and B can be discontinuous depending on the properties of 
the surface. Equations (2.7) plus the BC listed above will be the framework in which 
this thesis will be grounded on. 


= 0 

js 


(2.9) 
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2.2 Fields, operators and symmetry groups 

In the previous section, I have introduced the monochromatic Maxwell Equations for 
linear, non-dispersive, homogeneous, isotropic and source-free media (see equation 

(2.7) ). I will impose that all the theoretical developments done in this thesis fulhl 
equations (2.7). That is, I will require that any meaningful EM field fulhls equations 

(2.7) , BCs (2.9), it decays properly at infinity and it is finite everywhere else. 

Another framework to do optical physics is the paraxial equation [45, 53], which is 
an approximation of equations (2.7) but it is much simpler to use. In paraxial optics, a 
complex-amplitude electric held E can be expressed as E(r) = i7(r)u, with E{r) being 
a complex-wave amplitude scalar function that must fulhl the paraxial equation 

{dl + d^ + 2ikd,)E{r) = 0 (2.10) 

and u being an arbitrary polarization vector. Then, it is clear that in paraxial optics 
the vectorial character of the held (given by u) is decoupled from its spatial properties 
(given by E{r)). The importance of the paraxial equation is that collimated beams, 
which are very frequently used in the laboratory, can be accurately described within 
this framework. In those beams, if the propagation direction of the beam is the 
axis, a plane wave decomposition of the wave yields partial waves whose kp -C fcz, and 
kz ^ k (see 2.3.1 to check the dehnitions of kp and kz)- This implies that the complex 
amplitude function E{v) can be expressed as E{y) = A(r) exp(i/c 2 ;), where y4(r) is a 
magnitude called the envelope of the wave which varies slowly with z: 

BA 

-^<^kA ( 2 . 11 ) 

dz 

However, the paraxial approximation does not properly work in nano-optics. In nano¬ 
optics, light beams are either tightly focused or interact with sub-wavelength objects. 
In those cases, the polarization degrees of freedom of the beams cannot be decoupled 
from their spatial properties [49, 54]. Then, imposing that all EM helds must be so¬ 
lutions of the full Maxwell equations rules out all those beams that fulfil the paraxial 
equation, but are not solutions of the full Maxwell equations given by equation (2.7) 
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[53, 55, 56], 

Now, I will define an operator as an entity that picks up any EM field fulfilling 
equations 2.7 and transforms it into a different one that still fulfils equations 2.7. Fi¬ 
nally, I will define a symmetry group as a set of operators with group structure. A 
group G is defined as a set of elements {a,6, c,...} along with an inner product. The 
inner product (■) is such that the product of two arbitrary elements of the group yields 
an element of the group [39]. Furthermore, the product must satisfy the following 
conditions [39]: 

• It must be associative, i.e. a ■ {h ■ c) = {a ■ b) ■ c, Va, b,c & G. 

• There must be an element e G G, called the identity, which is such that a - e = a, 
Va G G. 

• For each a E G, there must exist an element a~^ G G, called the inverse of a, 
such that a ■ a~^ = e. 

Groups can be subdivided into two categories - discrete, and continuous. In physics, 
most of continuous groups of interest are part of a subdivision called ‘linear Lie groups’ 
[39]. Linear Lie groups can be labelled by one or more continuous variables. In fact, 
each label is related to a generator. The generators of the group define the local be¬ 
haviour of the group near the identity, and any element of the group can be obtained as 
a smooth function of the generator. For example, for the S'0(2) group of 2-dimensional 
rotations, any rotation R{(j)) can be obtained as R{(j)) = exp(—f0J), with J the gen¬ 
erator of the group. Note that this implies that if the element of the group is an 
operator, so must be the generator. In fact, group generators are very good candidates 
to compose a complete set of commuting operators. A complete set of commuting op¬ 
erators of a vector space is a set of operators that commute with each other and whose 
eigenvalues completely specify the state of a system [57, 58]. The idea is that if an EM 
field is an eigenvector of a generator of a symmetry, then it is also invariant under the 
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symmetry transformation generated by that generator, and vice versa^ [39, 59]. Then, 
nsing generators to make a set of commnting operators will be nsefnl as the helds will 
be characterized in terms of their symmetries. Now, I will differentiate between EM 
helds and EM modes. EM helds are solntions of Maxwell eqnations (2.7), whereas 
EM Modes are elements of a basis of solntions of Maxwell eqnations (2.7). Thns, in 
general, any EM held can be decomposed into EM modes. In addition, EM modes will 
be chosen so that they are eigenstates of some generators, i.e. they will be invariant 
nnder certain symmetries. Next, I will list the most common operators that I will nse 
to describe EM modes and their diherential representation will be given. Some of them 
are generators of symmetries, so in those cases I will also comment on the symmetries 
that they give rise to. 

d 

i) Hamiltonian , H = The hamiltonian is the generator of temporal transla¬ 

tions (or time evolntion) [40]. That is, the transformation Ta* = exp {—iHAt) 
applied to an EM held -(/^(r, t) the state evolve from t = Iq to t' = Iq + At: 

Ta* [i/’(r,f)] = V^(r,fAt) (2.12) 

ii) Linear momentnm operator, P = —iV. This is a vectorial operator, i.e. it is 
composed of three operators in three diherent directions: P = PxX+Pyj+PzZ. Its 
three components commnte with each other [Pj, Pj] = 0. Each of the components 
generates a translation in its corresponding axis. Hence, the vectorial addition 
of all of them is the generator of linear translations in any direction of space: 
PAr = exp (—iP ■ Ar) [39]. That is, 

^Ar [V’(r,t)] = ^(r + Ar,f) (2.13) 

iii) Angnlar momentnm operator, J = L + S, with L = —i(r x V) and S = 

— i (cinmX -|- e 2 nmy + ^ 3 nm^) respectively, where einm is the total anti-symmetric 
tensor and l,m,n = 1,2,3. It is also composed of three operators in three 

^This statement is a consequence of Noether’s theorem, which states that every continuous sym¬ 
metry of a physical system corresponds to a conservation law of the system. 
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components, J = Jx^ + Jyy + but they do not commute with each other: 
[Jj, Jj\ = ieijkJk- The reason why the three components of the angular momen¬ 
tum (AM) J do not commute is that they are the generators of rotations on three 
different axis, and rotations do not commute (whereas spatial translations do, see 
above) [51]. A general rotation of an angle ip along an axis on the direction n 
is given by: Rn{ip) = exp {—iJ ■ inp). Then, the application of a rotation onto a 
general EM held yields: 

RM ['0(r, t)] = Mfi((p) • '0(Mri((p) ■ r, t) (2.14) 


where Mn((p) is a rotation matrix [39] (see Appendix A). Finally, note that even 
though the L and S operators used to dehne J are usually called orbital and 
spin AM respectively, they are not proper AM momentum operator on their own 
[54, 60-62]. None of them generates proper rotations, only their addition does. 
Furthermore, they are not proper operators, as their application to an EM held 
fulhlling Maxwell equations yields a held that does not fulhl Maxwell equations 
[54, 60-62]. 


iv) Angular momentum squared, + Jy + J^, where R have been dehned in 

iii). is called a Casimir operator, and commutes with all the operators. Even 
though it does not generate any meaningful transformation, it is very useful for 
problems with spherical symmetry. 


v) Helicity operator, A = 


J P 


The helicity operator is dehned as the projection 


of the AM onto the direction of the linear momentum. It was proven in [63, 64] 
that A generates generalized duality transformations [41, 42]: = exp {—iRip). 

Generalized duality transformations mix the electric and magnetic parts of an 
EM held. That is, if we express RirR) as V’(r,f) = (E(r, f), H(r, f))"^, then 


[^(r,t)] = 


E(r,t) 

H(r,t) 


E(r, t) cos ip — H(r, t) sin ip 
E(r, t) sin ip + H(r, t) cos ip 


(2.15) 


Its diherential representation for monochromatic helds is A 


Vx 


Unlike all the 


k 
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rest of symmetries listed hereby, duality symmetry is a non-geometrical symme¬ 
try. It can be proven that its preservation only depends on the material properties 
of the medium in consideration [41]. A very thorough study of duality symmetry 
and its generator (helicity) can be found in [65]. 

vi) Parity operator, If. Parity transformations (r —>■ — r) are discrete, therefore they 


are not generated by any other operator. The application of If onto a Maxwell 
held yields [40]: 




(2.16) 


vii) Mirror symmetry operator, A mirror transformation of the kind indi¬ 

cates that the spatial coordinates are rehected off a plane given by its normal 
vector h^. Actually, it can be re-written as = ni?n(7r), i.e. a mirror sym¬ 
metry operator can be obtained by multiplying a parity and rotation operators. 
Like parity, it is a discrete transformation, so it is not generated by any other 
operator. Its application onto gives [40]: 


Mn f)] = V’(i‘ ~ 2n(h ■ r),t) 


(2.17) 


Last but not least, I will give the most useful commutation rules between the operators 
listed above [40]. These commutation rules are crucial in order to dehne a set of 
commuting operators. In fact, it is known that only four commuting operators are 
needed to describe EM helds, as they are elements of the Poincare group [39, 66]: 

• T^t commutes with all the other operators of the list above. 

• P commutes with Taj, J^, A, and Mn. It does not commute with J, their 
commutator is [Jj, Pj] = ietjkPk- It anti-commutes with parity {Pi, 11} = 0. 

• J commutes with Taj, J^, A, If and M^ (when h is parallel to J). It does not 
commute with M{n}, in fact it anti-commutes: JjM{n} = —M{ni,Jj. 

^Throughout this thesis, another notation will also be used for mirror transformations: . The 

difference is that in the reflection is off a plane whose normal vector is n, whereas Mjflj. indicates 
that the reflection is done upon a plane that contains the axis n. 
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• commutes with all the rest of operators, as it is a Casimir operator. 


• A commutes with Taj, P, J and J^. It anti-commutes both with If and Mn. 


• n commutes with T^t, and M^. It anti-commutes with A and P. 


2.3 Basis of solutions of free-source Maxwell equa¬ 
tions 

The aim of this chapter is to establish the basis on which the works presented on 
this thesis are grounded. In section 2.1, I have introduced Maxwell equations. Then, 
in section 2.2, I have discussed that I will require all the EM helds to fulhl Maxwell 
equations in its form given by equations (2.7). Furthermore, some operators and gen¬ 
erators have been introduced. In this section and the next one, six different sets of 
EM modes that stem from six different sets of complete commuting operators will be 
described. Each of these sets of EM modes will be especially suited to describe some 
particular EM interactions. In fact, since generators of symmetries will be used as 
the operators in a set, a set of modes will be chosen over another depending on the 
symmetries of the interaction. There are different ways of constructing these sets of 
EM modes. From a symmetries’ point of view, probably the clearest technique is ‘the 
induced representation method’. Wu-Ki Tung describes it in chapter 9 of his book 
[39]. As in most of the cases, the use of symmetries simplifies the amount of opera¬ 
tions tremendously. Unfortunately, this method is not very common in optics. Hence, 
another equivalent method much better-known in the optics community will be used 
here. This method is based on some mathematical theorems about the solutions of 
the vectorial Helmholtz equation (2.8). The method is described in [43, 46, 47] among 
others. It is based on the following property. Given a general solution of the scalar 
Helmholtz equation in a system of coordinates where the solution is in separated vari¬ 
ables 4/(xi, 0 : 2 , 2 : 3 ) = Xi{xi)X 2 {x 2 )X^{x^), then it can be proven that if a is a unitary 
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vector in a privileged direction^, 

M™ = V X (aT(Ti,a;2,a;3)) 

yJi. -1-0j 

M™ = {(Vx)/fc}[V X (aV>(a;i,X2,a;3))] 

are general solutions of the vectorial Helmholtz equation and fulhl the divergence equa¬ 
tions V ■ E = V ■ H = 0 [43, 47] . That is, any E field can be obtained as a superposition 
of the modes and dehned in equation (2.18). Then, the magnetic field H 

can be derived from Maxwell equations (2.7)^. Note that I have differentiated the 
behaviour of the two general solutions by its TE or TM character. The notation of 
TE/TM behaviour is a bit arbitrary, as both the electric and magnetic helds are ‘trans¬ 
verse’ in the Maxwell sense, i.e. V-E = V-H = 0. In fact, this notation gets even more 
cumbersome when spherical coordinates are used, as TE modes are usually referred to 
as ‘magnetic’ and TM as ’electric’. Here, I will stick to the most common notation 
and denote the hrst solution in equation (2.18) as TE and the second one as TM for 
cartesian and cylindrical coordinates [43, 47, 48]. Then, in spherical coordinates I will 
refer to the hrst solution as magnetic and the second one as electric [42, 50, 51], as 
that is the dominant notation. 


2.3.1 Plane waves 


It can be shown that T(a;, y, z) = exp(ik-r) is a general solution of the scalar Helmholtz 
equation for cartesian coordinates [43], with r = (t, y, z) being the position vector and 
k = {kx,ky,kz) the wave-vector, whose modulus |k| = k = co/c is the wave-number. 
Then, using a = z as the symmetry direction, the two following vectorial solutions 
arise: 


M™ = ?exp( 2 k-r)= — (fc^x — fc^jy) exp(ik ■ r) 

fhp 


= pexp(*k-r) = 


-kz {kx^ + kyy) + k^z 
kkn 


exp( 2 k ■ r) 


(2.19) 

( 2 . 20 ) 


^Not any vector a is valid. Some conditions apply. These conditions can be found in [47]. 

®As H can always be derived from E, I will only focus on the properties of E, and derive H using 

equations (2.7). Besides, any symmetry considerations regarding E will also apply to H, except for 
parity and mirror inversion. It can be shown that when E is odd, H is even and viceversa [42] 
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where {x, y,z} are the polarization vectors in cartesian directions. These TE and TM 
modes are very well-known in optics and are referred to as s' and p plane waves. In fact, 
they dehne a trihedron such that p x s' = k/|k|. As previously mentioned, all the helds 
considered here are monochromatic, therefore both s and p waves will be eigenvectors 
of the H operator. In addition, it can be checked that they are eigenvectors of the 
following operators: 


Fa; ['s'exp(*k ■ r)] = 's'exp(*k ■ r) F^, [pexp(ik • r)] = fca, pexp( 2 k ■ r) 

Fy ['s'exp(ik • r)] = 's'exp(*k • r) F^ [pexp(*k ■ r)] = pexp(ik ■ r) 

['s'exp(*k ■ r)] = —'sexp(ik • r) [pexp(ik ■ r)] = -|-pexp(ik ■ r) 

( 2 . 21 ) 

Then, ^ and p are EM modes invariant under temporal translations, spatial translations 
along X and y axis, and mirror symmetry transformations. Hence, they are eigenvectors 
of F, Fa;, Py, and M-§. Now, if we look at equations (2.18), we can see that the M™ 
mode (p wave) is obtained by applying (Vx)/fc to the M™ mode ('s wave). In section 
2.2, I have shown that the operation (Vx)/fc is actually the differential form of the 
helicity operator for monochromatic helds: 


A 


J P 
P 


Vx 

~k~ 


Then, it can immediately be seen that: 


( 2 . 22 ) 


A ['s'exp(ik • r)] = pexp(ik-r) (2.23) 

A[p exp(ik ■ r)] = s'exp(ik-r) (2.24) 


Thus, additions and subtractions of s' and p waves yield states with a well-dehned 
helicity®: 


A 

A 


s -h p 

.~7T 

s-p 

. V2 


exp(ik ■ r) 
exp(ik ■ r) 


= + 


s -h p 

.~VT 

s-p 

. 72 


exp(*k • r) 
exp( 2 k ■ r) 


(2.25) 

(2.26) 


® Saying that a state has a well-defined property (helicity, for example) means that the state is an 
eigenstate of the operator representing that property. 
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I will denote these combinations as = (s'-|-pp)/\/2, with p = ±1 depending on the 
value of helicity: 
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(2.27) 

(2.28) 


where kp and (f)k are the wave-vector components in cylindrical coordinates, i.e. k^ = 
kpCos(j)k, ky = /cpsin0fc and {&+,&-} are the circular polarizations vectors: (j+ = 
(x -|- iy)/\/2 is the left circular polarization and d- = (x — iy)/\/2 is the right one. 
Then, a new basis of EM modes can be created with e_|_ and e_. This will be symmetric 
under temporal translations, spatial translations along x and y axis, and generalized 
duality transformations: 

Px [e+ exp(*k ■ r)] = kx e+ exp(ik ■ r) Px [e_ exp(ik • r)] = kx e_ exp(ik ■ r) 

Py [e+ exp(ik ■ r)] = ky e+ exp(ik • r) Py [e_ exp( 2 k ■ r)] = ky e_ exp(*k • r) 

A [e+exp(ik ■ r)] = -|-e+exp( 2 k ■ r) A [e_ exp(ik ■ r)] = — e_exp(ik-r) 

(2.29) 

Finally, it can be checked that the two sets of EM modes given here are normalized in 
the following way: 


('sexp(*k • r))* ■ s'exp(ik'■ r)d^r = 5(k — k') 

(pexp(*k ■ r))* ■ pexp(ik'■ r)d^r = 5(k - k') 

(bp exp(*k ■ r))* ■ e^/ exp(ik' ■ r)d^r = 5(k — k') 5,^ 


(2.30) 


where (5(k — k') is a Dirac delta, and 6ppi is a Kronecker delta [47]. 


2.3.2 Bessel beams 

Cylindrically symmetric interactions will play a very signihcant role in this thesis. 
Therefore, a basis that exploits this symmetry needs to be found. A way of attacking 
the problem is to use equations (2.18) in cylindrical coordinates with a = z as well. 
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This method, used by [43, 47, 56, 67] and several others, naturally gives rise to the 
so-called TE/TM Bessel beams. Nevertheless, the symmetries of the problem are a bit 
blurred by the methodology. Then, in order to highlight the symmetries of the Bessel 
beams, I will use an alternative method to derive them. This method is based on 
the fact that any field can be decomposed into plane waves. Then, the superpositions 
are chosen so that the symmetry requirements are fulfilled. A Bessel beam can be 
defined as an EM mode symmetric under temporal translations, spatial translations 
and rotations along a symmetry axis (usually, the 2 ; is chosen). This implies that 
Bessel beams need to be eigenstates of H, P^, and Jz- As mentioned earlier, since I am 
only considering monochromatic fields, they are immediately eigenstates of H. Then, 
having a set u fixes a relation between kz and kp\ = k'^ = k^ + k^. Therefore, 

the z component and the transverse component of the linear momentum are linked via 
the frequency. Thus, having a well-defined Pz is analogous of having a well-defined Pp. 
Consequently, adding plane waves on a circular cone symmetric under the z axis will 
yield a field with a well-defined Pp. Then, I will be able to find the conditions on the 
superposition coefficients by imposing that the beam is an eigenvector of Jz- To start, 
let me re-write the four plane waves derived in the previous sub-section in terms of 
rotations (see Appendix A): 


s'exp( 2 k ■ r) 

= -K{(()k,Ok)iyexp{ikz) 

(2.31) 

pexp(*k ■ r) 

= -R(0fc,6'fc)xexp(ifc2;) 

(2.32) 

e+ exp(*k ■ r) 

= -R{(l)k,Ok)a+exp{ikz) 

(2.33) 

e_ exp(ik ■ r) 

= -R{(f)k, Ok)a_ exp{ikz) 

(2.34) 


where R(0fc,0fc) is the operator of rotations defined in Appendix A. Then, due to the 
symmetry considerations given by equations (2.21, 2.29), a superposition of s' waves 
will yield a mode with a well-defined frequency and a defined TE/TM character. In 
general, it will no longer be an eigenstate of P^ and Py, but if we only add waves 
on a cone, as mentioned in the paragraph above, we can retrieve a mode with a well 
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defined Pz (or Pp). This can be given by the general snperposition; 


- 27 r 
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sm6kd6k d(j)ka{<Pk,6k)sexp{ik-r) 


r- 27 r 


(2.35) 


= - sin6kd6k / d4>k<y{4>k,0k)'R-{4>k,0k)iy 

Jo Jo 


where a(0fc, 6k) = f{(j)k)S{dk — 6'^). By constrnction, the field B™ is an eigenstate of H 
and Pp. In order to make it a Bessel beam, a right fnnction /(0fc) needs to be chosen 
so that B™ is an eigenstate of It can be proven that /(0fe) = ^exp(im0fc) does 
exactly this, and it also assnres the nnitarity of the transformation: 
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S{dk - 6'j^)R{(j)k, 6k)y exp{ikz) = mBj 


TE 
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(2.36) 


where an integration by parts has been carried ont on the last step. It is straightfor¬ 
ward to see that similar Bessel beams conld have been obtained if instead of rotating 
iy exp{ikz), the other three plane waves listed on eqnations (2.32-2.34) had been nsed. 
This fact allows for the dehnition of two sets of EM modes. I will denote the hrst set 
as B^^^, with (y) = TE/TM, or Bessel beams^. Besides H, they are eigenstates of the 
following operators: 
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*mkf, 

(y) 

mkn 


= kr, B, 
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= m B 
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(2.37) 


mkn 


Their expressions as fnnctions of ^ and p waves are: 


r- 27 r 
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^mkn 


= — sin^fcd^fc / d(j)kS{9k - 9'^)exp{im(j)k)sexp{ik-r) 
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r- 27 r 
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d(l)kS{9k — 9'p) exp(im0fc)p exp(ik ■ r) 


(2.38) 


^Due to the fact that these are the Bessel beams that the community generally uses, I denote them 
as Bessel beams, instead of Bessel beams with a TE/TM character. 
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Whilst their expressions in the real space are: 
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(2.39) 


where kp = k sin 9[ and the following eqnality has been used [43] : 

exp(ik ■ r) = i^Jn{kpp) exp(m(0 — 0^)) Qyip{ikzz) (2.40) 

n 

with Jn{kpp) being a Bessel function of the first kind [68, 69]. An analogous procedure 
can be followed to derive the Bessel beams with a well-dehned helicity, Their 

expressions as functions of the waves can be compacted into: 


r-27r 


Bmfcp = ^ / sm9kd9k I d(j)k5{9k-9'f,)exp{im(j)k)epexp{ik-ir) (2.41) 


Which yields the following expression for the real space helds: 
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(2.42) 


As for the symmetries of B!* . , they are eigenstates of hf, and the following operators: 
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It is clear, then, that the relation between the parity and helicity modes is the following 
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As with the plane waves, the two sets of Bessel beams are also normalized: 


(B"„. )• ■ b;;:' d’r = s(k, - k'Mk, - 


(2.45) 


Finally, for completeness, the decomposition of plane waves into Bessel beams is also 
given: 


s'exp(ik ■ r) 
p exp(ik ■ r) 
Bp exp(ik ■ r) 


m=—oo 
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siiioyo' ^ 
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siliOidO' ^ 
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Ok) exp{-im(j)k)B'^^^ 
Ok) exp(-im0fc)B™ 
0fc)exp(-im0fc)B^fc^ 


(2.46) 

(2.47) 

(2.48) 


2.4 Multipolar fields 

The mnltipolar helds® are EM modes specially useful to describe EM interactions with 
spherical symmetry. They are widely used in antenna theory [70], Mie Theory (see 
sections (3-5)), nuclear [71], atomic and molecular physics [72-74], astrophysics [75], 
and many more fields. The fact that they are used in such vast variety of fields makes 
a unique notation almost impossible. Sometimes, even within the same field, different 
notations are used, making it hard to reconcile calculations from different sources [76]. 
My notation is very similar to Rose’s [50, 51], and it is based on the symmetries of the 
modes. I will denote a multipolar held as where {y) = (m), (e) accounts for the 

parity of the held®. The notation becomes clear when we look at the symmetries of the 

®Here, I will call multipolar fields to those multipolar fields with a well-defined parity. I will show 
that multipolar fields can also have a well-defined helicity, and those will be denoted as multipolar 
fields with a well-defined helicity. 

^Remember that a TE parity (or mode) is also called ‘magnetic’; while a TM parity can also be 


called ’electric’. 
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modes. Besides the hamiltonian if, these modes are also symmetric under: 



aM 

jmz 



a(-) 

= j(j + i)a;i 

Jz 

A M 

A (fP) 

= 

Jz 

a(-) 

jPPz 

A (^) 

n 

aM 

= (-IpA'S 

n 

A^) 

jrriz 

= (-ir«A;i 


(2.49) 


Like Bessel beams, they can also be expressed as general superpositions of plane waves. 
Their expressions can be found in [39, 61, 77]: 
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(2.50) 


where D{^^p{(j)k.,9k) is the Wigner rotation matrix [39, 51]. Note that the symme¬ 
try requirements dehned by equations (2.49) do not imply anything about the radial 
dependence of the multipoles. Nonetheless, in the plane-wave decomposition given 
by equations (2.50), the radial function has been chosen so that these multipoles are 
regular at the origin^®. Thus, they have the following expressions in the real space 
[46, 50, 76, 78]: 
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(2.51) 


where jjikr) is a spherical Bessel function [68, 69], P'^^irf) is a generalized Legendre 
polynomial [47, 68], rj = cos 9, |f,0, are the polarization vectors in spherical co¬ 
ordinates, and Cjm,, is a constant so that the modes are normalized. Its expression is 
given by: 


^2 ^ 47r(2j + l)(j + m^)! 

2k^j{j + l){j-m,)\ 


^°Notice that the radial dependence determines the behaviour of the multipolar field at infinity. In 
the next couple of pages, another radial function will be introduced. This radial function, known as 
Hankel function, is frequently used to model emission or scattering problems. 
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Indeed, with Cjm^ the multipolar helds are an orthonormal basis: 


/ (A‘1)* • = S(k - 

7r3 

Similarly, to the plane waves and Bessel beams, a new set of multipolar helds 
dehned by adding and subtracting the electric and magnetic modes: 






jmz jm^ 


y /2 ^2 

This implies that their plane wave decomposition is [39]: 
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which yields the following real space expression: 
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Besides of the H operator, the multipolar helds with well-dehned helicity A' 
eigenstate of the following three operators: 


— j(j + 1) 

[AjmJ = P 


Finally, by construction, A^^^ are also normalized: 



6(^k k )5jj!S^^Yn'^^pp' 


(2.53) 
can be 
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are 


(2.57) 


(2.58) 


Now, as it will be clear in chapter 3, there is a need to dehne other multipolar helds. 
Indeed, the two sets of multipolar helds dehned by equations (2.51, 2.56) have a far- 
held {kr —)• oo) behaviour which is not desired for scattering problems [42, 46, 48]. 
















2.4 Multipolar fields 


29 


This behaviour stems from the spherical Bessel function: 



(2.59) 


That is, the multipolar fields given by equations (2.51, 2.56) behave as stationary spher¬ 
ical waves in the far-held. Nevertheless, lots of EM problems with spherical symmetry 
have an emitting nature, e.g. radiation emitted from antennas [79], quantum dots 
[80, 81], nitrogen-vacancies [82], electrons [83], among others. One of these problems, 
which lays in the core of this thesis, is the Mie scattering problem, where a sphere em¬ 
bedded in a homogeneous medium scatters light that is illuminated on it [46, 48, 78, 84]. 
In all these problems, the spherical Bessel function of hrst kind must be replaced by 
a Hankel function of hrst kind h^p{krY^ [42, 69]. Because hp\kr) will not be used 
throughout this thesis, I denote the Hankel function of the hrst kind simply as hj{kr). 
In the far-held, the Hankel function produces a desired response, i.e. it behaves like 
an outgoing spherical wave: 



(2.60) 


The real expressions of these new multipolar helds is analogous to the ones presented 
by equations (2.51, 2.56), where hj{kr) replaces jj{kr)\ 



0 § 


(2.61) 


the harmonic dependence of all the fields in the problem had been chosen to be instead of 
the Hankel function of second kind hp\kr) should have been chosen [48]. 
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However, the plane wave decomposition of A^;^^ changes dramatically. 

The Hankel fnnction is singnlar in the origin, lim \hj{kr) \ —>■ oo, and conseqnently it 

kr^oo 

cannot be expanded only as a snperposition of propagating plane waves. Evanescent 
waves must appear in the superposition. Moreover, it is necessary to split the domain 
into two subspaces, z ^0. The decomposition is given by the following equation [85]: 
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where k is given by: 
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and the spherical momentum variables 0^ and 6k must be computed as functions of 
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(2.64) 


The expression for the semi-plane z < 0 can be obtained substituting the factor exp(ik- 
r)/{Kk) by exp(—zk ■ r)/{Kk) [85]. Equations (2.62) have large implications in nano¬ 
optics. Evanescent waves play a key role in obtaining information beyond the diffraction 
limit of light [49]. Here, it is shown that evanescent waves are an integral part of Hankel 
multipolar helds. Furthermore, it is seen that, like evanescent waves, Hankel multipolar 
helds are intrinsically linked to interfaces, as they are not well-dehned on the whole 
space. This fact is largely exploited in the excitation of WGMs (see section 4.4). As 
described by Oraevsky in [86] , given a sphere of radius R, Hankel multipolar helds have 
an evanescent dependence for a short range of r’s outside the spherical surface: 
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(2.65) 
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Last but not least, some authors have managed to link the evanescent behaviour of 
plane waves to complex rotations [87, 88]. In fact, this is what is seen in the expressions 
(2.62), where complex angles would be needed if the double integral f f dk;j;dky was to 
be expressed in spherical coordinates. 

2.5 Scattering control 

The work presented in this thesis belongs to the held of nano-optics. That is, most of 
the work deals with light-matter interactions where the light is at optical wavelengths 
(A G [390 — 780] nm [45]) and the matter is a nano-structure. I will be specially 
interested in single nano-structures, therefore it will be almost mandatory to tightly 
focus the beam onto the structure, so that the light-matter interaction is stronger. 
However, this is not the typical approach in nanophotonics or metamaterials, where 
structures are mostly designed to interact with paraxial Gaussian beams or plane waves 
[89-92]. In general, the approach followed in the design of nano-circuits or nano¬ 
materials is the following one. The nano-structure is characterized by its scattering 
matrix S. This scattering matrix is a function of many geometrical and material 
properties of the system, S(g, m), where g and m are general sets of variables describing 
the geometrical and material properties of the structure. However, the properties of 
the structure do not depend on the incident held. Then, the response of the structure 
to an incoming held E“ can be cast as a convolution of S(g, m) with E™; 

E°"*(r)= (S(g,m)*E“)(r) (2.66) 

Since the properties of the structure do not depend on E™ and the incoming held 
is well-known (a plane wave), the light-matter interaction is reduced to a complete 
characterization of the scattering matrix S(g,m). Then, adjusting the geometry or 
material of the structure, a controlled interaction can be carried out. In this thesis, 
I will look at these light-matter interactions from a diherent perspective [93]. I will 
study single highly symmetric nano-structures whose scattering matrix are analytical 
or easy to compute. Then, I will use diherent plane waves decompositions to control 
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the scattering. That is, instead of controlling with the geometry and materials of 
the structure, the interaction will be controlled with the incoming held E“: 

E°"*(r) = J d^k'i exp(ik ■ r) (S(g, m) * E“(k'i)) (r) (2.67) 

where k'; is each of the different plane waves that take part in the superposition, and 
E“(k'i) is the Fourier transform of the incident held, which modulates the plane wave 
decomposition. Thus, it will be crucial to know the analytical expression and the 
symmetries of the incoming beam. In addition, as previously mentioned, due to the 
fact that single structures will be used, tightly focused light beams will be used to 
strengthen the light-matter interaction. Tightly focused beams can be analytically 
described with the aplanatic model of a lens. This is explained in the following section. 


2.6 Aplanatic lens model 

In this section, I study one of the most frequently used models to describe the focusing 
of light beams, i.e. the aplanatic lens model [49, 52, 94-96]. In particular, I will be 
especially interested in the symmetries of this model. I will show that the aplanatic 
lens model does not change the helicity nor the AM content of the beam as long as 
some conditions are met. When those are met, the focusing of a Bessel beam with 
a well-dehned helicity such as can be easily modelled as = J . 

Thus, will still be an eigenvector of A and J^. But before getting there, the 

aplanatic model of a lens is explained. The variables playing a role in the model are 
the collimated incident beam E“, the focal distance of the lens /, and the refractive 
indexes of the media on the two sides of the lens n, rii. The geometrical representation 
of the aplanatic model is given by Figure 2.1. As it can be seen there, the collimated 
beam E“ travels in a medium whose refractive index is n and hits the back-aperture 
of a lens. The analytical description of the collimated beam at the a ^ plane of the 
back-aperture of the lens is done in the paraxial approximation, i.e. E“ = E^'^u with 
u being the polarisation vector. The incident beam gets refracted and its expression 
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n Ui 
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Figure 2.1: Schematics of the aplanatic lens model. E™ is the incident collimated electric 
field; is the focnsed field; n, ni are the refractive indexes of the media embedding the lens 
before and after the field is focnsed; |/3, and are polarization vectors in cylindrical 

and spherical coordinates. 


Eoo(0fc,6*fc) is obtained as: 


Eqo(0A:) ^k) 


r{ci>k,9k) ( e -. 


^ + tp{(Pk,ek) (E“-p)0 



( 2 . 68 ) 


where t^{(j)k,9k),t^{(pk,9k) are the effective Fresnel coefficients of the lens for s' and 
p waves respectively [49, 52]; |p, 0, z| and |f,0, are the polarization vectors in 
cylindrical and spherical coordinates (see Figure 2.1). Note that three assumptions have 
been made here. First, the sine condition of geometrical optics has been considered. 
The sine condition states that any ray emerging or converging to the focus of an 
aplanatic system meets its conjugate ray at the surface of a sphere of radius /, at a 
distance of the optical axis h = f sin 9, where 9 is the divergence angle [49]. Second, the 
intensity law has also been assumed. The intensity law implies that the energy incident 
on the aplanatic lens equals the energy leaving it. It can be expressed as = 

|E“|'\/n/r/(?7,i/ii) cos^/^ 9 [49]. Finally, it is imposed that the real space-expression of 
E™ on that z plane gives rise to the angular spectrum of the focused beam. That 
is, when introducing the expression of E“ in equation (2.68), the following change of 
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variables needs to be done: 


E“(p,0) = E“(/sin0fc>fc) 


(2.69) 


where no z dependence is snpposed as the expression is taken on a plane with a constant 
z. With all this, the electric held at the focns of the lens is obtained as: 



(2.70) 


with / being the focal distance of the lens, and 9^ = arcsin(NA) with NA the numerical 
aperture of the lens. 


2.6.1 A preservation 

The fact that the aplanatic model preserves the helicity of light is not trivial nor a 
casualty. In fact, simple lenses made of glass do not preserve helicity. However, it is 
interesting for different applications that the polarisation of light is preserved under 
a microscope set-up. In order to achieve that (which implies helicity preservation), 
microscope objective manufacturers need to apply a special coating to microscope 
objectives [49? ]. This fact is translated into a condition relating the effective Fresnel 
coefficients In order to show this condition, a general plane wave with a well- 

dehned helicity p will be focused down using the aplanatic lens model. It will be 
seen that the helicity of the beam is preserved provided a certain condition is met. 
Consider a general incident paraxial beam with a well-defined helicity p within the 
paraxial approximation and a exp(ik ■ r) dependence: 


E^" = E*A0A:,/sin0fc)dp 


(2.71) 


where cf)) is the amplitude of the electric held on the plane of the back-aperture 

of the lens. Then, the computation of Eoo(0fc,6*fc) yields: 



(2.72) 
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where the following equalities have been used: 


(2.73) 





(2.74) 


Equation (2.72) can be expressed in terms of s' and p; 



(2.75) 


where the two following relations have been used: 


p = -e 


s 


(2.76) 

(2.77) 


Hence, examining equation (2.75), it is straightforward to conclude that the aplanatic 
lens model preserves helicity if and only if 


Ok) = Ok) = Ok) W(j)k, Ok 


(2.78) 


I will consider that equation (2.78) is valid for the rest of my theoretical calculations. 

2.6.2 Jz preservation 

Looking at Figure 2.1, one notices that the geometrical representation of the aplanatic 
model is cylindrically symmetric around the axis. Thanks to Noether’s theorem, 
in principle it follows that the model preserves J^- Nevertheless, in the same way 
as it happened in the previous subsection, the Fresnel coefficients play an important 
role. Indeed, next it will be proven that either the Fresnel coefficients fulhl a certain 
condition, or the model is not cylindrically symmetric and therefore does not preserves 
Jz- Using equation (2.70), the expression for the electric held at the focus of the 
aplanatic lens when an incoming beam given by equation (2.71) can be computed: 



(2.79) 
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It is interesting to see that equation (2.79) can be expressed as a superposition of 
rotated plane waves with a well-defined helicity. The way to construct such a super¬ 
position is intuitive and reflects very well the definition of the helicity of a beam. We 
take a circularly polarized plane wave propagating along the z axis with helicity p. 
We rotate it in all directions and we sum all the different contributions with a certain 
weight g{4>k,(^k)- This rotation does not change the polarization in the system of ref¬ 
erence of the plane wave (it maintains transversality), therefore the helicity has still a 
well-defined value p. 

P7T p27T 

AP= sm{ek)d6k / d(j)k g{(j)k,Ok)R{(j)k,Ok)d-pexp{ikz) (2.80) 

Jo Jo 

where is a general held with a well-dehned helicity. Then, using equations (2.33, 
2.34), can be re-written as: 

p7Z p27V 

AP = - sm{9k)d9k d(j)k g{(l)k,0k)epexp{ik ■ v) (2.81) 

Jo Jo 

At this point, a direct comparison between equations (2.81) and (2.79) yields a function 
g{4>k, 9k) that makes both equations equal: 

g{(f>k,0k) = - — -£^*''(/sin6'fc,0fc)f(0fc,6'fc)e*^'^'' {l- H{9k- 9^)) J — {cos9k)^^^ 

(2.82) 

where H{9k — 9^) is the Heaviside step function. With g{(j)k, 9k) given by equation 
(2.82), the held at the focus of a lens can be given by: 

niv /*27V 

= / sm{9k)d9k d(j)k gi(j)k,9k)'R.i(j)k,0k)dpexp (ikz) (2.83) 

Jo Jo 

Now, I will show that has a well-dehned if the incoming held does, too. That 
is, I will show that the aplanatic model of a lens preserves by verifying the following 
condition: 

If J^[E“]=mE“ ^ J,,[E^°"] =mE^°" (2-84) 

A paraxial beam with a well-dehned helicity p and an value of Jz = is given by 
any expression of the kind: 

gpar _ E^’^^{p)e^^‘f^exp{ikz)dp (2.85) 

^^This is an abuse of language that will be used throughout the whole thesis. It means that the 
beam is an eigenstate of Jz with eigenvalue m. 
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If this expression is substituted in equation (2.82), the following expression for g{(j)k, Ok) 
follows: 

giOkAk) = sin ek)e^^^^^'^^H{ct>k, Ok) (l - - 0^)) 

( 2 . 86 ) 

Now, in order to prove that the aplanatic lens model does not change the content 
of EP^'', I apply Jz to E^°'^ when g{(j)k,0k) is given by equation (2.86). The fact that 
E^°^ can be expressed as a rotation helps to compute the result. As seen by equation 
(2.36), Jz can be applied as a partial derivative with respect to (j)k under the integral 
operation. Then, it can be seen that: 

sm{9k)d9k [ d0fc R(0fc, 9k)id^ exp {ikz) (2.87) 

do do 

with g{(l)k,9k) only depending on via two terms, and t{(j)k,9k)- As one 

could expect, it is easy to see that the aplanatic model preserves of E™ if and only 
if t{(j)k,9k) = t{9k)- Throughout this thesis, I will consider that all lenses (as long as 
they are centred with the beam axis) preserve and therefore t{9k) = t^{9k) = t^{9k)- 

2.7 Overview 

In this chapter, different theoretical methods to describe light-matter interactions have 
been explained. They will be used in the remaining of this thesis, both to understand 
the experimental results, as well as to carry out different theoretical developments. 
Next, I present a summary of the methods that have been described above as well as 
their applications in the remaining of the thesis. 

In section 2.1, Maxwell Equations for a linear, non-dispersive, homogeneous, isotropic 
and source-free medium have been introduced. These equations will be used in chap¬ 
ters 3, 4 and 5 to deepen in understanding of the interaction of vortex beams with 
spherical particles. 


In section 2.2, the paraxial approximation, which is a very relevant limit of Maxwell 
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equations, has been explained. This approximation is going to be used in all the ex¬ 
perimental chapters (6, 7, and 8). Then, the concepts of EM held, operator, and 
symmetry have been given. Of particular interest is the concept of symmetry, which is 
used throughout this thesis to predict theoretical results or interpret experimental re¬ 
sults. A complete list of symmetries with their respective generators can also be found 
in section 2.2. Many of these generators will be used to characterize the light beams 
as well as the samples that will be used in this thesis. In particular, chapters 4, 5, 
7, 8 describe different phenomena that arise when beams with duality and cylindrical 
symmetry interact with samples with mirror and cylindrical symmetry. 

In section 2.3, a general method to construct solutions of the vectorial Helmholtz 
equation is sketched. Then, plane waves are constructed using this method. Plane 
waves are characterised in terms of their symmetries, and they are used to construct 
Bessel beams, which are symmetric under some different symmetries. Both plane waves 
and Bessel beams are later used in section 6.2. 

Section 2.4 introduces two different basis of multipolar helds with two different radial 
functions each. The understanding of the symmetries of these beams will be crucial 
in chapters 3, 4 and 5, where they are used to solve Maxwell equations in a spherical 
domain. 

Section 2.5 presents in a formal way the key idea behind the work done in this thesis, 
which was introduced in 1. That is, controlling and/or characterising the scattering 
of nano-structures using symmetric light. This method is exploited both theoretically 
and experimentally in chapters 4, 5, 7 and 8. It is seen that symmetric beams can be 
used to effectively turn a non-dual particle into dual, or inducing circular dichroism in 
a non-chiral smaple. 

Finally, in section 2.6, the description of the aplanatic model of a lens is given. The 
aplanatic model is used in chapter 4 to theoretically control the multipolar content 
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of a beam. Furthermore, the symmetry conditions proven in 2.6.1 and 2.6.2 are then 
assumed in chapters 7 and 8 to understand the experimental results. 
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"/ learned very early the difference between knowing the name of 
something and knowing something." 


Richard P. Feynman 


3 

Generalized Lorenz-Mie Theory 


3.1 Introduction 

In this chapter, I will use the tools described in the preceding chapter to solve the 
electromagnetic problem described in Figure 3.1. An arbitrary electromagnetic field 
excites a single homogeneous, isotropic sphere embedded in a homogeneous, isotropic 
and lossless medium^. As mentioned in the introduction, the scattering of a sphere 
illuminated by a plane wave was independently solved by Ludvig Lorenz in 1890 and 
Gustav Mie in 1908 [97, 98]. Unfortunately for Lorenz, Mie’s name has stuck in the 
literature and nowadays the problem is mostly known as the Mie scattering problem 
[99]. A very complete review about the history of this problem can be found in [48]. 

^Note that denoting the incident field as E‘ is an abuse of language. The incident electromagnetic 
field has also a magnetic component H', and it is indeed used to solve the Maxwell’s equations. 
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Figure 3.1: Schematics of the Generalized Lorenz-Mie problem. An arbitrary EM field 
E‘ (red) propagates through a lossless, isotropic and homogeneous medium with optical 
constants The held interacts with an isotropic, homogeneous sphere with optical 

constants (e*, /Us). As a result, the incident held is scattered in all directions giving rise to E® 
(blue). 


Before the invention of lasers in 1958 [11, 12], Mie theory had already been successfully 
applied to a very wide range of helds, e.g. chemistry, material science or atmospheric 
physics [43, 84]. However, once Mie theory started being tested with lasers, it became 
clear that a more general formulation of the problem had to be established. Although 
in many cases the light beam produced by a laser could be modelled as a plane wave 
with a good approximation, a more general theory that took into account the spatial 
non-uniformity of light beams was needed. In this context, the Generalized Lorenz-Mie 
Theory (GLMT) was created in the 1980’s [100]. The GLMT development was espe¬ 
cially led by Gerard Gouesbet and Gerard Grehan, both Professors at the University 
of Rouen, France. In its hrst ten years, the progress of GLMT was mainly focused on 
the theoretical part. In 1985, Gouesbet et al. managed to solve the scattering problem 
of a free propagating Gaussian beam interacting with an on-axis arbitrary sphere using 
Bromwich potentials [101]. And then, in 1988, a seminal paper was also published 
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by Gouesbet and co-workers where the scattering of a sphere located in an arbitrary 
position with respect to the axis of propagation of a Ganssian beam was solved [102], 
This paper laid the fonndations of the GLMT, as it solved the scattering problem of a 
single sphere excited with an arbitrary beam of light [103, 104]. Since then, as the the¬ 
oretical framework had already been laid, the development of the theory happened in 
the applications and compntation side. Dne to compntation limitations at that time, a 
lot of effort was pnt into hnding ways of efficiently compnting the so-called beam shape 
coefficients of an arbitrary beam [105]. In addition to that, lots of applications were 
found, including phase-Doppler instruments, imaging, optical characterization and op¬ 
tical manipulation, among others [100]. In this chapter, I will summarize the main 
theoretical results of GLMT and interpret them from a symmetry-based point of view. 

3.2 Symmetry considerations 

The wide variety of applications of GLMT is a proof of the success of the theory. How¬ 
ever, due to the typical method of solving the problem via vector potentials (Bromwich, 
Debye or Hertz) [43, 48], the symmetries of the system are usually blurred. In order to 
gain a bit of insight about the symmetries in the problem, let’s take a look at Figure 
3.2, where the different symmetries of the Mie problem are presented. I have ruled 
out the incident beam of any consideration, as I consider it as an excitation to the 
system, but not part of the system itself. The system is formed by the sphere and 
the surrounding medium. The origin of the frame of reference in which the problem 
is described is placed at the centre of the sphere. This is necessary for the three last 
symmetry considerations to hold: 

• Temporal translations, Because both the sphere and the medium surround¬ 
ing it are assumed to be linear media, the excitation of the sphere is independent 
of the instant when the sphere is excited. 

• Rotations around the z axis, R^. A rotation around the 2 ; axis leaves the system 
invariant. This will also be referred to as cylindrical symmetry. 
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b) 





Figure 3.2: Symmetries in Mie Theory, a) Temporal translational symmetry. Two exact 
plane waves interact with the sphere at different instants {t for the blue wave, t* for the red 
one). The interaction is the same one for both, b) Rotations around the z axis. The system 
remains invariant when it is rotated around the z axis, c) Independence of orientation of 
reference frame. The problem is cylindrically symmetric around the z axis independently of 
the orientation of the reference frame, d) Parity. A parity r — ?• —r change of coordinates 
leaves the problem invariant. 


• Independence of orientation of the reference frame. As long as the origin is the 
center of the sphere, the orientation of the axis will be irrelevant. That is, any 
axis is an axis of cylindrical symmetry. 


• Parity, If. The system is invariant under parity transformations, i.e. changes of 

r —y —r. 


These four transformations commute with each other. Thus, the natural modes of the 
system will be eigenmodes of these four transformations. Now, due to the fact that 
continuous symmetries can be expressed as a function of their generators, I will use 
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the generators of the three hrst symmetries to classify the natural modes of the system 
(see section 2.2): 

• Hamiltonian, H. Temporal translations are generated by the Hamiltonian Ta* = 
exp{-iHAt). 

• Projection of the AM in the z axis, J^- The z component of the AM generates 
R^{d) = exp(-zJ^6'). 

• The AM squared, J^. The transformation Tj{6) = exp{—iJ‘^6) is meaningless. 
Nevertheless, the independence of orientation of axis implies that the eigenmodes 
of the system need to be eigenstates of J^. 

Parity (H) is not a continuous transformation, therefore I will use the transformation 
itself to classify the eigenmodes of the system. Now, in 2.4 we have seen that one of the 
complete basis of solutions of source-free Maxwell equations are the multipolar modes 
As shown in there, these modes are eigenmodes of H, J^, J^, and H 
respectively. Thus, the natural modes of the system depicted in Figures 3.1, 3.2 will be 
the multipolar modes, as they fulhl all the symmetry requirements. A harmonic time 
dependence will be assumed throughout the rest of the chapter, as non-linear 
phenomena are not going to be studied. 

3.3 The scattering problem 

In this section I will solve the problem sketched in Figure 3.1. That is, I will hnd 
the EM held (E, H) in all space for a given arbitrary incident held E5 In order to 
achieve that, I will describe all the EM helds as a superposition of multipolar modes. 
This will simplify all the calculations, as they are the natural modes of the system. 
The procedure will be the following one. First, I will break down the domain into two 
sub-domains which will be separated by the surface of the sphere. The EM held inside 
the sphere will be given by E‘“* and the interior electric and magnetic helds. The 
expression of the EM held outside the sphere will be given by the total held, which 
can be computed as the addition of the scattered held (E®, H®) and the incident held 
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(E', H'), i.e. = (E® + E‘,H® + ff). Then, I will impose that the tangential 

components of the helds inside and ontside the sphere are continnons [46]. To make 
the problem as general as possible, I will snppose that the incident electric held is the 
most general snperposition of mnltipolar modes: 


E‘ 


E E 

j=i 


(m) A (m) 
^ jrriz'^jmz 


yjmz 


+ A 


(e) 

jm^ 


mz=-j 


(3.1) 


where and g^^^^ are some general coefficients modnlating the mnltipolar snperpo¬ 
sition. Then, the magnetic held can be obtained nsing the Maxwell equations (2.7) 
and the relations between multipolar modes given in [51]: 




V X E‘ _ -1 

OJfl 


iojgik 


E E 

,j=l m^=-j 


M A h) 
yjmz jrriz 


(e) aM 

^jmz jrriz 


(3.2) 


The scattered and interior helds are also supposed to be a general superposition of 
multipolar helds, whose expression are given by: 


E" 

OO 

= E 

3 

E 


(3.3) 


j=i 

rtz=-j 




-1 

UJH 

OO 3 

h- -a- A^”"^ 

j=l m^=-j 

(3.4) 

j^int 

OO 

= E 

3 

E 

r- 4-d- 

(3.5) 


j=i 

nz=-j 



Hint 

-1 

cn/ii 

OO 

E 

-i=i 

c- -d- A^”"^ 

mz=-j 

(3.6) 


where jj^i and /i are the magnetic permeabilities of the sphere and the surrounding 
medium respectively. However, as mentioned in 2.4, even though we use the same 
notation for the multipolar modes in scattering and free propagation, their radial func¬ 
tions are diherent. For scattering, the radial function is a Hankel function of hrst kind. 
Whereas for the interior and incident helds, the radial function is a Bessel func¬ 
tion, jj [48, 68]. Once the general expression of the helds has been found, we apply the 
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tangential Maxwell bonndary conditions, which yield the following fonr eqnalities: 


rpint 


H't 


— El + E\ 


= m + Hi 


ET = 


El + E\ 


(3.7) 


rrint 


= m + H\ 


Eqnations (3.7) can be solved with Cramer’s rnle, yielding the following results: 






(e) 


(m) 

Cjdjrn. 






I, (™) 

^jSjniz 


(3.8) 


where {aj, bj,Cj, dj} are the so-called Mie coefficients, whose expressions are given by: 


ttj — 


b,= 


^,nljj{nrx)[xjj{x)\l - ^,ijj{x)[nrxjj{nrx)\ 
^nlj j{nrx)[xh^P {x)]' — {x)[nrxjj{nrx)]' 

^iijj{nrx)[xjj{x)'\l - ^lij{x)[nrXjj{nrX)] 

jjLii j{nrX)[xh^P {x)]' — jJih^P {x)[nrXj j{nrX)]' 


(3.9) 


liijj{x)[xhf\x)]' - ^iih^^\x)[xjj{x)]' 
^ij j{nrx)[xh^P {x)]' — {x)[nrxjj{nrxy\' 


(1)/ 


C = 


^iinrjj{x)[xhP{x)]' - ^iinrhp{x)[xjj{x)] 
jjLnli j{nrx)[xhP [x)]' — jjLihP {x)[nrxjj{nrx)]' 


(b/ 


dj 


where x = 2 tiR/\ is the so-called size parameter, with R being the radius of the 
sphere in consideration and A the excitation wavelength, inversely proportional to 
the frequency ca; Ur is the relative refractive index of the sphere with respect to the 
surrounding medium, It can be observed from equations (3.8) that once 

the coefficients have been found for the incident field, the solution of the 

scattering problem is straightforward: the expression for the scattered and interior 
helds are a mere formal copy of the incident held, where each multipole is modulated 
by a Mie coefficient. Hence, the key point to solve the Generalized Lorenz-Mie problem 
is hnding the coefficients Both and are usually known as beam 

shape coefficients, except for a subtlety. In GLMT, the beam shape coefficients are 
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defined as and and they are equal to and divided by a factor [48]: 


(m) 




(e) 
3jm^ 




(3.10) 


where are the plane-wave coefficients: 


vw ^ ^ / -v? ‘^3 + ^ 


(3.11) 


Now, given Eh we can hnd them applying the techniques shown in (2.3). That is, due 
to the orthonormality relations between and we can hnd g^^^ and g^^^^ 

projecting E' onto the multipolar modes, yielding a unique result. Thus, from equation 
(3.1), we obtain: 


(m) 

9jm^ 


9jm^ 


>1 r^TT 


-IJO 
1 /■2tt 


(A;”'J-.E' d(cos9)d.#. 
(A'l)--E' d(cos9)d.#, 


(3.12) 

(3.13) 


J -1 Jo 

Equations (3.12, 3.13) give us the formal solution of the problem. Hence, solving the 
GLMT problem basically consists in computing two sets of two-dimensional angular 
integrals. With our current computational power, and g^^^ can be computed using 
standard numerical integration techniques, without too much of a struggle. Neverthe¬ 
less, as mentioned in 3.1, that was not the case when GLMT was hrstly formulated, and 
so a lot of effort was put into computing the beam shape coefficients. A very complete 
review of all the different available techniques to compute the beam shape coefficients 
can be found in [104]. I will not explicitly comment further on those techniques, as in 
the proceeding sections I will simplify (3.12, 3.13) thanks to some symmetry consider¬ 
ations. 


3.4 Cross sections 

In the previous section I have given the expressions of the EM field in all points of 
space. Nonetheless, it is difficult to measure the EM field in experiments at optical 
frequencies. Therefore, we need to derive some quantities that we can experimentally 
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determine, even though they might provide a less complete description of the scattering 
process. One of these quantities is the EM power carried by the radiation, since many 
light detectors are sensitive to it. Thus, we need to mathematically characterize how 
much power is carried by the three fields in the problem. In order to do that, we first 
define the time-averaged Poynting vector^ [78]: 

(S*) = ^Re [E* X (H')*] (3.14) 

Now, the total fields are obtained adding the scattered and incident fields, he. (E*, H*) = 
(E® -h E', -|- H'). Hence, we can split the time-averaged Poynting vector in three 

parts: 


(S*) = 

(S') + (S^) + (S""*) 

(3.15) 

(S^) = 

^Re [E^ X (HT] 

(3.16) 

(S^) = 

^Re [E'' X (ff )*] 

(3.17) 

^gext) _ 

^Re [E' X (ff)* + E" X (H')*] 

(3.18) 


The flux of the Poynting vector across a surface gives us the energy flux per unit of 
area and unit of time, which accounts for the loss of energy inside the surface. We will 
denote this flux as minus the absorption cross-section: 

-Cabs= [ (S*)-ndH (3.19) 

Js 

where dH is a differential surface element, with n being the normal vector to this 
element at every point. Due to the distributive property of integrals, this can be 
also split in three terms, to which we will refer as incident, scattering and extinction^ 

^Remember that a harmonic exp(—iwt) dependence is supposed 

^The extinction cross section is defined as energy loss as well, that is why define it with a minus 


sign. 
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cross-sections: 


Cabs 

— Ci -\- Cg -\- Cext 

(3.20) 

a 

= [ (S') ■ h dn 

Js 

(3.21) 

C*. 

= [ (S") ■ h dfi 

Js 

(3.22) 

Cgxt 

= [ (S""*) ■ n dn 

Js 

(3.23) 


Clearly, if we consider the flnx of the Poynting vector across a spherical snrface, Ci = 
0, as the snrronnding medinm is non-absorbing. The choice of a spherical snrface 
simplihes the calcnlations withont any loss of generality, therefore I will nse it for 
the rest of the calcnlations. This will directly imply that the snrface element dfl = 
d (cos 6) d0, and the normal vector n will point in the radial direction f. I will start by 
compnting Cg- In order to do this calcnlation, the following orthonormality relations 
of the mnltipolar helds will be considered [50]: 


/ 



j (a“. X (a;1)-) . f dSJ = 0 


(3.24) 

(3.25) 


where F{r) is 


F(r) 


jlO+iP + (j + i)IO-iP 

2j + l 


for y = m 
for y = e 


(3.26) 


and is a Bessel fnnction for the incident and interior held and a Hankel fnnction of 
the hrst kind for the scattered held. Then, it can be proven that 

oo j 

= (l“>41t+l''i9Sl^) (3-27) 

j=l m^=-j 

where the two following relations have been nsed [42, 50]: 
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Note that equation (3.27) differs from the expression given in [48] in some factors due 
to the fact that the normalisation of the multipolar modes given by equation 2.52 is 
not considered in [48] using the Bromwich potentials. Very similarly, the following 
expression can be obtained for Cext [48]: 


Cp.T.t — 2 Re 


+ hj\g 

J=l m^=-j 


(m) 12 

jmz I 


(3.30) 


Consequently, Cabs can be obtained subtracting expressions (3.30) and (3.27): Cabs = 
Cext ~ Cs- 


3.5 Efficiency factors 

As I have shown above, the expressions for the interior and scattered field depend on 
the incident field. In particular, their expressions depend on the intensity carried by 
Eb This fact naturally translates into the cross sections, i.e. the cross sections depend 
on the intensity of the incident beam. In order to remove the intensity dependence 
from the cross-section expression, we normalize the cross sections over the intensity of 
the beam across a transverse section [48, 106]: 

where v can take three values v = I will not give general expressions for the 

efficiencies. Instead, I will either explicitly compute them (see chapters 4 and 5), or 
measure them (see chapter 8). A subtlety about the measurements in chapter 8 is in 
order at this point. Despite normalising the scattering intensity over the power of the 
incident beam (Jii°™)^ the result is modified by the limited NA of the detector. That 
is, the definitions of cross section and efficiency factors given by equations (3.19-3.23) 
and (3.31) suppose that the scattering in all angles (0,0) is collected. However, this is 
not the case in the experiments done in chapter 8. There, only scattered light in the 
backward semi-space is collected. Now, even though backward cross sections do not 
generally have the same A dependence as proper cross sections [46] , the features that I 
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will be especially interested in will still remain. Namely, as it will be shown in section 
4.4, the excitation of spheres with plane waves hides some resonant behavionr that can 
be unveiled using vortex beams. 



“I think it only makes sense to seek out and identify structures of 
authority, hierarchy, and domination in every aspect of life, and to 
challenge them; unless a justification for them can be given, they 
are illegitimate, and should be dismantled, to increase the scope of 
human freedom." 

Noam Chomsky 


4 

Excitation of single multipolar resonances 


4.1 GLMT with cylindrically symmetric beams 

In the previous chapter, I have solved the Generalized Lorenz-Mie problem. That is, I 
have found the EM helds in the whole space when a homogeneous and isotropic sphere 
is excited by an arbitrary monochromatic light beam Eb In the present chapter, I will 
particularize the solution given by equations (3.3-3.13) for an incident monochromatic 
cylindrically symmetric beam. The dehnition of a general cylindrically symmetric beam 
stems from the dehnition of cylindrical symmetry. As it has been described in chapters 
2-3, a light beam is symmetric under a transformation T when its EM held is an 
eigenstate of the transformation generated by T: Tu = c„u. I dehne a cylindrically 
symmetric beam as an eigenstate of Rz{0), the operator of rotations around the ^ axis. 
Now, because Rz{9) can be expressed as Rz{9) = exp{—iJz9), a cylindrically symmetric 
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beam can be alternatively defined as an eigenstate of Jz- As we have seen in chapter 2, 
there are different families of monochromatic modes that have a well-defined Jz- Both 
Bessel beams and mnltipolar fields are basis of modes that are eigenstates of Jz with 
eigenvalnes m and respectively (see eqnations (2.37, 2.49)): 


Jz^mkp 

( 4 . 1 ) 

Jz-^jmz '^z-^jiUz 

( 4 . 2 ) 


An arbitrary snperposition of them (keeping m or ruz constant) will also be an eigen¬ 
state of Jz- Now, both basis are valid to describe the problem, bnt I will nse the 
mnltipolar fields as they are the normal modes of a sphere (see section 3.2). In this 
basis, the expression for the electric field of a general cylindrically symmetric beam is: 
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where and are proportional to the beam shape coefficients, (see section 3.3). 
Note that a general snperposition of mnltipoles is given by eqnation (3.1). That de¬ 
composition contains two snmmations, on j = l,..,cx) and rriz = respectively. 

In contrast, expression (4.3) only involves a snmmation in j = m*,..,cxD. First, rriz 
needs to be fixed to m* so that is still an eigenstate of Jz- And secondly, dne to the 
fact that J^ > Jz, rn* bonnds j to valnes j > m*. 

Last bnt not least, I will also impose that the cylindrically symmetric beams are eigen¬ 
states of the helicity operator, A. In this way, they will be symmetric nnder dnality 
transformations and it will be easier to stndy the dnality properties of the system (see 
chapter 5). Imposing that the incident field E^ has a well-defined helicity has farther 
implications on expression (4.3). As it has been demonstrated in chapter 2, a mode 
of light with a well-defined helicity can be expressed as a snperposition of electric and 
magnetic modes. Conseqnently, the two families of beam shape coefficients can be 
rednced to one: 


r - - 

yjrriz 


(e) 


(4.4) 


where Cjm^p is the only family of beam shape coefficients remaining and p = ±1 are 
the helicity valnes for the beam. Thns, the mnltipolar decomposition of a cylindrically 
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symmetric beam with a well-defined helicity is: 

OO 

E' = E (a;™ 1 + ipA^.) (4.5) 

j=\m*\ 

Equation (4.5) will be the starting point for this chapter. I will use it to solve the 
GLMT problem with paraxial cylindrically symmetric beams in 4.2 and with general 
non-paraxial cylindrically symmetric beams in 4.3. Then, in 4.4, I will show how the 
use of these beams of light along with the proper choice of ratio between the particle 
size and the excitation wavelength is enough to excite WGMs [107]. 

4.2 Paraxial excitation 

In this section I solve the GLMT when the incident held E' is a Laguerre-Gaussian 
beam, LG;^g^. The LG modes are the normal modes of a cylindrical cavity and can 
also be obtained as superpositions of Hermite-Gaussian modes, which are the normal 
modes of spherical-mirror resonator^ [45] . They are a basis of solutions of the paraxial 
equation in cylindrical coordinates [12, 108], i.e. any paraxial beam can be described 
as a superposition of LG modes [109]. The LG modes had been known by the optics 
community for a long time, but their use highly escalated after a seminal paper by 
Allen and co-workers in 1992 [110]. In that paper, it was established that there was 
a relation between the phase and the AM content of these modes. This laid out the 
grounds of the held of the angular momentum of light. Since then, the AM of light has 
been used in many diverse helds such as quantum optics [31, 111], optical manipulation 
[112-115], optical communications [116, 117] or astrophysics [118, 119]. The general 

^LG beams are also named as vortex or doughnut beams. However, there are some subtleties. 
The ‘vortex beam’ term usually refers to beams of light whose phase dependence is of the kind 
exp [ilcj)]. This phase dependence generally gives rise to an intensity profile that resembles that of a 
doughnut. Nevertheless, a beam can have a doughnut-like intensity profile without having a exp [ilfj)] 
phase dependence. The LG modes are a well-defined set of modes (see equation (4.6)). Their phase 
dependence is that of a vortex beam, but their intensity profile is not always doughnut-like. In fact. 


their intensity profile is only doughnut-like when the radial parameter <7 = 0. 
^In fact, LG modes are also normal modes of this kind of resonator. 
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expression for a LG modes is given by: 

2p2 


= Ni^q exp 


-P 


w{zY^ 


Jt’- 

P 


_ 

w‘^{z) ^ I ta(2:) j 


/+! 


2 (z 2 !j./ 2 ) +(2g+^+I) tan ^{z/zq)+1(I> 


(4.6) 


where g > 0 is the radial index, and I the azimuthal one. The azimuthal number I 
describes how the phase winds around the center of the beam, whereas q has to do 
with the number of intensity nodes in the radial direction, (p, 0, z) are the cylindrical 
coordinates; are the generalized-Laguerre functions; w{z) is the waist of the beam 


at a plane of z constant; Ni^g = 


1/2 


is a normalization constant; and k 


[(7r(/ + g)!_ 

is the wave-vector. Note that due to the fact that a LG mode is a paraxial beam, 
the polarization is added to the beam independently of its spatial shape (see section 
2.2) [53]. Next, I will decompose a circularly polarized LG beam into multipoles. In 
fact, a circularly polarized paraxial beam has a well-dehned helicity within the paraxial 
approximation [54, 120]. That is, if we apply the helicity operator A to a circularly 
polarized LG mode, and we only consider paraxial terms, we obtain the same mode 
times the helicity value of the beam: 


A [LGi^gdp] ^ p [LG;,g(Tp] (4.7) 

with dp being a general circularly polarized vector defined as dp = {k + ipy) /\/2, and x 
and y the horizontal and vertical polarization vectors. The differential form of A is given 
by the expression A = (Vx)/fc [40, 54]. Then, the multipolar decomposition of LG^^gdp 
is given by equation (4.5) and only one set of beam shape coefficients is needed, Cjm^p. 
Before computing them, one comment needs to be made. A beam such as LG^^^dp does 
not fulhl Maxwell equations. However, its decomposition into multipoles using equation 
(4.5) does. This is a known paradox in optics that was solved by Lax in 1975 [48, 53]. 
The idea is that the Maxwell held given by the multipolar decomposition is equal to 
the paraxial beam only in the paraxial limit. Furthermore, Lax also showed that a 
paraxial held gives rise to many diherent helds fulhlling the Maxwell equations. The 
general formalism to hnd one of the possible multipolar decompositions of a circularly 
polarized LG mode was given by my supervisor in [120]. Then, in [121], we applied 
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that formalism to give an analytical formnla for the beam shaping coefficients Cjm^p- 
Onr hndings are shown next. If = LG/^gdp, then we can re-express it as: 


OO 



j=\m*\ 


( 4 . 8 ) 



where m* = I + p stands for value of The reduced rotation matrix d?rn%pi.^k) can 
be found in [51], and the function fi^q{kr) is related to the Fourier transform of the 
incident held, with kp being the transverse momentum, i.e. kp = ^k^ + ky = ksinOk- 
Note that the function Cjm^p is still dehned as an integral, although it is no longer 
a double one, unlike (3.12, 3.13). This is a consequence of the fact that E' and the 
multipolar helds are eigenstates of Jz with value m*, therefore projections are only 
needed on values. Also, note that the dehnition of Cjm^p in equation (4.8) is slightly 
different from its dehnition in (4.5): a factor P(2j -|- 1 )^A been made explicit. The 
dehnition of Cjm*p given by equation (4.8) will be the one used for the rest of the thesis. 

Actually, equation (4.8) is valid for any beam whose Jz is well-dehned. For the specihc 
case of a beam, the integral in equation (4.8) can be solved analytically thanks to 
the expression of fi^q and provided the integral 4.8 is extended to inhnity. The following 
expression is obtained: 


c,m:r = (-i)‘ l(i + P)'U - P)'U + * + P)'U - i - pY]'''^ 



( 4 . 9 ) 


s 


^Remember that as described in chapter 2, the AM content of a paraxial beam can be obtained by 
adding the topological charge of the beam {I for a beam) and the value of the helicity p. This is 

what is usually described as AM = 0AM + SAM, i.e. the total angular momentum can be obtained 
adding its orbital and spin parts 
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where the number Ma{fi,q) = kpdkpkpfi^g{kp) is the momentum of order a of the 
function fi^g. The expression for fi,q{kp) when the incident beam is a l^Gi^q is: 

( 2 I \ 1^1 

27r(|/°+g)! j (woK/\/2j L\^^{wlkl/2)ex^{-wlkl/4) (4.10) 

with wq = = 0) being the beam width of the LG mode in real space. Finally, I 

give the expression of Ma{fi^q) when the function /; g is given by the expression (4.10). 
As it is shown, in this particular case, the integral is analytical. Also, I particularize 
the result for g = 0: 




X 


(4.11) 


M,+2,(/,,„) = (s + (4,12) 

where r(z) is the Gamma function and 2 ^ 1 ( 0 , 6, c : z) is the hypergeometric function. 
The problem, then, is solved. Using equations (4.8, 4.9) and (4.11) the decomposition 
of the incident beam is found. Then, using the GLMT formulation introduced in the 
preceding chapter, the scattered and the interior fields are obtained (see section 3.3): 


E" = 






(4.13) 


j^int _ 


■!i(2j + l)‘-'=C,. 


jmlp 








As already mentioned in chapter 3, it is clear that the solution of the scattering prob¬ 
lem only depends on the beam shape coefficients Cjm*p- The beam shape coefficients 
depend on the parameters of the incident LG modes, i.e. the azimuthal and radial 
numbers {l,q), the helicity p, and the waist of the beam wq. That means that once 
the LG mode and its polarization have been chosen, the only way of tailoring the 
multipolar content of the beam is with the width of the beam, wq. However, wq can¬ 
not get arbitrary values, as the beam still needs to fulfil the paraxial equation [53] . As 
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a consequence, the multipolar content of the beam is almost hxed for every single mode. 

In Figure 4.1, the multipolar content for four different LG modes is depicted. It can 
be observed that over 40 multipolar modes of each parity are needed to describe a LG 
mode of an arbitrary order^. This fact will be further explored in the next section. I 
advance that using high NA lenses will enable us to describe a cylindrically symmet¬ 
ric beam with a handful of multipoles. This fact will be a key point to excite single 
multipolar resonances in section 4.4. 

4.3 Non-paraxial excitation 

In this section, I will extend the results of the previous section to the non-paraxial case. 
That is, in this case the incident held E' will fulhl Maxwell equations. Mathematically, 
the only difference is the way the beam shape coefficients Cjm^p are computed. How¬ 
ever, physically there is a very big difference: the use of an aplanatic lens enables us to 
control the multipolar decomposition at the focal plane with a much higher precision. 
This fact is of great interest to control multipolar resonances both in nanostructures 
and atoms [73, 107, 122-124]. I will further develop this idea on section 4.4. In addi¬ 
tion, I will show that by adjusting the focal distance of a lens, the magnitude of the 
scattering efficiency can be increased by orders of magnitude. 

To start, I will show how the aplanatic model of a lens tailors the multipolar con¬ 
tent of a beam. In section 2.6, I have demonstrated the conditions under which the 
aplanatic lens model does not change the helicity and the angular momentum con¬ 
tent {Jz) of the incident helds. In order to do that, I have expressed the aplanatic lens 
transformation [49, 95] as a rotation in the momentum space. Next, I will use the same 
mathematical technique to retrieve the beam shape coefficients Cjm^p as a function of 
the focal distance of the lens, the NA, and the media surrounding the lens. On one 

^Multipolar modes of both parities are excited as Cjm*p gives rise to the two families of beam 
shape coefficients (see equation (4.4)). 
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Figure 4.1: \Cjm*p\‘^ for different beams, a) LGo,o b) LG 2 ,o c) LG 4 ^o d) LG6,o- 

The relation between the width (wq) of the beam and the wavelength is w/X ~ 4. The 
helicity is p = 1 for the four different cases. The insets are intensity plots of the multipolar 
decomposition for each of the four cases. The intensity of each of the LG modes is then 
retrieved. 


hand, I will consider a general snperposition of circnlarly polarized propagating plane 
waves with the same helicity: 


El = 


2n 


sinOhdOk / d(j)k g{0k,(j)k)'R.{0kAk)d-pe 


ikz 
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0 


(4.14) 
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On the other hand, the aplanatic model can be modelled as the integral (2.70) of a 
collimated (or paraxial) beam Einc over the surface of a sphere [49]: 


Eo — 


ikf ex.p{—ikf) 


r-27r 


27r 


sinOk d6kd(j)k'E^{0k,(pk)e 


{ikz COS 9) ^{ikp sin 9 cos((/)— 


ikf exp{—ikf) . 


c2tt 


27r 


sinOk d6kd(j)k'E,^{6k,4>k)exp{ik-r) 


(4.15) 


In chapter 2, I proved that the polarization transformation produced by equations 

(4.15) and (4.14) is the same one provided the two transmissivities of the lens are 
equal, k^{0k) = fP^Ok) (see section 2.6). I also gave the value of the function g{(j)k,dk) 
that makes Ei = E 2 , which is (see equation 2.82): 

9{0k,4>k) = --E-(/sin0fc,0fe)t(4)e'^^'“ (l - 77(0^ - C)) 

(4.16) 

where / is the focal of the lens; 9^ = arcsin(NA), with NA the numerical aperture 
of the lens; rii and n are the index of refraction of the lens and the medium sur¬ 
rounding it respectively; H{6k — 9^) is the Heaviside step function; and 77*” is the 
value of the paraxial incident field at the back-aperture of the lens. Note 77*”(p, 0) = 
77*”(/sin6'fc, 0fc). That is, there is a match between the cylindrical coordinates in the 
real space of the paraxial incident beam, and the spherical coordinates of its angular 
spectrum in the momentum space. As mentioned in section 2.6, this fact is a key 
factor to prove that the aplanatic model does not change the helicity nor the AM of 
the incident hend. Now, using the properties of the multipolar fields under rotations 
[51], the expression (4.14) can be re-written as: 


E = / sin 


(*27T 

in^fcd^fc / d0fc g{9k, (pk) ^ *^(2j + 9k) + ipAj' 

7o , _ L 






- 27 r 


= / sm 


'0 


in9kd9k I dPk g{9kAk)J2^^{2j + 

j,mz 

+ ip A) 


aJ-^I + ipA'). 


(e) 

jmz 




]m^p 






Orriz 


J,mz 


( 4 . 17 ) 
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with 


pTi ^27r 

Cjm.p= / sin^fcd^fc / d(j)k g{0k,(pk)di,,p{0k)e~^'^'^'‘ (4.18) 

Jo Jo 

If the expression for g{9k, (j)k) found in (4.16) is used in the previous expression (4.18), 
the value for the Cjm^p is the following one: 

r^k fru , 

C,m.p = / sm9kd0j-—. -cos^/^9kdl^p{0k) / d0fce-*™'^'=e*^’'^'=E-(/sin0fc, 0^) 

Jo V ^^2 Jo 

(4.19) 

As I mentioned in the beginning of the chapter, I will be interested in cylindrically 
symmetric beams. The helicity has been set to p already, so the only other requirement 
on the incident paraxial field EP°'^{p,(j)) = sm9k,4>k) is having an azimuthal 

dependence of the kind In general, this can be achieved by a superposition of LG 
beams with a fixed azimuthal number /: 


Er{fsm9k,cPk) = 5^c,LQ,,(/sin0fc)e*''^'= (4.20) 

9=1 


However, for simplicity, the calculations will be done for a single LG/^g mode. Now, 
using the following orthogonality property of the exponentials. 


27 r 

dx = 6mm' (4.21) 

it can be seen that the second integral in (4.19) can be analytically solved, setting 
rriz = m* = I + p- Then, Cjm^p can be expressed as a 1-dimensional integral. Its 
expression, as well as the one for the incident field Eb is given by the following equation: 
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j=\m*J 




jmlp 
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jmlp 
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(4.22) 


exp 


-(/Bin4)n ( V2 


w{zy 


^+1 


w‘^{z) \w{z) 


+ (2g + / + l)tan ^{z/zQ)j 

Next, I show the influence of the NA of a lens in the beam shape coefficients. In 
Figure 4.2, I have plotted the multipolar decomposition of the same four LG beams 
used in Figure 4.1. The four LG modes are plotted for two different focusing scenarios, 
one of them being focused with NA= 0.25 (yellow), and the other one with NA= 0.9 
(red). In all cases I have used LG modes with a well-defined helicity p = 1. The width 
of all beams in the paraxial approximation is chosen so that the entrance pupil of the 
lens is filled. The results presented in Figure 4.2 show that the higher the NA, the 
narrower the distribution of multipolar modes is. For some cases such as Figure 4.2(a), 
the beam shape coefficients Cjm*p can be described with only a handful of multipoles. 
In those cases, the value of \Cjmip\‘^ increases considerably, as the incident field needs 
to be normalised on a transverse plane [120, 121]. As a consequence, the following 
condition needs to be fulfilled: 


exp 


-k{f sin 9kf‘z 
2(Z + 


^(2j + 1)|G,>„*p|2 = 2 (4.23) 

j 

The other parameter to control the distribution of multipoles is the Z of the incident 
beam. As mentioned earlier, the distribution of multipolar amplitudes is zero for values 
of j < |/ -l-p|. This second degree of freedom can be also used in the paraxial case, as 
displayed in Figure 4.1, but its influence is less important. Also, comparing Figures 4.1 
and 4.2, one can see that in the paraxial case more multipoles are needed to describe 
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Figure 4.2: Beam shape coefficients (\Cjm*p\‘^) for different focusing scenarios. The insets 
represent the intensity plots of the paraxial LG modes at entrance pupil of the lens. These 
paraxial modes are focused with two lenses of NA= 0.9 and NA= 0.25, and their multipolar 
decomposition given by Cjm%p is found at the focal plane. The red coloured bars indicate 
NA= 0.9, and the yellow ones NA= 0.25. The paraxial beams used are a) LGo, 0 ) b) LG 2 ,o, 
c) LG 4 ,o, and d) LGe.o- Note that Cjm^p can be described with very few multipoles when a 
high NA lens is used. 


the same LG modes, even if a lens with small NA= 0.25 has been nsed. Now, similarly 
to the previons section, once Cjm%p has been calcnlated, the fields in the whole space 
are obtained applying the GLMT, and expressions (4.13) are re-obtained: 
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(4.24) 


































4.4 Excitation of Whispering Gallery Modes 


65 


where Cjm*p is given by expression (4.22). Finally, I want to acknowledge that before 
I published these results in [107], based on the work done by my supervisor in [120], 
N.M. Mojarad et al. had observed the influence of the NA of a lens in focusing a plane 
wave onto plasmonic particles [125, 126]. However, their study only dealt with a single 
focused plane waves, and therefore it could not be applied to any paraxial held. Also, 
they never applied their formalism to more elaborated beams such as vortex beams. 

4.4 Excitation of Whispering Gallery Modes 

In this section, I will look into one of the applications of GLMT - the excitation of 
WGMs, or Morphological-Dependent Resonances (MDRs) [86, 127]. WGM and MDR 
are interchangeable expressions, however the term MDR has been typically only used 
in the context of GLMT, whereas WGM has been used in almost any other branch 
of physics. Hence, I will use the term WGM throughout this thesis to refer to MDR. 
In the next pages, I will use paraxial and non-paraxial cylindrically symmetric beams 
to excite single WGMs in dielectric spheres. I will show that provided an adequate 
ratio of R/X is obtained the excitation beam, a single WGM can be excited with direct 
light, without the need of going to evanescent coupling. In order to do that, I will 
need to compute the scattering efficiency factors defined in section 3.5 for an incident 
field given by expression (4.8) and a scattered and interior fields given by equations 
(4.13). As described in chapter 3, the first step is computing the temporal average of 
the Poynting vector, (S®). Then, the flux of it across a spherical surface is computed 
to obtain the scattering cross section Cg- The general result was given in chapter 3 by 
equation (3.27): 

oo j 

= + \bj9jZl\^) (4-25) 

j=l m^=-j 

Then, if the result is particularized for = F (2j -|- l)^/^Gjm*p and = R^^p{2j + 
(see equations (4.4, 4.8)), the scattering efficiency factor Qg can be computed 

[107, 121]: 

j=|m*| 


(4.26) 
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where m* = Z + p , p being the helicity of the incident beam, and I the topological 
charge of the LG mode used to produce the focused beam; and x = 27rR/X is the 
size parameter of the problem. It is insightful to compare equation (4.26) with the 
expression for Qg in Mie theory. This expression is [46]: 

= + R27) 

J=1 

Expressions (4.26) and (4.27) give rise to very different scattering efficiency as a func¬ 
tion of the wavelength. I showed these differences in [121] using paraxial LG beams 
as E7 That is, I used the equation (4.9) to compute Cjm^p- The differences between 
the scattering efficiency using equations (4.26) and (4.27) are summarised in Figures 
4.3, 4.4, 4.5. Firstly, I use equation (4.27) to compute the scattering efficiency for 
Mie theory. The result is shown in Figure 4.3. It can be observed that high order 
multipoles only give rise to ripples in ths scattering structure. That implies that the 
scattering is dominated by the low order modes. Secondly, I use equation (4.26) to plot 



Figure 4.3: (blue curve, on top) for a single sphere with rir = 1.33 and x G [50, 52]. 

At the bottom, four Mie coefficients ( 055 , 655 , a^Q and 650 ) are plotted. It is seen that the 
ripples in the scattering efficiency are caused by high order Mie resonances. 
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the scattering efficiency for the same incident four LG modes whose multipolar decom¬ 
position has been given in Figure 4.1, i.e. LGq.o, LG 2 ,o, LG4,05 LG6 ,o- As seen in Figure 
4.4, an increase in the AM of the incident beam is translated into an enhancement of 
the ripple structure with respect to the background produced by the low order modes. 
This phenomenon is clear when we compare Figure 4.4(a) with (4.4(d)). In Figure 
4.4(a), it is observed that the scattering is dominated by the low multipolar orders, 
and the high ones only contribute to the scattering as ripples, very similarly to Figure 
4.3®. Nevertheless, in Figure 4.4(d) the ripple structure is largely enhanced and the 
scattering is dominated by the ripple structure itself. In order to highlight this effect 


o 


(/) 

o 





Figure 4.4: Qg as a function of wavelength, A G [450,800] nm, for different incident LG 
beams, a) LGo,o b) LG 2 ,o c) LG 4 _o d) LGe^o- All the rest of parameters do not vary from a) 
tod): nr = 1.5,p = l,i?=1.3 pm. The multipolar decomposition of each of the LG beams 
is given in Figure 4.1. 


®The underlying reason is that both a plane wave and Gaussian beam are very similar when the 
sphere is much smaller than the beam waist of the Gaussian beam wq. Furthemore, both of them 
have the same AM content, Jz = P- 
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and quantify it, the influence of the same four excitation beams (LGo,o, LG 2 ,o, LG 4,05 
LG6 ,o) on a single ripple is plotted in Figure 4.5. The maximum value of the scattering 





Figure 4.5; Variation of a resonance in Qg as a function of wavelength, A G [728.6, 744.6] 
nm, and AM of the incident beam. The excitation beam is a) LGo,o b) LG 2 ,o c) LG 4 ^o d) 
LGg^o- All the rest of parameters do not vary: = 1.5, p = 1, R = 1.3 pm. The multi¬ 

polar decomposition of each of the LG beams is given in Figure 4.1. The chosen resonance 
corresponds to 047 , that is, an electric multipole with j = 17. 


efficiency, Qmax, which peaks at 736.6 nm, and the minimum value are computed for 
each plot. The minimum value is computed as the average of two values 8 nm apart 
in wavelength from the position of the maximum. Then, the ratio 6 can be defined as 
the relative quotient between the subtraction of the maximum and the minimum, and 
the minimum. 
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Qmax Qr, 

Qmin 


(4.28) 
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Table 4.1: Enhancement of the ripple structure using different LG beams. 



LGo,o 

o 

1-1 

O 

o 

1-1 

LG 6,0 

5(%) 

8.184 

61.30 

153.9 

218.9 


where and Q™"" = + QP®'®)/2. The results are presented in table 

(4.1). It can be seen that the background is highly reduced. When a LGo,o is used 
to excite the sphere, the peak only represents 8 % of the background. Nevertheless, 
when the incident beam is a LGg^o, the peak contribution to Qg is two times larger 
than the addition of all the other modes. This is easily understood looking at equation 
(4.26). Due to the rotational symmetry of the system, the AM needs to be preserved. 
That is, if E' is an eigenstate of has with value m*, its multipolar decomposition 
is such that the hrst (m* — 1 ) multipolar modes are not included in the scattering. 
This way of enhancing the ripple structure by cancelling the contribution of the low 
order modes has a severe drawback: it reduces a lot the scattered power. In order to 
make up for it, a higher NA lens can be used. This is shown in Figure 4.6, where the 
same resonance presented in Figure 4.5 is excited with a LG 6 ,o focused with different 
NAs. It can be observed that there is a difference in more than 10 orders of magnitude 
between the scattering produced by a paraxial beam (or NA^i 0.25, see section 4.3) 
and a highly non-paraxial one. Thus, highly non-paraxial beams with a well-dehned 
AM and helicity are a good choice to excite single resonances, and in particular WGMs. 

WGMs are widely used in physics. Their incredibly high Q factors make them very 
useful to probe any sort of disturbance in the environment [ 86 , 128-131]. Their deh- 
nition can be found in a very good review done by Oraevsky [ 86 ]. He dehnes a WGM 
as a single multipolar mode, with j large, = j and ^ = 1 , where ^ is dehned in 
the next lines®. Even though dehning something as ‘large’ is not very mathematical, 

® Oraevsky does not state anything about the parity of the multipolar mode, nor the frequency. 
Thus, a WGM can have both magnetic or electric parity, and can also be excited at any wavelength. 
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Figure 4.6: Effect of the NA of the lens on the excitation of a single resonance at A = 736.6 
nm. A LGg^o with helicity p = 1 is used for all the cases, but different lenses are used to 
focus it down to the particle. The particle has a radius R = 1.3 fim, and a relative index of 
refraction = 1.5. In all the cases, Qs is plotted in logarithmic units. 


it makes sense in this context. The WGM resonances are computed as the roots of the 
Mie coefficients of order j. Each Mie coefficient has multiple roots, which are usually 
listed by the root number The Q factor of the resonance decreases with the root 
number .^ = 1,2,.., therefore given a WGM mode {i.e., given j) the narrowest WGM 
resonance happens for = 1 [86]. Due to losses, it is experimentally very challenging 
to obtain Q factors larger than 10®, therefore a large j can be understood as such that 
makes the Q factor of a resonance larger than 10®. There is not a quantitative answer 
to this, because the Q factor not only depend on j, but it also depends on the relative 
refractive index of the sphere with respect to the surrounding medium, 


(v) 

as long as the other conditions are fulfilled. I will denote them as Ajl'k 
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Now, using beams of light with a well-dehned Jz makes it straightforward to excite 
modes of light with a single m* = j, where j is large. The control of rriz is done 
with the azimuthal number of a beam and its helicity p: m* = I + p. However, 
following Oraevsky’s dehnition, this is not enough to excite a WGM - only one j needs 
to be excited, too. 

To achieve that, the idea is to use the properties of the Mie coefficients for dielec¬ 
tric particles^. In this case, the Mie coefficients, aj and bj, that is the multipolar 
moments of the sphere, are complex and their absolute values are bounded between 
zero and one, 0 < {|aj|, |^j|} < 1- Both aj and bj are very close to zero for small 
X = 27iR/X, and they start to grow for a value of x which is proportional to the order 
of the mode j. The proportionality value is a function f{nr) that depends on For 
example, when Ur = 1.5, it can be computationally verihed that the Mie coefficients 
approximately start to grow when x ~ 4j/5 (see Figure 4.7) [132]. Then, if a: < 4j/5 
both Qj and bj are negligible. For x > 4:j/b their absolute values oscillate between 0 
and 1. This is depicted in Figure 4.7, where aj and bj are plotted for j = 1,10,20. 
It can be observed that ah the Mie coefficients follow the pattern described above: 
their absolute value is 0 until x ~ 4j/5, and for x > 4j/5 they oscillate between 0 
and 1. Hence, it always exists an interval around x ~ 4j/5 where aj and bj start 
growing and the higher Mie coefficients are approximately zero (as they start growing 
for X ~ 4(j -|- l)/5 ). Consequently, if the incident beam has a. Jz = m*® and we pick 
a particle whose size R is such that kR ^ f{nr) ■ ml, only modes with j = ml will 
be excited. This condition can be achieved by tuning the wavelength of the incident 
beam, A = 271/k. In summary, a single WGM with indexes j and = j (A^-^^) can 
be excited by doing the following two steps process: 

1. We create a cylindricahy symmetric beam with a large value of Jz = m* = I + p 
and we select a MO with a large NA to tightly focus it. The wavelength of the 

^Remember that a dielectric particle has a real index of refraction. 

® As mentioned in chapter 2, this is an abuse of language. It means that the beam is an eigenstate 


of Jz with eigenvalue ml. 
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Figure 4.7: Norm of the Mie coefficients \aj\ (in red) and \bj\ (in blue) for j = 1,10, 20 in 
a), b) and c) respectively, as a function of the size parameter x = ‘I'kt jX . The relative index 
of refraction is = 1.5. It can be seen that all of them start being significantly different 
from zero when x ~ 4j/5. 

laser is supposed fixed. 

2. We choose a sphere whose radius R is such that the following condition is fulhlled: 
X - ml ■ f{nr) / (27r) 

Next, I will illustrate how precise and flexible this method of exciting WGMs is by 
giving three examples of excitation of WGMs. In the three examples, A will be hxed 
at A = 532 nm: 

• Example 1 : I want to excite a WGM with j = 15, Hence, a beam 

of light with ml = 15 needs to be used. That can be a LGi 4 ^o beam with a 
helicity p = 1. Now, the size of the sphere needs to be selected so that only 
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a single multipolar mode is excited. All the spheres are surrounded by air and 
made of a material with ni = rir = 1.5. It can be computationally proved that 
Ur = 1.5 —)■ /{rir) ~ 0.8 [132]. As I want to excite a mode with j = 15, then 
I need that kR/f{nr) ~ 15 —)■ kR ~ 12. Given the excitation wavelength of 
A = 532 nm, the condition is fulfilled by a particle with i? = 1.0 pm. It can be 
checked that this multipolar resonance has a magnetic parity and a Q factor of 
Q = 70. 

• Example 2 : A WGM with j = 40 needs to be excited and the available spheres 
and surrounding medium are still such that = 1-5. A LGsg^o beam with p = 1 
is created to excite the resonance. Then, the radius of the sphere needs to be 
such that kR/firir) ~ 40 —)■ kR ^ 31, which is fulfilled hy R = 2.6 pm. In this 
case, the Q factor is Q = 5 ■ 10^, and the parity is still magnetic. 

• Example 3 : The material which the sphere is made of is changed and now 
Ur = 2. It can be computationally checked that this implies that f{nr) ~ 0.6. 
I want to excite the same multipolar modes as in the previous case, i.e. 
Therefore, I need to used a LGag^o mode with p = 1. The condition on the radius 
of the sphere is the same one, but the result is different due to the change of 
material: kR//(rir) ~ 40 —?■ kR cs 23, which yields R = 1.9 pm and a Q factor 
Q = 1 ■ 10®. 

WGMs or multipolar modes are normal modes of the sphere, and as such they are 
symmetric under parity transformations (see sections 2.4 and 3.2). It is interesting to 
see what happens when we decompose the scattering into its two orthogonal helicity 
components. For an electric or magnetic resonance, the scattering is split equally into 
two orthogonal helicity components [54, 107]. This effect is put forward in Figures 4.8, 
4.9. 

First, Figure 4.8 shows the background-suppression effect previously shown in Fig¬ 
ures 4.4, 4.5. That is, in Figure 4.8(a) I show the scattering efficiency produced by 
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(a) (b) 



k (nm) 



Figure 4.8: Suppression of background in the scattering efficiency (Qs)- Input beams 
are a) LGo,o and b) LGi 8 , 0 ) both with p = 1. In a) an aplanatic lens with NA= 0.25 is 
used, whereas in b) an aplanatic lens with a NA=0.9 is used. The rest of parameters are 
kept constant for both plots, i.e. R = 1.3 pm, rir = 1.5. The insets indicate the typical 
profile of the beam used to excite the sphere. The scattering efficiency is represented with 
a blue continuous line. The Mie coefficients 620 and 022 are plotted with a red and a green 
line respectively. In a) we have indicated with a dashed line the position of two particular 
resonances for these modes. Note that the ordinate axis in a) and b) are different. The y-axis 
scale in both plots is given by Qs- Then, 620 and 022 are re-scaled to match the values of Qs- 


a Gaussian beam on a particle whose radius is i? = 1.3 pm and whose relative index 
of refraction is = 1.5. Two resonances given by 022 and 620 are singled out. Their 
contribution to the scattering is seen to be a single ripple, due to the large background 
produced by all the other modes. However, in (4.8(b)), where a focused LGig^o is used, 
the scattering efficiency is completely dominated by the ripple structure, as the back¬ 
ground produced by the low order multipolar modes is not excited due to the large 
value of Jz- 


Then, in Figure 4.9, the intensity of the electric held (|Ep) in all the space is plotted 
at the 620 resonant condition, A = 503 nm, and 2 nm out of resonances at A = 505 nm. 
The particle and the surrounding medium are still the same, i.e. R = 1.3 pm, Ur = 1.5. 
|Ep is plotted for two excitation cases: E^ = LGo.odp and E‘ = LGig^odp, both of them 
with p = 1. Then, E‘ = LGo,odp is focused with a NA= 0.25 and E^ = LGi 8 ,odp 
with a NA= 0.9. The resonant behaviour of the scattered held is especially clear when 
the opposite component of the helicity A“ is examined, since the contribution due 
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to the incident field is 0. Figures 4.9(a)-(d) depict the fact that the influence of the 
620 resonance on the scattering is minimum, as the same behaviour is observed for on 
and off-resonance cases. However, in Figures 4.9(e)-(h) the situation is very different. 
There are two orders of magnitude of difference in intensity between the cross helicity 
component in Figure 4.9(f) and 4.9(h). Finally, it can be observed that unlike Figure 
4.9(e-h), where the field can be described with a single mode, many multipoles are 
needed to describe the highly complex pattern seen in Figure 4.9(a-d). In fact, the 
intensity profile observed in Figure 4.9(a-d) has drawn a lot of attention since 2004, 
where it was first characterized as photonic nanojet [133]. Since then, nanojets have 
been extensively studied and applied in many different fields [134-137]. 

Previous works on the excitation of metallic particles by LG beams were carried out 
in [138]. There, the use of purely numerical techniques allowed the authors to distin¬ 
guish certain properties of the scattered field. However, due to the large absorption of 
silver, the ripple structure could not be enhanced in a similar fashion to what I have 
shown. On the contrary, the control of the AM of LG modes was used to achieve optical 
transparency in dielectric particles in [139]. LG modes were also used in some recent 
scattering experiments with silica spheres [140]. The theory developed here allows for 
a better understanding of the experiments. In fact, the same calculations could also 
help to expand the computations of the optical forces when the incident beam is not 
a plane wave [141] and when the particle is not much smaller than the wavelength 
[142-144]. Last but not least, the technique presented in the present chapter to control 
the multipolar content of cylindrically symmetric beams, as well as its application to 
excite single multipolar modes could open up new possibilities in the forefront research 
of magnetic resonances [145-149]. 
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Figure 4.9: Projection into the helicity basis of the intensity of the interior and total field 
from a sphere with parameters {7? = 1.3 /rm, = 1.5}. The helicity of the incident beam 
is always p = 1. The intensity is plotted in logarithmic scale, where 0 corresponds to the 
maximum of intensity in a). Images a), b), e) and f) have been simulated for a Aresonance = 503 
nm, whereas a 2 nm shift has been introduced in the off-resonance wavelength for the others, 
i.e. Aoff-resonance = 505 nm. The simulations for the LGo,o beam have been executed with a 
NA=0.25 lens in order to be able to excite the resonance of order 20. Note that the scattering 
is dominated by the off-resonance modes, since a change of 2nm in A leaves the field invariant. 
Calculations with higher numerical apertures give similar results. On the contrary, a lens with 
a NA=0.9 has been used for the excitation with a LGops- In both cases, the entrance pupil 
of the lens was filled. Note that in contrast with b), the mode of order 20 is greatly enhanced 
when we tune the wavelength to the resonance. 


























“Physics is iike sex: sure, it may give some practicai resuits, but 
that's not why we do it" 


Richard P. Feynman 
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5.1 GLMT problem in helicity basis 

In the previous chapter, I have used the GLMT developed in chapter 3 to find the scat¬ 
tering off a sphere illuminated with a cylindrically symmetric beam with a well-defined 
helicity. As explained in chapter 3, the problem is described with the multipolar fields 
since they are the normal modes of the system. Consequently, and as we 
have seen in section 4.4, the resonances of the system have a well-defined parity (either 
electric or magnetic). In the current chapter, I will show that a wealth of information 
about the sample can be easily obtained by changing the basis in which the problem 
is described. The same information could be obtained with the expressions derived 
in chapter 3, but the change of basis will unveil this information in a straightforward 
manner. I will describe the problem with the multipolar fields with a well-defined 
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helicity (see section 2.4). Their relation with the multipolar helds (with well-dehned 
parity)^ is the following one [60, 150]: 


A"*" = 

jrriz 


■^jrriz ' ^'^jmz 


= 


jrriz ^ jrriz 


V2 

It can be checked that these modes have a well-dehned helicity, i.e. they are eigenstates 
of the A operator. In order to check that, the following relations between the multipolar 
helds need to be taken into account [51, 107]: 


AA'S = 


(e) 


A A = -?■ A 

jrriz ^-^jm 


(m) 


(6.2) 


Then, it is straightforward to check that: 


- +-^j7n. 




(5.3) 


As discussed in chapter 3, duality transformations are generated by the helicity operator 
[63, 64]. It is well-known that Maxwell equations are not symmetric with respect 
to electric and magnetic helds in the presence of charges [42]. This is due to the 
lack of magnetic monopoles in the universe [151]. However, it was proven in [41] 
that the macroscopic Maxwell equations for isotropic and homogeneous media can 
be dual-symmetric if some conditions are fulhlled. Microscopically, duality symmetry 
is still broken, but the collective ehect of all the charges and currents in the medium 
restores the symmetry in the macroscopic approximation. In [41], the non-conservation 
of helicity was carefully quantihed. In the same way as it happens with any other 
generator of a symmetry, if the helicity of a light beam is preserved upon interaction 
with a material medium, this necessarily implies that the system is symmetric under 
its associated duality symmetry. I will refer to these sort of media as ‘dual’, i.e. media 
that preserve the helicity content of the incident beam upon scattering. In order to 
study the dual behaviour of spheres, I will use beams of light with a well-dehned helicity 
with value p. Its decomposition into multipolar modes of well-dehned helicity will be 
the following one: 

= (5T 

3,'nnz 

^Due to their common use, the multipolar fields with well-defined parity will be simply referred to 


as multipolar fields. 
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where amplitudes that determine the multipolar content of the incident 

beam^. Then, the scattered and interior held can be obtained following the procedure 
shown in section 3.3. The following helds are obtained: 


j^int 


j,mz 

/ y jrriz 


bj + a 


2 jm^ ' 
Cj + dj p 
2 





(5.5) 


J^rriz 

where {aj, bj, cj, dj} are the Mie coefficients [46]. The rest of quantities dehned for the 
GLMT in chapter 3 could be equally derived, but I will not compute them as I will not 
need them for my purposes. 


5.2 Kerker conditions 

In 1983, Kerker et al. demonstrated that a sphere with /x 7 ^ 1 could have zero scattering 
in both back and forward directions [152]. These two anomalous conditions, nowadays 
known as hrst and second Kerker conditions, arguably passed unnoticed until 2006, 
when Mehta and co-workers published experimental evidence of zero forward and back- 
scattering by magnetic spheres [153, 154]. Since then, a lot of work has been done in 
this held, as researchers seek to independently control both the electric and magnetic 
resonances of diherent structures and control its directionality [145, 146, 155, 156]. Re¬ 
cently, diherent research groups have managed to measure the two Kerker conditions 
in diherent regimes. Gehrin and co-workers measured these two conditions in the mi¬ 
crowave regime [157]. Whereas Person et al. and Fu et al. measured the hrst Kerker 
condition in the dipolar approximation for optical wavelengths [158, 159]. Nonethe¬ 
less, the measurement of the Kerker conditions in the optical regime for arbitrary large 
particles still remains an unsolved problem. 

In this chapter, I will extend the Kerker conditions to systems with cylindrical sym¬ 
metry. We will see that this extension naturally arises from equations (5.5). Indeed, 

will not call the coefficients beam shape coefficients as I reserve that name for the coefficients 

(m) (e) 

and that multiply the multipolar fields in the GLMT described in chapter 3 
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looking at equations (5.5), it can be seen that the scattered and interior helds do not 
generally preserve the helicity of the incident held (see Figure 5.1(a) and 5.1(b)). That 
is, in general, an incident beam with a well-dehned helicity such as E' impinging on an 
isotropic and homogeneous sphere will give rise to scattered and interior helds which 
will not be eigenvectors of the helicity operator A. This is evident from equations 
(5.5), as multipolar modes with helicity —p (AJ^J are created upon scattering. The 
underlying reason for this is that the two pairs of Mie coefficients {aj, bj} and {cj, dj} 
are not generally equal and therefore the amplitudes modulating AJ^^ will be diherent 
from zero. 

However, as it was proven by Kerker et al, the Mie coefficients can be equal for some 
cases, i.e. aj{x) = bj{x) Vj,x when = Pr [152]. This is the so-called hrst Kerker 
condition (Kl), and it implies zero-backscattering from the sphere in consideration. 
Actually, observing equations (5.5), it can be seen that the Kl also implies that the 
sphere in consideration will preserve the helicity of the incident field upon scatter¬ 
ing and therefore be dual (see section 5.1 to refresh the dehnition of a dual system). 
Indeed, for the scattered held, the amplitudes of the components with opposite helicity 
to the incident one are given by a^^^{bj — aj)/2. Hence, when aj = bj they are all zero, 
granting the scattered held with the same helicity as the incident one. This feature 
has been overlooked in the past, and it will enable me to extend Kl. 

The Kl is a particular case of a more general condition that restores the EM dual¬ 
ity symmetry in material media [41]. In general, 

— = const A conservation (5.6) 

A* 

for any medium made of an arbitrary number of isotropic and homogeneous sub-media 
i, Cj and pi being its electric permittivity and magnetic permeability. That is, when 
an arbitrary light beam impinges on a medium such that the electric and magnetic 
helds behave symmetrically, i.e. condition (5.6) is met, the helicity of this beam is 
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not changed regardless of the geometry of this medium. Now, even though the zero- 
backscattering condition or K1 was introduced using spheres, it can be demonstrated 
that it can be relaxed to any system with cylindrical symmetry [150]. Suppose we have 
a dual scatterer and impose that it must also be cylindrically symmetric around the z 
axis. Then, besides helicity, the Jz of the incident beam must also be preserved. Under 
these conditions, a plane wave incident along the symmetry axis will not backscatter. 
The proof comes from the dehnition of helicity A = J ■ P/|P| (see section 2.2). Imagine 
that a plane wave is directed in the ^ axis (its linear momentum is P = P^z) and it has 
a well dehned helicity A = Jz ■ Pz/\Pz\ = Jz = 1- A plane wave with P = — P^z must 




Figure 5.1: Sketch of relations between Kerker conditions and cylindrical and Duality 
symmetry. The large red helix pointing the scatterer represents an incident plane wave with 
a well defined helicity. The other helices represent scattered plane waves. The small red 
helix has the same helicity as the incoming wave, while the small blue helix has the opposite 
helicity. a) A general scattering process with neither cylindrical nor dual symmetries is 
shown. In general, both back and forward scattering have contributions from both helicities. 
b) Scattering from a cylindrically symmetric object. In this case, symmetry imposes that 
forward scattering is of the same helicity than the incident held, and backward scattering 
of opposite helicity. c) A dual and cylindrically symmetric system. Only forward scattering 
is allowed with the same helicity. d) An anti-dual system with cylindrical symmetry. Only 
back scattering is allowed with hipped helicity. 

necessarily have K = — Jz- If the system preserves Jz and A, such plane wave can never 
exist. Thus, it is very natural to extend K1 from a sphere such that er = to a dual 
system with cylindrical symmetry, i.e. a system that preserves both helicity and the 
z component of the AM (see Figure 5.1(c)), as both entities have zero backscattering 
[41, 150]. 
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The second Kerker condition (K2) can also be extended in a similar fashion. K2 pre¬ 
dicts zero forward scattering for dipolar magnetic particles. Here, I will relax its nse to 
cylindrically symmetric systems of any size. In order to do this, I am going to define 
what an “anti-dnal” scatterer is. It is a scatterer whose scattered light flips the helicity 
state of the incident beam. That is, if E' = ^ then E" = 

Looking at eqnations (5.5), the conditions for an anti-dnal sphere can be fonnd. For 
E®, the amplitndes of the components with the same helicity as the incident held are 
given by Thus, a necessary and sufficient condition for a sphere to be 

anti-dual is that: 

aj{x) =—bj{x) Wj,x (5.7) 

If this condition is expressed in the dipolar approximation, i. e only terms with j = 1 
are considered, it yields Oi = —bi. This is the so-called K2. 

With this new perspective, the zero-forward scattering condition for spheres (K2) can 
be cast again as a particular case of a more general system, a system which is cylindri¬ 
cally symmetric and anti-dual. Once again, only considering the symmetries of such a 
system, it is possible to derive that the forward scattering must be zero. To prove that, 
a plane wave travelling in the positive z axis, whose helicity is A = Jz ■ Pz/\Pz\ = 1 
is used as Eh Now, as the system is anti-dual, the helicity in all the directions (in 
particular in forward) have to flip, giving A = —1. However, cannot change since 
the system is cylindrically symmetric, and Pz/\Pz\ = 1 in the forward direction. Con¬ 
sequently, this plane wave cannot exist (see Figure 5.1(d)). 

At this point, it is important to note that hnding an exact anti-dual sphere is chal¬ 
lenging. Some authors have noticed that dielectric, isotropic and homogeneous spheres 
cannot behave as anti-dual materials in the strict sense, as they would contradict the 
Optical Theorem^ [156, 160, 161]. But a possible alternative would be to use particles 

^Their argument is only valid for plane wave excitation, though. As it is explained in [48], the 
optical theorem only applies to Mie theory, and fails to provide meaningful results in GLMT, where 
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with gain, as the Mie coefficients can get valnes greater than 1 [132, 150, 156]. 
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Eqnations (5.4, 5.5) give us the expression of the electric field in all space. Then, the 
magnetic field can be obtained using the Maxwell equations given by equations (2.7). 
The following expressions for the H field are obtained: 

(5.8) 

(5.9) 

(5.10) 
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The energy density of the EM held can be computed as [42, 162]: 

1 


w = 


167r 


e|E|2 + /i|H|2) 


(5.11) 


Then, if w is integrated on a surface, the energy of the EM held per unit of length is 
obtained: 

W= / — (e|E|2 + /i|Hndfl (5.12) 

Jn IStt 

To hnd out more information about the scattering process, I will particularize equation 
(5.12) for the scattered held, as it contains all the information about the scatterer. 
Given equations (5.5, 5.9), the orthogonality relations given in section 2.3, the dehnition 
of multipoles with well-dehned helicity (5.1), the following relations yield: 



1 P 

\(y- 

\ jrriz 

r ( n 

(5.13) 

j,mz 




= -V 

> 4 

1 P 
\Oi- 
\ jruz 


(5.14) 

j,'mz 




= -V 

-P 4 

1 P 

Q: • 

\ jmz 


(5.15) 


j,mz 


more complicated beam of lights can be used. 
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represents the total energy of the scattered held. Due to the orthogonality of the 
multipolar helds, the scattered energy can be split into two orthogonal components, 
and these are W® and W^p. They account for the energy scattered in modes with the 
same helicity as the incident held (Wp and with the opposite helicity (Wip). 

As described in the previous section, note that in general the scattered held will al¬ 
ways carry energy in modes with the opposite helicity, he. W^p 7 ^ 0. However, if 
aj{x) = bj(x) Wj,x, the particle only scatters energy with the same helicity as the 
incident beam, conserving the helicity of the EM held and therefore behaving as a dual 
medium. Note also that if aj{x) = —bj{x) Wj,x, then W® = 0. That is, the scattered 
held has the opposite helicity to the incident one and therefore the sphere behaves as 
an anti-dual scatterer. It can be proven that a spherical particle will only be dual if 
fir = Cr [41, 152]. Nonetheless, here I will relax this dehnition in order to address some 
experimentally relevant cases [163-165]. I will consider that a particle of radius R 
behaves in a dual or anti-dual manner at a given wavelength if all the scattered energy 
goes into the same helicity component or the opposite one at that wavelength. 


These two facts have been experimentally verihed in the dipolar regime [157-159], 
experimentally achieving oi = bi and ai ~ —bi. I will show how to extend these con¬ 
ditions to other regimes and shapes. As a result, it will be clear that even a dielectric 
sphere can approximately fulhl the anti-dual condition or extended K2. In order to 
achieve these purposes, I will dehne a transfer function Tp{R,nr)- Its dehnition de¬ 
pends on the helicity value of Eh the radius R of the sphere in consideration and the 
relative refractive index rir'. 
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(5.16) 


Tp{R, Ur) is the ratio between the scattered light going into modes with opposite helicity 
with respect to the incident light (Wi^), and the scattered light going into modes with 
the same helicity (Wp, for a given helicity (p) of the incident beam. Thus, Tp{R,nr) 
varies from 0 to inhnity. When Tp{R,nr) tends to zero, the particle is dual and all 
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Figure 5.2: Plot of the log of the Transfer function Tp{R,nr) for p = 1, i.e. 

log {Ti{R, Ur)), as a function of the radius R of the particle (horizontal axis) and the rel¬ 
ative index of refraction (vertical axis). The radius of the particle is varied from lOOnm 
(left) to 300nm (right) and the relative index of refraction goes from 1.2 < < 3. 


the scattered light has helicity p: in other words, it fulhls the extended K1 [132, 150]. 
On the contrary, when Tp{R,nr) tends to inhnity, the particle is anti-dual and all the 
scattering is transferred to the cross helicity, —p: it fulhls the extended K2 [132, 150]. 

Given an incident held E' and a sphere with radius R and relative refractive index 
Tir, Tp{R,nr) can be computed and the combinations of {A,i?, that make the par¬ 
ticle behave as dual or anti-dual can be studied [132]. For example. Figure 5.2 depicts 
the transfer function of sphere excited with a Gaussian beam at A = 780 nm. It is seen 
that the sphere can behave as dual for many resonant combinations of {R,nr}- This 
way of looking at duality behaviour and Kerker conditions represents a considerable 
improvement over the methods used before, where only the hrst order Mie coefhciens 
were used [145, 155, 156]. As it can be observed in Figure 5.2, not all the points sitting 
on the curve Oi = bi are equally dual. Generally speaking, it is due to the fact that the 
particle has a size (or an index of refraction) for which the dipolar approximation does 
no longer apply. The transfer function dehned in (5.16) exactly quantihes this fact, as 
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Figure 5.3: Plot of |ai — 6i| as a function of the radius of the sphere R (horizontal axis) 
and the relative index of refraction (vertical axis). It can be observed that there are three 
major regions where |ai — 6i| = 0. The wavelength is set to A = 780nm. 


it accounts for the contributions of all the higher multipolar orders present in scattering. 

It is also interesting to see how the induced duality strongly depends on the rela¬ 
tive index of refraction. It is apparent from Figure 5.2 that the dielectric sphere gets 
closer to the dual condition when gets larger. That means that the helicity of the 
incoming beam will be better preserved when particles with a high refractive index 
embedded in a low refractive index medium are used. 

In addition to the duality considerations mentioned above, Figure 5.2 also depicts 
the anti-dual behaviour of the particle. Two main features can be observed. First, as 
it can be deduced from the colorbar, the anti-dual condition is not achieved as hnely as 
the dual is. That is, the extension of K2 is more difficult to achieve. Hence, the forward 
scattering is never reduced as much as the backward is. This feature had already been 
observed and analysed in [160]. Moreover, it is consistent with the few experiments 
done until now [157, 159]. Secondly, if the relative index of refraction is maintained. 
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the anti-dual condition is held for larger particles than the dual one [145, 166]. 


It is worth mentioning that the same method for characterizing the dual behaviour 
of a medium can be undertaken even though its geometry is not spherical. In that 
case, the analytical expression of the transfer function cannot be given, but the idea is 
the same one. The scattered energy can still be split into two orthogonal components 
with a well-defined helicity content, i.e. W® = W® -|- Then, the only difference lays 
in the calculation of W® and WL^. Whereas this calculation for a sphere only depends 
on the Mie coefficients, for a general medium is achieved in the following way [41]: 

w;= ||E® + ^pH®||2 (5.17) 


where the norm is computed as a surface integral. Therefore, the equation for the 
transfer function in general is: 


Tp(A, Ur) 


E® — ipH®||^ 

E" 


(5.18) 


It can be checked that equations (5.14, 5.15) are re-obtained if the expressions (5.5, 
5.9) for E® and H® are introduced in equation (5.17). For a more general case though, 
numerical methods are needed to compute E® and H®. This method has been applied 
both in [41] and in [163] to compute the helicity content of the scattered field. As 
mentioned in chapter 3, the calculation of E® and H® depends on E^ but an eigenmode 
calculation could retrieve calculations which are independent of the excitation beam. 


5.4 Angular momentum-induced helicity transfer 

In this section, I will particularise the equations (5.4, 5.5) for the case of E' being 
cylindrically symmetric (and with a well-dehned helicity). As in the previous chapter, 
I will consider that its eigenvalue for is m* = / -f p, where I is the topological charge 
of the beam in the paraxial approximation. Thus, equations (5.4, 5.5) will be simplified 
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Figure 5.4: Plot of log (|a 5 — 65 !)) as a function of the radius R of the particle (horizontal 
axis) and the relative index of refraction rir (vertical axis). 
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Then, the transfer fnnction can be re-written as: 


, , Mi, - aj\ 

w- ^ Er-wi KPh-+“t 


(6.19) 

(5.20) 

(5.21) 


(5.22) 


Given an incident cylindrically symmetric beam, I will show how the dual and anti-dual 
properties of the particle previously depicted in Figure 5.2 dramatically change when a 
higher angular momentum mode is used as Eh To exemplify this, a cylindrically sym¬ 
metric beam with m* = 5 and p = 1 is used. In particular, the beam is a left circularly 
polarized LG 4 ,o focused with a lens whose NA= 0.9. Its decomposition into multipoles 
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Figure 5.5: Plot of the log of the transfer function Tm*p{r,nr) for rriz = 5 and p = 1, 
log Ur)), as a function of R and Ur- 


is given by Figure 4.2 in section 4.3. If we compare equations (5.19) and (4.22), it 
can be observed that the relation between the beam shape coefficients Cjm%p and 
is R{2j + lY/^Cj^*pV2 = Then, equation (5.22) can be used to compute the 

transfer function for this new excitation scenario. 

In Figure 5.4, the shape of log (jos — & 5 I) as a function of R and Ur is displayed (A = 780 
nm), where the log function has been used to make the plot more readable. Although 
Figure 5.4 has a similar shape to Figure 5.3, where |ai — 61 1 was plotted, some dif¬ 
ferences can be found. The range of sizes for the particles to achieve the condition 
as = 65 has increased for the same interval of rir = [1.2, 3]. Before, it spanned 200nm, 
whereas now it spans 500nm. Nonetheless, as previously discussed, this plot does not 
give any information regarding the higher multipolar orders. Therefore, I have plotted 
log (T 5 i(i?, Rr)) in Figure 5.5 to capture this behaviour. 

As in the j = 1 case, an increase in is linked to an increase of the dual proper¬ 
ties of the sphere. Also, it can be inferred from a comparison between Figure 5.5 and 
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Figure 5.2, that the dual and anti-dual conditions are fulhlled with a better approxi¬ 
mation in this new occasion. Indeed, the minimum value of the colorbar drops almost 
an order of magnitude more, and the maximum value increases more than two orders 
of magnitude. That is, for some certain combinations of R and the energy of the 
scattered held in the modes of opposite helicity is 10, 000 times smaller than the en¬ 
ergy going to the original helicity of the incident held; whereas for some other certain 
conditions, the scattered held energy is dominated by modes with the opposite helicity 
with a ratio of 500 to 1. As in the incident Gaussian beam case, an increase in Ur is 
linked to an increase of the dual and anti-dual properties of the sphere. To summarize, 
some very general conclusions can be reached after a careful look into these results: 

• The larger the relative refraction index of the particle Ur is, the more accurately 
the two extended Kerker conditions shown in [150] can be achieved. 

• Fixing the index of refraction and size of the particle, we can always approxi¬ 
mately induce the duality symmetry by choosing a right combination of AM and 
optical frequency of our laser, regardless of how big the sphere is with respect to 
the wavelength. 

Even though I have not included other plots for other excitation beams, I have con- 
hrmed these conclusions by doing calculations similar to the ones presented with dif¬ 
ferent sizes, wavelengths, refractive indexes, and excitation beams. The results have 
always been consistent with the conclusion above. Finally, note that in these simula¬ 
tions Ur < 3 and A = 780nm. With these two conditions, the smallest a particle can 
be to induce duality symmetry is 120nm. This size could be reduced down to 81nm if 
a A = 532nm was used. 


5.5 Possible experimental implementations 

In the previous section, I have compared different scenarios where the EM duality sym¬ 
metry can be induced with dielectric spheres and arbitrarily high AM modes. I have 
shown that it is possible to induce dual and anti-dual behaviours for certain incident 
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beams E‘, regardless of the nature of the particle. In this section, two hypothetical 
experiments in the laboratory will be quantihed where the duality condition can be 
achieved. 

I will consider spheres made of Silica and Alumina. Their respective refractive in¬ 
dexes at A = 780nm are nsi 02 = 1-54 [167] and nApoa = [168]. I will suppose that 

they are embedded in water, therefore their respective relative refractive index will be 
= 1.16 and = 1.32. The necessary equipment in order to induce this dual 

behaviour will be a tunable laser, a set of wave-plates and linear polarisers to create 
well-dehned helicity states [169], an element to modify the spatial properties of light 
such as a Spatial Light Modulator, an imaging system (two microscope objectives), 
and a dielectric sphere embedded in a homogeneous, lossless and isotropic medium. 
Then, the way to proceed to induce helicity conservation will be the following. I con¬ 
sider the particle and its embedding medium as a given system. It will be clear that 
the helicity of light can be conserved regardless of the nature of the particle (size and 
index of refraction), as long as it is approximately spherical, the surrounding medium 
is homogeneous, lossless and isotropic, and considering a tunable laser with a broad 
enough wavelength modulation. 

Once the two parameters dehning the particle R and are known, the range of 
size parameters x that could be achieved with a tunable laser can always be com¬ 
puted. Supposing that the tunable laser can offer wavelengths spanning from 700nm 
to lOOOnm (that would be the case of a TiiSapphire laser, for example), the achievable 
X will belong to the interval {6.28r(/rm), 8.98r(/im)}. Now, as it has been proven in 
the previous section, there also exist a large number of radii R* (and consequently, x*) 
for which the dual condition is achieved. This is a consequence of the fact that given rir, 
there exist different radii of particles R that make the particle dual provided adequate 
excitation beams with angular momentum are used. Hence it is highly probable 
that regardless of R and rir the dual condition Tm*p{R,nr) ~ 0 can be achieved. In 
fact, this statement is even more true inasmuch as R gets bigger. To be more specihc. 
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Table 5.1: Wavelengths at which the dual condition is achieved depending on the AM 
of the incident beam. The bold wavelengths are those at which the dual condition could be 
achieved with the range of wavelengths available in a Ti:Sap laser. The dual conditions are 
achieved with a minimum precision of Tm^p{r,nr) = 2% for the four different cases. 


A (nm) 

mz = l 

TIlz = 5 

^Si02 

859 

330 

C)Si02 

K 2 

1860 

710 

RAI 2 O 3 

986 

377 

RAI 2 O 3 

2122 

812 


suppose that we have four different spheres. Two of them are made of Silica and the 
other two are made of Alumina. For each of the materials, suppose that the radii of the 
spheres are Ri = 325nm and i ?2 = 700nm. Now, given these sizes and their respective 
index of refraction, the R dependence on the horizontal axis of Figure 5.2 and Figure 
5.5 can be turned into a. x = 2 ttR/\ dependence and obtain the wavelength for which 
the dual condition will be achieved. 

The results are presented in Table 5.1 for the different four combinations of mate¬ 
rials and sizes and for the two different excitation beams 

considered in sections 5.3 and 5.4, i.e. m* = 1 and m* = 5. It can be observed that 
regardless of the size and the material, the duality condition can be achieved if a wave¬ 
length is chosen. Moreover, it is seen that when the particle is larger, higher order 
beams are needed to reach the duality condition, for a hxed A range. Finally, it can 
also be observed that when the relative refractive index is increased, the duality 
condition is pushed to longer wavelengths. 

All this evidence, enables me to conclude that dual systems are easily realisable in 
the laboratory if an arbitrary dielectric sphere is properly matched with a proper light 
beam. The results for aggregates of particles or metamaterials are different, though. In 
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those cases, each of the individual cells is typically smaller than the wavelength, there¬ 
fore the results about high order AM modes do not apply there, as their behaviour is 
dipolar due to their size. However, Figure 5.2 can still be used to engineer some dual 
individual cells [164, 170]. 
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“It doesn't matter how beautiful your theory is, it doesn't matter 
how smart you are. If it doesn't agree with experiment, it's wrong." 

Richard P. Feynman 


6 

Experimental Techniques 


6.1 Holography principles 

The aim of this chapter is to introduce the reader to different experimental techniques 
that will be used to carry out the experiments in chapters 7 and 8. As shown in 
all the previous chapters, the precise characterization of both the spatial properties 
and polarization of light is crucial to control light-matter interactions at the nano¬ 
scale. In particular, a great emphasis has been put into using cylindrically symmetric 
beams, i.e beams of light which are eigenstates of J^- Here, it will be shown how 
to experimentally create these beams in the laboratory with two different techniques: 
Computer Generated Holograms (CGH) and Spatial Light Modulators (SLM). Both 
techniques are based on holography, whose basis are grounded on Fourier optics [45, 
171]. In Fourier optics, the polarization does not play any role - the whole description of 
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light phenomena can be done with a scalar magnitude named complex wave-function 
U{r,t). Furthermore, as it has been done in all the previous chapters, a harmonic 
time dependence exp{—iut) will be considered, therefore U{r,t) = U{r) exp{—iut), 
where U{r) is referred to as complex amplitude. Then, U{r) can be retrieved with the 
intensity /(r) = \U{r)\‘^ and the phase arg{t/(r)} of the wave: 

U{r) = A//(r)exp (zarg{f/(r)}) (6.1) 


Only paraxial waves will be considered here [53]. Paraxial waves were described in 
section 2.2. For these waves, the complex amplitude can be expressed as U{r) = 
A(r) exp{ikz), where A(r) is the envelope of the wave which varies slowly with z. Then, 
a very important magnitude in wave optics is the complex amplitude transmittance of 
an optical element, t{x, y). This magnitude is dehned as the ratio between the complex 
amplitude of the wave after and before the optical element. That is, 

U{x,y,Zf) 


t{x,y) = 


U{x,y,Zi) 


( 6 . 2 ) 


where Zi and Zf are the two z planes where the optical element begins and hnishes. 
In Figure 6.1(a), a sketch of a thin layer of glass is depicted, with Zf — Zi = do- The 
excitation wave is paraxial and propagating in a direction normal to the interface of 
the material. In this case, t{x,y) is naturally dehned as U{r) at the outer face of the 
layer over U{r) at the inner face. Then, if the refractive index of the layer is n, 


U{x, y, Zf) ~ exp [—inkdo] U{x, y, Zi) t{x, y) ~ exp [—inkdo] (6.3) 

where exp(—mfcdo) is fhe phase shift introduced by the optical element^. However, 
when the optical element is not plane-parallel, the dehnition of t{x, y) is a bit more 
complex. If the element lies between Zi and zj, we can dehne a function d{x,y) which 
gives the thickness of the optical element as a function of its transverse coordinates 
{x,y). The air occupies the rest of it, i.e. the thickness of the air is do — d{x,y) (see 
Figure 6.1(b)). Considering the result given by equation (6.3), the transmittance of 
this variable-thickness optical element is: 


t{x, y) ~ exp [—inkd{x, y)] exp [—ik{do — d{x, y))] ~ ho exp [—i{n — l)kd{x, y)] (6.4) 


^Here, reflection at the interfaces and absorption in the material are neglected. 
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Figure 6.1: Sketch of two different transmission plates, a) Thin layer of glass with 
constant thickness Zf — Zi = do- b) Thin layer of glass with variable thickness d{x,y). Both 
figures have been copied from [45] with the permission of Prof. M.C. Teich. 


with ho = exp [—ikd^ being a constant phase factor. Now, a hologram is an optical 
element (a transparency) that contains a coded record of an optical wave. That is, it 
contains the amplitude and the phase of the wave that has been recorded, allowing for 
a complete reconstruction of the held. For example, suppose that a transparency was 
made with t{x,y) = U{x,y), with U{x,y) being the wave to be reconstructed. If this 
transparency was illuminated with a wave U{x,y,Zi) = 1, then the transmitted wave 
would be U{x,y,Zf) = t{x,y) ■ 1 = U{x,y). However, it is challenging to encode both 
amplitude and phase on a transparency. Next, a method to reconstruct t{x,y) only 
modulating the intensity of the beam is explained. Then, in sections 6.2 and 6.3, two 
different techniques to carry out phase modulation holograms will be discussed. 


6.1.1 Intensity modulation holograms 

There are different techniques to make intensity (or amplitude) modulation holograms. 
Since the hrst holograms made by Gabor in 1948 [172], a lot of work was done by 
many researchers [173-175]. An intensity modulation hologram encodes the intensity 
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(b) Reconstruction 


Figure 6.2: a) The intensity interference pattern of a reference wave Ur and the object 
Uo is recorded on a holographic plate, b) The object wave is recovered after the interference 
pattern is shone onto the plate. Both hgures have been copied from [45] with the permission 
of Prof. M.C. Teich. 


and the phase as a form of intensity pattern^. In order to do that, an interference 
pattern needs to be recorded so that the information of the phase of the wave is carried 
through. Then, if the wave that wants to be reconstructed is named object wave 
Uo and the wave that will impinge on the hologram is called reference wave Ur, the 
transmittance needs to be proportional to their interference: 

tOc\Uo + Ur\^ = Ir + Io + U^Uo + UrU^ (6.5) 

where Uo, Ur, and consequently all the rest of terms can take different values for 
different transverse coordinates {x, y). When the reference wave illuminates the trans¬ 
parency, the following wave is obtained as a result: 

U = tUr(X Ur{Ir + lo) + UUo + U^U^ (6.6) 

Any wave can be chosen to illuminate the transparency. It is convenient to choose a 
wave such that U is a constant, i.e. Ir{x,y) = U- If that is the case, equation (6.6) 
can be re-formulated as: 

U{x, y) (X Ir + Io{x, y) + \/TrUo{x, y) -f \/TrU*o {x, y) (6.7) 

^The first holograms were made of photographic plates, which indeed are only sensitive to the 
intensity of light shone on them. 
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sin 

Figure 6.3: a) Recording of the hologram, b) Reconstruction of the hologram, giving rise 
to three components: a constant wave, Uoexp{ik sin 9x) and U* exp(—iA: sin 0x) Both figures 
have been copied from [45] with the permission of Prof. M.C. Teich. 

Now, to make sure that the four terms of the equations are properly separated in space, 
a oblique phase can be added to Uo- That is, if U^{x,y) = Uo{x,y) exp{ikx sin 9), then 
equation (6.7) can be expressed as: 

U'{x, y) oc Ir+Ioix, y) + \/TrUo{x, y) exY>{ikx sin 9)+ ^/TrUl (x, y) exp(—iA:xsin6*) (6.8) 

The hrst two terms of equation (6.8) are waves that travel in the axis. However, 
the third and fourth terms of the equations are waves travelling in an oblique manner, 
as shown in Figure 6.3. The third term represents the object wave Uo travelling at an 
oblique angle 9 with respect to the 2 ; axis and modulated by a constant factor ^/U. 
On the contrary, the fourth term represents the conjugate wave U* travelling at an 
angle —9 with respect to the z axis. Hence, the addition of an oblique plane wave 
allows for the total reconstruction and isolation of the object wave, which is the aim of 
any hologram. In fact, the addition of an oblique phase is widely used in holography, 
and will also be considered in sections 6.2 and 6.3 when phase-only holograms are 
introduced. 




6.2 Computer Generated Holograms 


The meaning of the word hologram changed with the advent of modern computers. 
Before then, a real object wave needed to be encoded in the hologram. The invention 
































































100 


Experimental Techniques 


of computers allowed for the creation of holograms representing optical waves coming 
from non-existing objects. Holography became an inverse scattering problem, where 
the object wave was given and the diffraction pattern was sought [176]. The hrst 
CGHs made by Brown and Lohmann [176] were a very specihc kind of amplitude- 
only modulated holograms: they were binary, i.e. a set of apertures on an opaque 
mask. More sophisticated methods to produce intensity-modulated CGH using gray- 
tones came up later, even though their fabrication remained challenging [177, 178]. 
A review of the history of intensity-only GGH can be found in [179]. Then, in the 
1970’s, another type of holograms was discovered [180]. Using phase-only holograms 
(which were already known), Kirk managed to encode amplitude information in a 
phase-only hologram. However, their use went almost unnoticed until 1999, when 
Davis and co-workers used SLMs (instead of photographic dims or transparencies) to 
encode amplitude and phase modulation in a phase-only hologram [181]. Now, even 
though phase-only GGH did not start being used to encode amplitude modulation until 
the 2000’s, they were very popular to change and encode phase shifts. The fabrication 
processes had been carefully characterized by Hariharan in the 1970’s [174, 182-187], 
and their use was extended in imaging [188], Fourier optics [189] and even acoustics, 
where they were hrst conceived [190, 191]. Here, I will describe an applications of phase- 
only GGH in the held of singular optics, and in particular in the creation of beams 
of light with a well-dehned AM. Light beams with phase singularities, also known as 
vortex or doughnut beams^, boomed after the seminal paper by Allen and co-workers 
in 1992 [110]. However, even when the connection between phase singularities and the 
AM of light had not been properly established, GGH had already been used to produce 
vortex beams [192-197]. In fact, they kept being mainly used to manipulate the AM 
of light up until the beginning of 2000, when the use of SLMs took over [198-201]. 
Next, I will describe why GGH are useful to manipulate the AM of light, and how they 
achieve it. 


^See section 4.2 to know a bit more about why they are called in this manner. 
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6.2.1 Phase singularities 

In chapter 2, Bessel beams and multipolar fields have been introduced. Both of them 
are EM modes, i.e. they are basis of solutions of the Maxwell equations (2.7). Also, 
they are both eigenstates of Jz- Therefore, following the nomenclature given in chapter 
4, they are cylindrically symmetric. Cylindrically symmetric beams have been used in 
chapters 4 and 5 to show many different effects - from excitation of WGMs, to inducing 
duality symmetry in non-dual materials. In chapters 4 and 5, a cylindrically symmet¬ 
ric beam has been decomposed as a superposition of multipolar fields, as these are the 
normal modes of a sphere. Unfortunately, multipolar fields are experimentally difficult 
to produce, as they propagate in all directions of space. 


Here, it will be considered that a cylindrically symmetric beam is a general super¬ 
position of Bessel beams (with either TE/TM character or well-dehned helicity): 


EL = 


kpdkp + 


h 


kpdkp ( ) (6.9) 


Now, because of their definition as a basis of solution of Maxwell equations, Bessel 
beams include paraxial and non-paraxial waves. However, as mentioned in section 6.1, 
only paraxial waves will be considered in this chapter. Therefore, equation (6.9) needs 
to be re-expressed in the paraxial approximation, i.e. kp -C k. Using equations (2.38, 
2.39, 2.42, 2.76, 2.77), it can be observed that the four different Bessel modes can be 
expressed in the paraxial limit as: 


B 


TE 
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^mkp 


B pm^^ik.z 


pJe{kpp) +Z0J5(/Cpp) 

pJs{kpp) +i^Je{kpp) 


B 


mkn 
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DkpJm-\-l(^kpp^C 


i(m+l)(l) ikzZ - 




( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 


with Akp,Bk^,Ckp,Dkp being functions of kp, and Jsikpp) = Jm-i{kpp) - Jm+i{kpp), 
■k{kpp) = Jm-i{kpp) + Jm+i{kpp). Taking into account the asymptotic properties of 
Bessel functions [69], which yield Jm{x) oc x™' when x —)■ 0, it can be seen that a 
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cylindrically symmetric beam in the paraxial approximation can be done with waves 
whose complex amplitude, when p —)■ 0, is of the kind: 

f/i(p,0,^) ocp'e*^V^^^ (6.14) 

and whose polarization has to be properly matched depending on the value of Jz that 
needs to be created. That is, if a value of Jz = m* is needed, then Ui{p,((),z) can be 
used if the beam is radially or azimuthally polarized, whereas hl(z-i) or 17(z+i) ought to 
be used for left or right circular polarization respectively. 

Polarization states such as (T+ and (T_ can be easily achieved using wave plates [45]. On 
the contrary, 0 and p polarization states are more challenging to obtain, as typically 
interferometry was needed [202, 203]. Nowadays, with the invention of liquid-crystal 
polarization converters [204, 205], radial and azimuthal polarization states are not as 
challenging to prepare as before. For simplicity though, only circular polarization states 
will be considered in all the experiments of this thesis concerning AM states. 

Independently of the polarization state though, it is clear looking at equations (6.10- 
6.13) that the beam must possess a phase dependence of the kind to be cylindrically 
symmetric beams. Figure 6.4 depicts the phase of Ui{x,y) on a transverse plane for 
I = —3,—1,1,4. It can be observed that the phase winds around the center of the 
transverse plane {x,y) = (0,0). In fact, {x,y) = (0,0) is a singular point, as the an¬ 
gle 0 is not well-dehned. Any point whose phase is not well-dehned is named phase 
singularity. A point which is singular in phase must have a zero intensity, so that 
the held still fulhls the Helmholtz equation. Then, the azimuthal number I gives the 
order of the singularity. Thus, the aim of the phase-only CGHs made for this project 
will be to create beams with a twisting phase such as In fact, as mentioned in 
section 6.1.1, an oblique phase of the kind exp^ik^x) will be added to the spiral phase 
gQ hologram separates the different diffraction orders coming out of the 

hologram. 
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Figure 6.4; Phase profile of a paraxial vortex beam on a transverse plane x — y: 
aig{Ui{x,y)} for a) I = —3, b) I = —1, c) Z = 1, d) f = 4. Black means y)mod27r] = 0, 
and white y)mod27r] = 27r 

6.2.2 Fabrication of phase-only CGH 

It is clear now that a hologram which imprints a phase dependence onto the beam 
is needed. Now, if a paraxial wave of the kind f/(r) = y4(r) exp(ffc^; 2 ) goes through 
a hologram with t{x, y) = the held after the hologram will be a vortex beam. 
Equation (6.4) gives a way to produce that kind of hologram. If a homogeneous and 
isotropic material is modihed in a way that its thickness has a 0 dependence, then 
t{x,y) can be made equal to e®^'^: 

t{x,y) ~ hoexp[—i{n — l)kd{x,y)] 
t{x,y) = exp(f/0) 

That is, a homogeneous material made of an index of refraction n can be used to create 
phase singularities of order I if its thickness varies according to the following expression: 


d{x,y) = d(0) 


(6.15) 


d(0) = 


l(j)X 


{n — 1)271 


(6.16) 
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Note that equation (6.16) is wavelength-dependent, what means that in theory the 
thickness d{(j)) to create a perfect phase dependence can only be achieved for a single 
wavelength. However, in practice, the efficiency of CGHs does not considerably vary 
for a broad range of wavelengths. The method to obtain d{(j)) is based on the following 
idea. As it was explained in section 6.1.1, photographic dims are only sensitive to 
the intensity of light. In a conventional black and white negative developing process, 
a him transforms the intensity captured by the camera into a grey scale where black 
means high exposure and white low exposure [206] . Then, the negative (which could be 
understood as an amplitude hologram) can be transformed into a phase hologram by a 
bleaching and hxing process. The bleaching and hxing processes remove all the silver 
salts and the unexposed gelatin, consequently giving different widths to the material 
depending on the exposure. More importantly, the efficiency of the hologram increases 
[207]. The fabrication process is detailed in Appendix B.l. 

6.2.3 Characterization 

The performance of CGH can be characterized by different parameters. One of the 
most relevant ones is their efficiency. Here, I dehne what the efficiency of a CGH is, 
and how to measure it. Then, in Appendix B.2, the characterization of the CGH fab¬ 
ricated throughout my project and their dependence on some controllable parameters 
is detailed. The wavelength of the laser used to perform the characterisation of the 
holograms is A = 632.8nm. 

Even though the efficiency is not the only parameter to characterize the performance 
of a hologram, it is certainly a very important one. The efficiency is a ratio between 
the power of the incident beam onto the hologram, and the power going to the hrst 
diffraction order of the hologram (that is, the object wave in Figure 6.3). Note that 
if the hologram absorbs some power, then the efficiency drops. As depicted by Figure 
6.5, the intensity is measured with a CCD camera. That is, a snapshot is taken of the 
incident beam when the hologram is removed from the optical path (Figure 6.5(a)). 
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Figure 6.5: a) CCD image of the reference beam when no hologram is on the optical 
path, b) Typical diffraction pattern created by a CGH when the efficiency is over 30%. The 
0th order is at the same position as the reference beam, whereas the 1st order has been tilted. 
The colour bar indicates goes from 0 to 1, where 1 is taken as the absolute intensity maximum 
of the reference beam. 


Next, the hologram is added in the optical path making sure that the beam almost 
hlls the hologram dimensions. The hologram tilts the incident beam and its effect is 
again captured by the CCD camera (Figure 6.5(b)). Due to the imperfections of the 
hologram, other diffraction orders arise giving rise to the pattern in Figure 6.5(b). Note 
that in Figure 6.5(b), there is light at the same position where the incident beam is 
in Figure 6.5(a). This is the so-called Oth-order, i.e. the incident beam going through 
the hologram without being modified by it. The tilt given by the hologram depends 
on the number of lines of the hologram, i.e. the number of complete 27r phase jumps. 
In fact, the angle of diffraction given by the hologram is given by: 




A 

Ax 


( 6 . 17 ) 


where ^dx is the angle that the 1st diffraction order forms with the direction of the 
initial beam (or the 0th order); Ax is the width of a 27r phase ramp in the hologram 
(see Figure B.l in Appendix B.l); and x is the direction in which the phase jump is 
given. In Appendix B.4, it will be determined how the efficiency critically depends on 
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the oblique phase exp(i/ca;x) added to the hologram, which is proportional to (Ax)“h 
Also, a detailed explanation of how to experimentally separate the different diffraction 
orders coming out of a hologram, i.e. how to obtain Figures 6.5, 6.8 and 6.9, is given 
in Appendix B.3. 

6.3 Spatial Light Modulators 

The technology to digitally control nematic liquid crystals boomed at the end of the 
1980’s giving rise to SLMs that could modify the phase, the amplitude and the polar¬ 
ization of beams [208-211]. A sketch of a transverse section of an SLM can be observed 
in Figure 6.6. Two electrodes close a volume that is filled with liquid crystals and a 
mirror. First, the incident light on the SLM is transmitted through the coverglass 
electrode. Afterwards, it goes through the volume enclosing the liquid crystals. Their 
orientation is determined by the difference of potential between the two electrodes, 
which is at the same time linked to the CGH displayed on the computer. Due to 
the anisotropy of the liquid crystals, their index of refraction will be different for two 
orthogonal incident linear polarizations. Then, controlled changes in their orientation 
will vary the index of refraction for one of polarization as a function of the transverse 
coordinates {x,y), giving rise to a computer-controlled transmittance t{x,y)‘^. Then, 
the light hits a mirror and is reflected back off. The mirror sits on top of a grid of 
pixel electrodes, which set the potential difference on the liquid crystal volume in con¬ 
junction with the hrst electrode. Thus, in most of cases, the SLM technology works in 
reflection, in contrast to the CGH holograms described in section 6.2.3, which work in 
transmission (see Figure 6.7). One decade after their boom, the ability of SLMs to in¬ 
teractively control holograms in real time made them very prolific in the field of optical 
micro-manipulation and optical tweezers [199, 213-217]. Nowadays, the holographic 
optical tweezers are used in a myriad of helds and it is even possible to hnd portable 
commercial optical tweezers set-ups [218, 219]. Now, as commented in section 6.2, the 

'^Unlike the CGHs, because of the anisotropy of the liquid crystals, t{x, y) is very polarisation- 
dependent. 
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Figure 6.6: Schematics of a transverse section of an SLM. Two electrodes and a mirror 
enclose a volume filled with liquid crystals. The potential difference applied by the electrodes 
sets the orientation of the liquid crystal, thus achieving the desired transmittance function 
t{x,y). Figure copied from [212] with their permission. 


use of SLMs to produce beams carrying phase singularities highly escalated after 2002 
[198-201]. Actually, the use of phase-only CGH for these purposes is practically obso¬ 
lete now. Nonetheless, in this section I will characterize the performance of SLMs and 
compare it with the performance of CGH. It will be seen that even though CGHs are 
less efficient and certainly less flexible, they can still be useful for certain applications. 

6.3.1 SLM characterization 

Here, similarly to what has been done in section 6.2.3, the performance of an SLM 
is characterized. The SLM used in this chapter as well as in chapter 7 and 8 is a 
Boulder Nonlinear Systems^ (BNS) XYseries optimized for wavelengths ranging from 
[630 — 800] nm. The SLM characterization is done again in terms of the SLM efficiency. 
Figure 6.7(a) shows a real picture of the set-up to test the efficiency of the SLM. The 
laser beam is given by a He-Ne laser, whose A = 632.8nm. The beam is expanded 
to properly fill the SLM chip (see Figure 6.7(b)). Then, its polarization is controlled 


® http : / / WWW. bnonlinear. com / 
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Figure 6.7: a) Picture of the experimental set-up to characterize the XY series SLM from 
BNS. A Gaussian beam from a He-Ne laser is expanded, then its polarization is controlled 
before it is sent over to the SLM. The SLM changes the phase and reflects off the beam, 
which is focused with a lens onto a CCD camera, b) Picture of the SLM. The incident beam 
almost hlls in all the SLM chip. 


by a linear polariser and a half-wave plate. Later, a mirror sends the beam to the 
SLM, which changes its phase. Finally, a lens Fourier-transforms the beam and a CCD 
camera records the image at the Fourier plane. Then, the efficiency is measured as the 
intensity of the laser beam hitting the SLM divided over the intensity going to the hrst 
diffraction order. Note that this efficiency measurement takes into account the losses 
of the SLM: imperfect reflection of the mirror, as well as diffraction losses due to the 
grid of pixels. Similarly to Figure 6.5, a typical diffraction pattern is captured with a 
CCD camera in Figure 6.8. Figure 6.8(a) shows the position of the beam on the CCD 
camera when the SLM works as a mirror (the phase shift is 0). Then, Figure 6.8(b) 
displays the diffraction pattern created by the SLM. Note that unlike the diffraction 
pattern created by the CGH shown in Figure 6.5, most of the intensity goes to the 
1st diffraction order. In Appendix B.4, the efficiency of the SLM is characterized as a 
function of different controllable parameters in the laboratory. 
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Figure 6.8: a) CCD image of the reference beam when the SLM works as a mirror, b) 
Typical diffraction pattern created by the SLM when the efficiency is over 80%. The 0th 
order is at the same position as the reference beam, whereas the 1st order has been tilted. 
The colour bar indicates goes from 0 to 1, where 1 is taken as the absolute intensity maximum 
of the reference beam. 

6.3.2 CGHs or SLMs? Pros and cons 

Due to their flexibility and interactive control, SLMs have taken over the role that 
CGHs used to play in the held of optical micro-manipulation, and in particular in the 
held of the AM of light. However, CGHs made of photographic hlms can still be useful 
for some applications. Thus, to conclude with this study of diherent methods to create 
phase-only holograms, the pros and cons of using each of the two methods will be 
discussed. 

• Price . Making CGH in the way described in section 6.2 is cheap. All the 
equipment and materials needed are easily available. On average, making 60 
holograms (which corresponds to two photographic hlms) costs less than $100. 
In contrast, the cheapest SLM in the market (from Cambridge Correlators®) is 
valued $1000. And then, the rest of them are in the range of $15000-$25000. 

• Efficiency. The CGH made and characterized in this thesis reached efficiencies 



® WWW. cambridgecorrelators. com / 
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Figure 6.9; CCD image of a typical diffraction pattern created by a CGH and a SLM. It 
is seen that the SLM sends most of the light to the first diffraction order, whereas the light 
coming out of the CGH is split between many diffraction orders. 

values of 45%. However, it is known that the efficiency can reach values of the 
order of 60%. The main reason why those values were not reached is because the 
chemical process is not optimised and some silver has remained in the sample, 
thus producing some absorption. In contrast, the efficiencies measured with the 
SLM from BNS are constantly above 80%. In order to visualize the difference 
of efficiency. Figure 6.9 displays the diffraction pattern from a CGH and the one 
from an SLM. It is seen that with the SLM, almost all the light goes into the 1st 
diffraction order, whilst with the CGH the 0th and 2nd order are also clearly seen. 
The cheap SLMs from Cambridge Correlators do not have such high efficiencies, 
though. Theirs is usually below 60%. 

• Polarization . Polarization does not play any role for CGHs. That is, the 
transmittance t{x, y) given by the hologram is completely independent of the 
polarization of the incoming beam. However, for SLMs, the polarization of the 
incident beam is crucial. Nematic liquid crystals form an anisotropic medium, 
therefore a change of polarization alters their response dramatically. In order to 
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Table 6.1: Summary features CGHs and SLMs. 



CGH 

SLM 

Price 

$100 to make 60 holograms. 

Cambridge Correlators: $1500. 
Others: $15000-$25000. 

Efficiency 

40% easily achievable. Chemistry 
needs to be rehned to go above 
50%. 

Cambridge Correlators: 55%. 

Others: > 60%, 85% with BNS. 

Polarization 

Does not play any role. 

Efficiency highly depends on polar¬ 
ization. 

Flexibility 

Cannot be changed once it is made. 

Controlled by computer. Can be 
changed interactively. 

A dependence 

Phase shift depends on A. Effi¬ 
ciency is not very affected by A. 

Performance depends a lot on A. 

Set-up 

Transmission. 

Reflection. 


achieve the desired phase shift, the incident beam must be linearly polarized in 
a certain direction given by the axis of the liquid crystals [208]. 

• Flexibility . CGH are optical elements that cannot be varied once they have 
been made. SLMs are dynamically controlled, therefore they can steer the beam 
interactively in real time. This makes them much more flexible. 

• A DEPENDENCE . The phase shift given by both CGHs and SLMs is wavelength 
dependent (see equation (6.4)). However, for the CGH, the efficiency remains 
almost constant for a large range of wavelengths. In contrast, the efficiency of an 
SLM is very dependent on the range of As for which it has been designed. When 
the SLM is tested outside that range, the efficiency drop considerably. One of 
the main reasons for the drop in efficiency is the drop in reflectivity by the mirror 
standing between the liquid crystals and the pixel electrodes. 

• Set-up . CGH work in transmission, whereas SLMs typically work in reflection. 


A summary of the pros and cons can be found in Table 6.1. 
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Figure 6.10: SEM images of two nano-apertures of diameter = 237nm (a) and b)) and 
= 424nm (c) and d)). None of the images is post-processed. Images b) and d) have been 
used to retrieve the values of 4>, Eccentricity and Ellipticity 


6.4 Fabrication of samples 

In order to carry out the experiments done in chapters 7 and 8, different samples 
have been needed. In this section, both the samples and the methods to obtain them 
are broadly described. Notice that this is not the main topic and this thesis and 
consequently the level of detail of this section is lower than some others. 

6.4.1 Nano-apertures on metallic films 

In chapter 7, some experiments using circular nano-apertures are presented. I assisted 
Alexander E. Minovich in the fabrication of these circular nano-apertures at the ACT 
Node of the Australian Nanofabrication Facility, in Canberra. First, a glass him of 
1mm thickness was cleaned using a plasma etching technique. Afterwards, a gold 
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layer of 200nm was deposited on top of the glass layer using thermal evaporation. 
Then, the gold layer on the glass him was milled using a FEI Helios NanoLab Focused 
Ion Beam system^. The circular nano-apertures have diameter sizes ranging from 
200-450 nm, and they are all 50 /im apart, so that the coupling of surface plasmons 
launched from one nanohole to the closest one is avoided. The nano-apertures have been 
imaged using a secondary electron Scanning Electron Microscope (SEM, JEOL JSM- 
6480) operated at lOKeV, at the Microscopy Unit of Faculty of Science in Macquarie 
University. Then, the images have been analyzed with Matlab, where the boundaries 
of the nano-apertures were determined by selecting the pixels whose intensity was 
below the 10% of the maximum. Figure 6.10(a) depicts a SEM image of a 237 nm 
nano-aperture and Figure 6.10(b) shows the same image with a boundary traced by a 
Matlab code to determine the diameter (4)), the eccentricity and the ellipticity of the 
hole. Likewise, Figure 6.10(c) and Figure 6.10(d) show the same features for a larger 
aperture of 4) = 424 nm. All the images on Figure 6.10 have not been post-processed. 
The results for the 12 milled nano-apertures are summarized in Table 6.2. 

6.4.2 Single spheres on a substrate 

The theory developed in chapters 4 and 5 is based on the assumption that an isolated 
sphere is embedded in an isotropic, homogeneous medium. It is not trivial to achieve 
these conditions experimentally. One option would be levitating a single sphere. This 
was done already in the pioneering experiments of Ashkin [220-223]. However, levi¬ 
tating particles of diameters smaller than 1 pm is experimentally hard. Even though 
objects of the order of 100 nm have been recently levitated [224, 225], the experimental 
set-ups typically use vacuum chambers for cooling and other equipment that was not 
available in my current set-up. Now, due to the fact that we wanted to be able to 
use small particles to reach the dipolar regime to reproduce results from chapter 5, 
another technique to have them embedded in a homogeneous medium was considered. 
This technique is spin-coating. 

^http://anff-act. anu.edu.au/HTM/detail/flagship_services.htm 
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Table 6.2; Characterization of the circular nano-apertures images taken with the SEM. 


$ (nm) 

Eccentricity 

Ellipticity 

212 

0.324 

0.0538 

214 

0.414 

0.0897 

237 

0.422 

0.0933 

253 

0.462 

0.113 

317 

0.295 

0.0446 

325 

0.279 

0.0396 

333 

0.118 

0.00693 

341 

0.360 

0.0672 

424 

0.260 

0.0346 

429 

0.388 

0.0783 

432 

0.273 

0.0381 

433 

0.310 

0.493 


Spin-coating is a simple technique that is used to uniformly deposit a layer of ma¬ 
terial (usually a polymer or an organic material) on top of a flat surface. A good 
review of this process can be found in [226]. The process starts by mixing particles 
with a resin. A droplet of the mix (with a controlled volume) is deposited on top 
of a planar surface, a microscope slide in this case. The microscope slide has been 
previously plasma-etched to remove any impurities that the glass surface could have. 
The microscope slide with a droplet on top is located on the spin-coater. Three pa¬ 
rameters need to be entered to control the spin-coater: the duration time (tsc), the 
maximum angular speed (cj^c), and the angular acceleration (osc) to reach Ugc- Then, 
the spin-coater spins for tgc seconds following a pattern assigned by age and Ugc- Once 
the process is hnished, the sample is removed and placed under the illumination of a 
UV lamp, which joins the monomers of the resin to form a solid polymer. In this way, 
after the last step, the particles hnd themselves trapped in a homogeneous medium. 
If the concentration of particles is right (neither too low nor too high), it will be easy 
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to find particles under the microscope and some of them will be isolated, i.e. single 
particles whose closest particles are at a distance of 10 As. 

As it has been explained in chapters 4 and 5, the relative index of refraction between 
the single spheres and the surrounding medium is a crucial parameter to observe many 
phenomena. In fact, the example given at section 4.4 and Figure 5.2 clearly show 
that high play an important role to either excite WGMs, or to preserve helicity. 
The refractive index of polymers varies from 1.3-1.7 [227]. The most common ones for 
lithographic processes, such as PMMA or SU-8, have n ~ 1.5 [228, 229]. That means 
that with these typical polymers, particles with n > 2.4 are needed in order to obtain 

Ur > 1 . 6 . 

Finding spheres with high refractive indexes is hard. Silicon is a very good candi¬ 
date working at telecommunication wavelengths with n > 3, but it absorbs too much 
in the optical region [168]. In the experiments carried out in chapter 8, particles made 
of Titania (Ti02) have been used. Titania can be found in different phases, amorphous, 
anatase and rutile. The three phases have different refractive indexes. At A = 500 nm, 
rutile has rim = 2.7 [230], anatase has rian = 2.3 [231], and the amorphous phase 
Ram = 1-8 [231]. Spheres made of these materials have been purchased from a Cana¬ 
dian company, MKnano®. 

In order to increase even further the contrast, a low-refractive index resin from an 
Israeli company called AlY Polymers Ldt.® was acquired. The refractive index is 
Tire = 1.33. However, mixing titania particles with the resins acquired from MY Poly¬ 
mers has turned out to be challenging. The resins are highly hydrophobic, that is, 
they repel water. Whilst titania particles are given in aqueous solutions as they are 
hydrophilic. The way to fix this problem is to mix the aqueous solution in methanol. 


® http: / / WWW. mknano .com/ 
® WWW. mypolymers .com/ 
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Figure 6.11: Dark-field image of ‘k = 950 nm titania particles deposited on a coverslip 
after being spun. The images have been obtained from directly above the coverslip. 

which is also hydrophobic. Xavier Vidal and I have seen that when the aqneous so- 
Intion of titania particles is mixed with methanol in a 1/10 or 1/100 water/methanol 
manner, then the mix of titania and resin give acceptable resnlts. That is, even thongh 
lots of particles still agglomerate, some of them can be isolated and snrronnded by 
an homogeneons and isotropic medinm. However, we have not been able to prodnce 
samples where the air-polymer interface is parallel to the coverslip snrface. Since we 
are interested in symmetric samples, we have dismissed them dne to these asymmetries. 

Instead, we have done something else. We have deposited the aqneons solntion of 
titania particles directly on top of a coverslip. The sample has been spnn nsing the 
spin-coater to redistribnte the particles across the sample and get some isolated single 
spheres. The resnlt is displayed by Fignre 6.11, which is a dark-held image of a certain 
region of one of onr samples. Even thongh these samples do not have spherical sym¬ 
metry, they still keep the cylindrical one. Thns, they are still good candidates to test 
the theories developed in chapters 4 and 5, as they must preserve J^. 



“For a succesful technology, reality must take precedence over public 
relations, for nature cannot be fooled." 

Richard P. Feynman 
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Experiments with single nano-apertures 


7.1 Introduction 

In this chapter, a way of controlling the scattering of a given nano-structure with the 
incident field is shown. As it was explained in section 2.6, the conventional way to 
control the scattering response of a structure is by characterizing its scattering matrix 
S(g,m), where g and m are general sets of variables describing the geometrical and 
material properties of the structure, respectively. Then, assuming a ^ or p plane wave- 
excitation, the geometrical and/or material properties of the structure are modified to 
control the optical response given the corresponding plane wave excitation. There is 
an alternative and less-explored way of controlling the response of structures, and that 
is by tuning the properties of the incident beam. That is, once the scattering matrix 
has been characterized, a certain response can be obtained as a weighted superposition 
of individual responses of different plane waves: 

E°"*(r) = (^j d^k'i exp(ik'i ■ r)S(g, m) * E”(k'i)^ (r) (7.1) 

where kl is each of the different plane waves that take part in the superposition, and 
E“(k'i) is the Fourier transform of the incident field, which modulates the plane wave 
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decomposition. Both approaches are equally challenging and complementary. For ex¬ 
ample, once a structure is characterized, the knowledge of the response to different 
input beams can be used to control its scattering. Alternatively, once the structure is 
chosen but not characterized, the incoming beams used to probe it can be modified 
to maximize the information retrieved from the sample. In this sense, exploiting the 
symmetries of the sample allows for identifying correlations between input and output 
fields. 

Here, the scattering of a non-chiral sample is controlled with vortex beams to ex¬ 
perimentally induce a huge CD. These results are presented in section 7.3, but before 
that an analysis of the experiment in terms of symmetries is carried out. 


7.2 Symmetries 

The interaction of light with nano-apertures is a problem that has been carefully 
studied by many scientists. In particular, the work of Ebbesen et al. [232] show¬ 
ing that nano-apertures could have an extraordinary transmission due to the coupling 
of light and SPPs opened up lots of possibilities among different fields in optics, micro¬ 
manipulation, biophysics and condensed-matter [233-236]. Lots of these studies have 
focused on how the nano-apertures couple with SPPs depending on many different 
parameters [237-240]. Some others have focused on the radiation diagram of these 
nano-apertures and have studied how the s' and p polarization states are transmitted 
through the apertures both theoretically and numerically [91, 241-245]. Finally, the 
excitation of surface plasmon with subwavelength nano-apertures has also been used 
to induce spin-orbit interactions [35-37, 246-248]. 

Here, it will be shown that a symmetry characterization of both samples and exci¬ 
tation beams can allow for a predictive description of light-matter interactions at the 
nano-scale. First, the symmetries of the nano-apertures shown in section 6.4.1 are 
characterized. As described in section 6.4.1, the sample is a 1mm glass film coated 



7.2 Symmetries 


119 


with a gold layer of 200nm. The sizes of the nano-apertures as well as their devia¬ 
tion from a circle are given by Table 6.2. Then, with a very good approximation, the 
circular nano-apertures are symmetric under rotations along the axis, where the 2 ; 
axis stands for the normal direction of the sample. With the same assumption, they 
are also symmetric under any mirror transformation that contains the z axis. Last 
but not least, since the samples behave as linear media (see section 3.2), they are also 
symmetric under time translations and therefore the frequency of light is preserved in 
the interaction. In contrast, spatial translations and duality symmetry are broken. All 
these symmetry considerations of the samples are inherited by the scattering matrix 
S(g,m). Thus, the scattering matrix commutes with the following three operators^: 

[S,R,] = = [S,TAt] = 0 (7.2) 

Now, even though the experiment will be done with a single nanohole, the CD study 
will be carried out with a more general system formed by two MOs and the sample 
(see Figure 7.1). This more general system will be referred to as target (T). The 
idea is that the microscope objectives are symmetric under the same transformations 
as the sample: rotations around the z axis, mirror symmetries containing the z axis 
and temporal translations (see section 2.6). Thus, target and sample are equivalent 
in terms of symmetries. This fact will be exploited in the study of CD in section 7.3. 
Moreover, it will be considered that the action of the two microscope objectives and 
the sample can be characterized by a linear integro-differential operator that inherits 
all the properties of the target. I will denote this operator as T. Similarly to the 
scattering matrix, T fulfils the following commutation rules: 

[T,R,] = = [T,TAt] = 0 (7.3) 


7.2.1 Duality symmetry 

One main difference between the MOs and the sample is that the MOs preserve helicity 
(see section 2.6), whereas the sample does not. In this subsection, the mechanism by 


^See section 2.2 to see the definition of these operators 
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Figure 7.1: Schematics of the Target = MOs + Sample. MOi focus light onto the sample 
and the scattering is collected by MO 2 . 

which duality symmetry is highly broken in the sample is explained. As it has been 
explained in previous chapters, the helicity operator A is the generator of duality trans¬ 
formations. Therefore, the helicity of light is not preserved upon interaction when the 
sample is not dual, and vice versa. 

It is known that the microscopic Maxwell equations for material media always break 
duality symmetry [41, 42]. However, it was proven in [41] that duality symmetry can 
be restored in the macroscopic Maxwell equations if there is an impedance matching 
between all the media where the studied EM interaction takes place. That is. 


— = const A conservation 


( 7 . 4 ) 


for any medium made of an arbitrary number of isotropic and homogeneous sub-media 
i, Cj and fii being its electric permittivity and magnetic permeability. This is not the 
case for the sample, made of gold and glass. Clearly, 


^glass / ^air / ^gold 


( 7 . 5 ) 



therefore duality symmetry must be broken by the sample. 


In contrast, as shown in section 2.6, MOs are made in a way such that helicity is 
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Figure 7.2: Schematics of a helicity projection scheme. A beam with a well-defined 
helicity is focused down to a sample that breaks duality symmetry. As a consequence, an 
additional helicity component is created. A quarter-wave plate and a linear polariser are used 
to separate the two helicity states. 



preserved^. That is, consider a paraxial beam of light with a well-dehned helicity p. 
Because the beam is paraxial, the polarization has to be circular (dp). Now, it is 
well-known that when a beam is tightly focused, the polarization changes so that the 
transversality of the held is maintained [38, 49, 54]. However, because MOs preserve 
helicity, the helicity of the beam at the focus is still p. 

Then, the interaction with the sample takes place. The sample is not dual-symmetric, 
therefore helicity is not preserved. That is, an opposite helicity component —p is cre¬ 
ated upon scattering. 

The scattered light is collected by the second MO, which also preserves helicity. If 
the collecting MO is placed at a focal distance of the sample, then the beam coming 
out of the objective is a paraxial held again. Hence, one component of the helicity is 
left circularly polarized (TCP) and the other one is right circularly polarized (RCP). 
TCP corresponds to d+ and RCP to d_. Then, by projecting the states back to a linear 
basis, the two diherent helicity components can be analysed separately. 

The method is depicted in Figure 7.2. Between MOi and QWP, the red helix rep- 

^As a consequence, the action of two MOs preserves the polarisation state. This statement has 
been experimentally verified in the laboratory using a Thorlabs polarimeter. The value of the third 
Stokes parameter S 3 remained the same when it was measured with or without the MOs. 
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resents a beam with a helicity p, and the blue helix represents a beam with a helicity 
—p. It can be seen that the beam has helicity p before it interacts with the sample. 
Once the interaction has been taken place, the beam is a combination of red and blue, 
as a —p helicity component has been created. The quarter-wave plate (QWP) changes 
the polarization state from LCP and RCP to horizontal and vertical, so that the linear 
polariser (P) can filter out one of the components. Note that the —p helicity compo¬ 
nent has an exceptional signal ratio. That is, all the light analysed with the —p helicity 
component exclusively comes from the sample and not from the incident held, except 
for the experimental errors introduced by the QWP and LP. 

The fact that a nano-aperture scatters light in the crossed helicity component is not 
surprising. Any macroscopic material which does not fulhl equation (7.4) breaks dual¬ 
ity symmetry and creates light with the crossed helicity component —p. As I showed 
in chapter 5, some very specihc conditions need to be fulhlled so that dielectric spheres 
behave in a dual way. However, the creation of modes with helicity —p with a circular 
nano-aperture is interesting from different perspectives. 

Firstly, due to the sub-wavelength nature of the nano-apertures, duality symmetry 
is highly broken. To analyse it, it is useful to dehne a ratio of light scattered into the 
crossed helicity component divided over the direct one, 7 . The dehnition of 7 is very 
similar to the transfer function Tp{R,nr) defined in section 5.3, except for a subtlety: 
Tp{R,nr) takes into account the intensity which is backscattered, whereas 7 is only 
dehned in transmission. If the 7 ratio is computed for the glass -|- gold multilayer sys¬ 
tem, it turns out to be of the order of 10“^ [169]. However, the nano-aperture highly 
breaks duality symmetry and then 7 is of the order of 1 [169]. This incredibly high 
helicity conversion ratio is due to the coupling of the scattered field to surface modes 
of the multilayer system, whose character is TE/TM [169, 243]. Thus, the study of the 
—p helicity component as well as its intensity ratio with respect to the direct one can 
allow for the characterization of the coupling of light to surface modes. 
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Secondly, the creation of modes with helicity —p can also be very useful to link ge¬ 
ometrical and material symmetries of the sample to spin-orbit interactions. In the 
literature, spin-orbit couplings are used to explain phenomena where the polarization 
(spin) and spatial properties of light (orbital) are linked. Due to the non-paraxial 
nature of the helicity formalism, spin-orbit interactions are naturally included. Thus, 
in section 7.4, it will be shown that the two helicity components give rise to different 
spatial modes with different spatial properties and polarizations. Last but not least, 
the spatial properties of the crossed helicity component have been also used in [249] to 
develop a highly stable nanopositioning method to center a sample in only three steps. 

7.3 Circular Dichroism 

In this section, the second method to control the scattering of a structure explained in 
section 7.1 is used to induce giant CD on a sample. Since its discovery in the 19th cen¬ 
tury, CD has been widely used in science. Defined as the differential absorption LCP 
and RCP [250], its uses are as diverse as protein spectroscopy, DNA studies and char¬ 
acterization of the electronic structure of samples [251]. In the advent of nano-photonic 
circuitry, a lot of work has been put recently into characterizing plasmonic components 
and optical metamaterials in terms of CD [252-260]. The approach in all these works 
is the same one. CD wants to be induced with a structure that is illuminated under 
a plane wave excitation. In order to achieve it, structures with a non-trivial geometry 
made of a smart combination of materials are engineered. In fact, in all the cases cited 
before, the structures engineered are chiral. A chiral structure is such that none of 
its mirror symmetric images can be superimposed the sample itself [250, 261]. One of 
the mains reasons why the design of the structures is chiral is an important theorem 
known in molecular physics: 

“We know today, in molecular terms, that the one necessary and sufficient 
condition for a substance to exhibit optical activity is that its molecular structure 
be such that it cannot be superimposed on its image obtained by reflection in 
a mirror.” 

David M. Bishop, Group theory and chemistry, (2012) page 20 
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The approach that it is taken here to achieve CD is completely different. The nano¬ 
structure that will be used to induce CD is a circular nano-apertures. Typically it was 
thought that a non-chiral sample could not induce CD. However, recent experiments 
with non-chiral plasmonic planar structures have shown that a non-zero CD can be 
obtained with a sample that lacks mirror symmetry with respect to the incident input 
beam, but it is non-chiral on its own [262-265]. Other attempts to create CD with a 
non-chiral sample have been carried out surrounding the sample with a chiral medium 
[266, 267]. Here, I will go a step further and induce giant CD in a non-chiral sample 
such as a circular nano-aperture using light at normal incidence. In contrast to the 
previous approaches, the mirror symmetry of the system is broken with an internal 
degree of freedom of the input beam: its angular momentum. The results presented 
in this section show that CD can be induced in a symmetric sample if the beam used 
to make the CD measurement is a vortex beam. The reasons behind this interesting 
phenomenon are two. Firstly, the two LCP and RCP beams used to carry out the 
CD measurement are not mirror symmetric. And secondly, the LCP and RCP vortex 
beams have an AM that differs in two units. 

A sketch of the experimental set-up used in the experiment is depicted in Figure 7.3. 
It can be divided into three parts: preparation of states, non-paraxial interaction with 
the sample, and measurement. For the state preparation, a HeNe laser working at 
A = 632.8nm was used. It produces a collimated, linearly polarized Gaussian beam. 
The beam is expanded with a beam expander (lenses L 1 -L 2 ) to match the dimensions 
of the chip of the SLM. Both lenses are achromatic with an anti-reflection coating for 
650 — 1050nm. Their focal lenses are fi = 200nm and /2 = 300nm respectively, making 
the magnitude of the beam expander equal to 1.5. Then, the polarization state of the 
beam is modified with a linear polariser (Pi) and a half-wave plate (HWP) to maximize 
the efficiency of the SLM. The linear polariser is a Glan-Taylor Galcite Polarizer with 
a 650 — 1050nm coating. The HWP is an achromatic zero-order wave plate from New¬ 
port (10RP02-28). Then, the SLM creates a vortex beam by displaying an optimized 
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c) 

d) 


1 Mm 


Figure 7.3: (a) Schematic view of the optical set-up in consideration. LCP and RCP 
vortex beams go through a sub-wavelength circular aperture and their transmissivities are 
recorded with a CCD camera, (b) A pitchfork-like hologram is used to create a vortex beam, 
(c) Intensity profile of one of the vortex beams used in the experiment, (d) SEM images of 
two isolated circular nano-apertures. 


pitchfork hologram [268] (see Figure 7.3(b)). Proper control of the pitchfork hologram 
allows for the creation of a phase singularity of order I in the center of the beam, i.e. 
the phase of the beam twists around its center from 0 to 27il radians in one revolution. 
Note that when I = 0, the SLM behaves simply as a mirror. As explained in section 
6.2.3 and appendix B, the SLM produces different diffraction orders. In order to filter 
the non-desired orders of diffraction, a modified 4-f filtering system with the lenses L 3 - 
L 4 and an iris (I) in the middle is used. Lens L 3 (achromatic AC254-150-B-ML from 
Thorlabs) Fourier-transforms the beam, and the iris selects the 1st diffraction order 
and filters out the rest. Then, lens L 4 (achromatic AC254-125-B-ML from Thorlabs) 
is used to match the size of the back-aperture of the microscope objective that will be 
used to focus the beam down to the sample. But before focusing, the preparation of 
the input beam is finished by setting its polarization to either LCP or RCP. This is 
done with a linear polariser (P 2 ), which is also a GTIO-B like Pi, and a quarter-wave 
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plate (QWPi). The QWPi is a 10RP04-28 zero-order quartz wave-plate from Newport. 
Since the beam is collimated, this change of polarization does not appreciably affect 
the spatial shape of the input beam. 

After this initial preparation, the light is focused down to interact with the sample 
(S), i.e. one of the single circular nano-apertures listed in Table 6.2. The focus¬ 
ing MOi is a water-immersion Olympus LUMFLN60XW^ with NA=1.1 and a long 
working distance (WD=1.5mm). The nano-apertures are centred with respect to the 
incident beam with a nano-positioning stage. The sample is mounted on a piezo elec¬ 
tric transducer on closed-loop with resolution below O.Snm (translation range 300 /rm 
with 20-bit USB interface and noise floor of tens of picometers). Then, the centring 
is done maximizing the intensity going through the aperture when a Gaussian beam 
(vortex beam with charge / = 0) is used. When a vortex with / 7 ^ 0 is used, the 
central position is found by minimizing the intensity going through the nano-aperture. 
The interaction of the light and the centred nano-aperture occurs in the non-paraxial 
regime. Typically, the nano-aperture only allows a small part of the incoming beam 
to be transmitted. The transmitted light is scattered in all directions. This facilitates 
the coupling of light to superhcial modes of light, as it has been described theoretically 
and experimentally by many authors [243, 245, 269]. The light transmitted through 
the sample is collected by another microscope objective (MO 2 ). MO 2 is an Olympus 
MPLELNlOOx"^, whose NA=0.9 and WD=lmm. MO 2 collects the transmitted light 
through the aperture and images it on a charged-couple device (CCD) camera, which 
sits more than Im apart of MO 2 . The CCD camera is used to capture the transmitted 
intensity. 

The procedure to measure CD is the following: First, a vortex beam of topological 
charge I is created with the SLM. Secondly, the QWPi is rotated to create a LCP 
state. Then the sample is centred with respect to the beam at the focal plane of MOi 


^http://www.olympusamerica.com/seg_section/uis2/seg_uis2. asp 
'^http: / / WWW. Olympus-ims. com/en/microscope/mplfln / 
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Table 7.1: Measurements of CD(%) for three different phase singularities I as a function 
of the diameter of the nano-aperture. CD is computed using equation (7.6). 


<F(nm) 

1 = -1 

1 = 0 

1 = 1 

di = 237 

-78 ±6 

-0.130 ±0.003 

68 ±6 

d2 = 212 

-91 ±6 

-1.843 ±0.004 

84 ±6 

ds = 325 

-10 ±4 

0.596 ± 0.004 

9.8 ± 1.1 

di = 317 

-7 ±5 

1.258 ±0.005 

12 ±4 

ds = 333 

-9 ±4 

0.841 ± 0.005 

13 ±2 

de = 432 

-22.3 ±0.8 

0.596 ± 0.003 

27.1 ±0.8 

dr = 429 

-34 ±2 

-1.055 ±0.007 

37.5 ± 1.3 

dg = 433 

-46.8 ±1.3 

0.629 ± 0.004 

45.9 ±0.8 


in the manner explained before. The transmitted intensity //" is measured with the 
CCD camera, where L stands for the polarization of the beam (LCP) and I for its 
topological charge. Later, the QWPi is rotated again to change the polarization state 
to RCP. The sample is re-centred. Finally, the transmitted intensity is measured. 
Given the two measured values //" and J/^, the CD measurement given an incident 
vortex beam with charge I can be computed as: 

CD,(%) = ■ 100 (7.6) 

The results of the experiment are presented in Table 7.1 and Figure 7.4. The hrst 
column of Table 7.1 shows the size of the nano-aperture in consideration. The rest of 
columns show the measured CD using beams with different phase singularities of order 
I = —1, 0,1, respectively. As expected, due to the fact that the circular nano-aperture 
is mirror symmetric and the incidence is normal, the CD is very close to zero when 
the incident beam is Gaussian (/ = 0). The residual GD can be attributed to the 
small asymmetries of the sample (characterized in section 6.4.1) or the incoming beam. 
When vortex beams with I = —1,1 are used the situation is very different. A non-null 
GD is obtained. This value is of the order of 90% for some nano-apertures, even though 
the incidence is normal and the nano-aperture is still mirror symmetric. Because the 
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Figure 7.4: Measurements of CD(%) for three different phase singularities I as a function 
of the diameter of the nano-aperture. CD is computed using equation (7.6). Black diamonds 
are used to represent the CD obtained with 1 = 1', green circles for I = 0 and magenta hexagon 
stars for / = — 1. 

theoretical maximum value for CD given the formula (7.6) is 100%, it is fair to say 
that for certain sizes a giant value of CD has been obtained. Furthermore, it can be 
observed that there is an underlying symmetry relating the value of CDi and CD_i. 
Indeed, CDi ~ —CD_i. Hereafter the results are discussed and the mechanisms by 
which CDi ~ —CD_i % 0 are unveiled. It will be seen that symmetries considerations 
are almost enough to predict these results. However, it will be shown that not only the 
symmetries of the sample need to be taken into account, but also the symmetries of 
light. The symmetries of the target have been analysed in section 7.2. T is symmetric 
under time translations, rotation around the axis, and mirror symmetries containing 
off a plane containing the 2 ; axis. Thus, T fulhls the commutation relations given by 
equation (7.3). Now we will turn our attention to the symmetries of the light beams. 
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The electric field of the light beam incident on T can be described within the paraxial 
approximation with the complex vector: 

(7.7) 

where dp = (x+ipy)/ \/2, x and y being the horizontal and vertical polarization vectors, 
A is a normalisation constant, k = 2'k/\ with A the wavelength in consideration, and 
p, 0, 2 ; are the cylindrical coordinates. An implicit harmonic exp(—icat) dependence is 
assumed, where u = 2110 /\ is the angular frequency of light, and c is the speed of light 
in vacuum. Note that p = 1 refers to LCP and p = — 1 to RCP. Now, if we apply a 
rotation around the axis by an angle 9 to this beam, the transformed field acquires 
a constant phase: 

R,(0)ES = exp(-z(p + /)0)ES (7.8) 

That is, Ep" is a cylindrically symmetric beam (see section 4.1). Thus, E))) is also an 
eigenstates of the generator of rotations around the z axis, i.e. the z component of the 
total AM Jz- 

J,ES = (p + OES (7.9) 

The Jz eigenvalue is m = p+/. Also note that a beam with p 7 ^ ±1 would be elliptically 
polarized (linearly when p = 0 ) and would no longer be an eigenstate of J^, since a 
rotation of the field would not leave it invariant except for a phase (see Appendix A for 
more information about rotations of EM fields). The helicity operator A^ can also be 
applied to E),". If only paraxial terms are considered, then the following result yields: 

AE^^pEj (7.10) 

Hence, in the paraxial approximation, a beam E)," with p = 1 is both LCP and also an 
eigenvector of A with value 1. In the same way, a beam E)," with p = —1 is RCP and 
its helicity equals to —1. It is important to emphasize that this simple relation between 
helicity and polarization is only valid because the light beam is being described in the 
paraxial regime. When the paraxial approximation does not hold, the polarization of 
the field and helicity can no longer be so simply related with helicity, as it is shown in 


®See section 2.2 to find its differential expression. 
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chapter 2 and [54]. Now, if a mirror transformation is applied to the incident beam, 
the following resnlt is obtained: 

= exp(!a)E“ (7.11) 

where a is a phase given by the specihc mirror transformation chosen®. That is, a 
beam with flipped values of Jz and A is obtained. As shown in section 2.2, this is 
a consequence of the fact that both Jz and A anti-commute with mirror symmetry 
transformations M{z} [40]: 

= -J. (7.12) 

A4,AAf,„ = -A (7.13) 

In order to explain the experimental results given by Table 7.1 and Figure 7.4, the field 
transmitted through the nano-apertures and collimated by MO 2 E* ^ is classified with 
the parameters p and I from the incident held: p = —1,1 and I = —1,0,1. Remember 
that for the incident held p is modihed with QWPi and I with the SLM (see 
Figure 7.3). Mathematically, the held E*^j can be obtained through the use of the 
target operator T, which can be found using the Green dyadic formalism and contains 
all the relevant information about target T [243, 270]. That is, E^^ = T{Ep"}, where 
the action of T on the incident held will in general be in the form of a convolution. As 
T is simply the mathematical description of T, T inherits the symmetries of T. Thus, 
due to the cylindrical and mirror symmetries of the target and given an incident held 
Ep", the following two statements hold: First, the transmitted held E* ^ will also be an 
eigenstate of Jz with the same eigenvalue of m = p + 1. Second, two incident beams 
which are mirror images of each other will produce two transmitted helds which will 
be mirror images. The proof of both statements is given next. Due to the bijective 
properties of the exponential function, the commutation relation between T and Rz 
can be extended to the generator of rotations Jz'- 

JITJz = T (7.14) 

^Remember that is any mirror symmetry that contains the z axis. Therefore, a rotation 

angle is needed to define a particular If is written as the unitary normal vector to 

the plane is n = cos ax -I- sin ay. 
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Then, given that is an eigenvector of Jz, must also be an eigenvector of Jz 
with the same eigenvalue: 

J. [EL] = [TIELJ] = T{J. [EL] } = (P + 0T{E‘”i} = (P + i)EL (7.15) 

Then, as explained in section 7.2, given a sample which is symmetric under mirror 
transformations, its scattering matrix must commute with M{z}, i.e. [T, = 0. 

Now, given two mirror symmetric beams such as EJ," and El"p _; (see equation (7.11)), 
it can be checked that their respective transmitted helds (E*^; and Elp_J are related 
with a mirror symmetry: 


= T{EJ)} = T{exp(ia)M{2|E‘j"p_;} = M{^}T{exp(ia)E‘Tp,_J 

= exp(m)M{2}Eip 


(7.16) 


That is, given two mirror symmetric beams such as E)," and E^^^ (see equation (7.11)), 
their transmitted beams E*; = T{Ep" } and Eip _; = T{E ^-p-i} are connected via a 
mirror symmetry: E*; = exp(m)M| 2 }Elp This last result is one of the key points 
to understand the CD measurements with vortex beams presented in Table 7.1 and 
Figure 7.4. Indeed, in equation (7.6), the intensities //" and can be obtained from 
the transmitted electric held as: 

/ CO poo 

/ (7.17) 

-CO J — oo 

where the integral is taken on the plane of the detector (in our case CCD chip of the 
camera). Then, for a mirror symmetric sample, it can be proven that: 


jL/R _ jR/L 
-'i “ -'-i 


(7.18) 


The proof of equation (7.18) is done in Appendix C. Next, equation (7.18) is applied 
to prove that CDq = 0 and that CD/ = —CD_;. When Z = 0, it is easy to see that 
Iq = /(f. Substituting this in the dehnition of CD, i.e. equation (7.6), yields CDq = 0. 
However, when I ^ 0, equation (7.18) leads to //" — if- = = —(/f; — I^i), 

which implies that CD; = CD_/, in very good agreement with our measurements h 

^Note that all the experimental values are consistent with CD; = CD_;, within the experimental 
error bars. 




132 


Experiments with single nano-apertures 


Powerful as these symmetry considerations are, they still cannot explain the measured 
quantitative results of CD, nor their variation with the diameter of the nano-aperture. 
They only indicate that when I ^ 0 the two opposite circular polarizations are not the 
mirror image of each other and then the associated CD/ does not have to be zero. Here, 
the AM of light plays a crucial role again. In general, it can be observed that CD mea¬ 
surements compare the differential ratio of electromagnetic fields with opposite circular 
polarization and a difference of AM of 2 units. For example, CD/=i relates |E* and 
lE^i whose respective AM values are mp=i = 1 -|- 1 = 2 and mp=_i = 1 — 1 = 0. It 
is then interesting to observe that even though CD is usually dehned as the differential 
absorption of circular polarization states, it also measures the differential absorption 
of AM states. In the case of the nano-aperture, this is the most probable cause of the 
giant value of CD obtained in the experiments. Even though the sample is cylindrically 
symmetric, thus preserving the AM of field, input beams with different values of AM 
have very different scattering amplitudes®. This was shown and used in chapter 4 to 
enhance the ripple structure and excite WGMs. 

Furthermore, using the model of the aplanatic lens [49, 95], one can check that the 
fields at the focal plane of MOi produced by El"^ ^ and E^^"; are completely different. 
Indeed, due to their different AM content, E]^^; has a vortex beam in the z polarization, 
whereas E‘Ti_/ does not (see Figure 7.5(b))®. Consequently, their multipolar decom¬ 
positions are also very different (see Figure 4.2), surely giving rise to very different 
coupling scenarios. 

This fact can be analytically proven for spheres, where different couplings can oc¬ 
cur not only depending on the geometry of the particles but also on the beams used 
to excite these particles (see chapter 4). But actually, very similar effects have been 

®If the sample was a sphere, the scattering amplitudes would be the Mie coefficients 

®The fields shown in Figure 7.5(b) can be obtained using the expression for the Bessel beams given 
by equation (2.42). 
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Figure 7.5: a) Sketch of the addition to the set-up in order to perform a projective 
measurement of the circular polarisation states. The projective measurement is done in the 
following manner. Firstly, the metallic sample is removed from the nano-positioning system. 
Secondly, a LCP beam is focused with MOi and collimated with M02. Then, a second 
quarter-wave plate (QWP 2 ) and a linear polarizer (P 3 ) are used to select either the LCP 
(direct component) or the RCP (crossed component) of the transmitted beam. At last, the 
sample is put back on the nano-positioning device, which places it at the focal plane of MOi, 
and re-centres it with the incident beam, b) Intensity and phase plots of the two beams used 
to carry out the measurement of CD;=i, i.e. and E’j^^ In the upper row, the intensity 
and phase of the beams are shown on the back-aperture plane of the MOi. Both the intensity 
and the phase can be described with equation (7.7). In contrast, the three rows below show 
the intensity and the phase of the same beams (Ej"j^ and E’jl^ j^) at the focal plane of MOi. 
As it can be observed, even though their paraxial intensities and phases are analogous, their 
structure is completely different at the focal plane. This is a direct consequence of the fact 
that the AM of both beams differ in two units. 
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experimentally observed for other plasmonic structures [36, 247, 271, 272], 

In conclusion, a new method to induce CD has been demonstrated. Instead of tweak¬ 
ing the scattering matrix S of the structure by to changing the geometry or material 
properties, vortex beams have been used carry out the measurement. As a result, it 
has been proven that a giant CD can be induced on a single non-chiral sub-wavelength 
sample, at normal incidence. To do so, the transmission of beams with three different 
phase singularities of order I = —1,0,1 has been studied. It has been seen that the 
results can be conveyed from a symmetry perspective. In particular, it has been proven 
that even for mirror symmetric systems, CD can be induced if the LCP and RCP beams 
are not connected via a mirror symmetry. In addition, it has been shown that the in¬ 
formation carried by CD measurements has two different contributions: the differential 
scattering of different circular polarization states but also the differential scattering of 
different angular momentum states. Thus, it is expected that vortex beam-induced 
CD will be able to unveil properties of the sample which are hidden to standard CD 
measurements. Finally, it is interesting to see that in other related phenomena, such 
as molecular optical activity, the interplay between the two symmetries associated to 
helicity and AM (electromagnetic duality and rotational symmetry) are also essential 
to fully understand the problem from first principles [163]. 

7.4 Angular momentum preservation 

In this section, in order to further explore the AM scattering from nano-apertures, an 
experimental characterization of the modes transmitted through the nano-apertures 
has been carried out. As mentioned earlier, due to the cylindrical symmetry of the 
nano-apertures, both E)," and ^ are eigenstates of with value (p+l). Nevertheless, 
because the sample is non-dual, helicity is not preserved upon interaction, therefore 
p and I may not be preserved individually. This phenomenon has been detailed in 
subsection 7.2.1: it is a consequence of the duality symmetry being highly broken by 
the nano-aperture and the multilayer system [41, 54, 132, 169, 273]. In order to measure 
it, two elements have been added to the set-up to perform a projective measurement of 
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the two helicity components, and at the same time determine the topological charge I 
of the corresponding mode (See Figure 7.5(a)). The two elements added in the set-up 
are a quarter-wave plate and a linear polarizer (QWP 2 and P 3 in Figure 7.5(a)). First, 
QWP 2 transforms the circular polarization into linear. Then, due to the fact that the 
two different helicity components are orthogonal, if one gets transformed into x, the 
other one goes to y. Therefore, the linear polariser P 3 serves as an analyser: it selects 
one of the two helicity components, and removes the other one. As mentioned earlier, 
due to the fact that the duality symmetry is highly broken by the sample, the helicity 
of the incident beam E)," is not preserved. Thus, the held E*; comprises two helicity 
components, one of polarization dp and another one of polarization d_p: 

Ep,z = Al i{x, y)ap + ?/)d_p (7.19) 

where Api{x,y) and Bpi{x,y) are the complex amplitudes of the two polarizations at 
the plane of the camera. Api{x, y) is denoted as the direct component, as it maintains 
the polarization state dp and consequently the topological charge I, as AM needs to 
be conserved. The other orthogonal component is Bpi{x,y), and it will be denoted 
as crossed component. The crossed component has a polarization state d_p when the 
incident state is dp. Due to the cylindrical symmetry, the value of the AM along the ^ 
axis must be preserved. Therefore, when p changes to —p, I goes to I + 2p. That is, the 
crossed component Bpi{x,y) is a vortex beam whose topological charge differs in two 
units with respect to the incident beam E),”. This is depicted in Figure 7.6. The six 
images in Figure 7.6 are the CCD recorded images for Bpi{x,y) when the values of p 
and I of the incident beam are p = —1,1 and I = —1,0,1. Hence, the image taken with 
(p = 1, / = 1) corresponds to B\ i{x,y), and consequently its vortex has an topological 
charge /' = / -|- 2p = 1 -|- 2 = 3. However, instead of observing a vortex of charge I' = 3, 
three singularities of charge I = 1 are observed. This occurs because higher order 
phase singularities are very unstable and prone to split into hrst order singularities 
[274-277]. Thus, in the current scenario, a phase singularity of order I' will split into 
|/'| singularities of order sign{l'). This is the principal feature that is captured in the 
snapshots of Figure 7.6. They can be experimentally used to verify that the AM along 
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Figure 7.6: Crossed component \B^i{x,y)\'^ for p = —1,1 and I = —1,0,1. The values 
of p and I are the ones carried by the incident beam Given a row p and a column I, 

the image represents i{x,y)\‘^, which is a mode of light with polarization —p and a phase 
singularity I + 2p. The six far-held images have been taken using the nano-aperture dg 


the z axis is conserved within the experimental errors. The preservation of AM in the 
z direction is also observed with the direct transmitted component Api{x,y). This 
is depicted in Figure 7.7 for the same six incident beams. The intensity prohles are 
not as surprising as the ones shown in Figure 7.6, as their intensity distributions are 
qualitatively the same ones as the incident beams whose expression is given by 
equation (7.7). The preservation of AM in a cylindrically symmetric system has also 
been reported by some other authors. In a very recent experiment, Gorodetski et al. 
observed similar phenomena to the ones reported here [248] , even though their samples 
were not always cylindrically symmetric. Instead, they managed to transfer AM from 
their sub-wavelength samples to the EM held and measure it in the far held. Also 
very recently, Chimento and co-workers observed the preservation of the z component 
of the AM using an incident Gaussian beam and a 20 pm circular aperture [278]. 
Finally, Vuong et al. described how phase singularities can arise and change with the 
polarization of the incident held, in the near-held of sub-wavelength apertures [35]. 
However, to the best of my knowledge, this is the hrst time that the transmission of 
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Figure 7.7: Direct component |Ap;(x,y)p for p = —1,1 and I = —1,0,1. The values of 
p and I are the ones carried by the incident beam In contrast to Figure 7.6, p and I 

for these images coincide with the values of p and I for the incident beam. All the snapshots 
have been taken in the far-field using the nano-aperture dg- 


high AM through circular sub-wavelength apertures is measured in the far held. 
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'We are trying to prove ourselves wrong as quickly as possible, be¬ 
cause only in that way can we find progress" 

Richard P. Feynman 


8 

Experiments with single spherical particles 


8.1 Scope 

The use of light to manipulate small particles has widely spread since the seminal 
works of Ashkin in the 1970’s [220-223]. Ashkin studied the levitation of particles of 
the order of 10 /im. He showed that the radiation pressure exerted on single particles 
of that size is enough to beat gravity and Brownian motion. His seminal works led 
to the development of what nowadays is known as optical tweezers [279]. Besides, 
optical levitation was also used to measure scattering (or absorption) cross sections 
as a function of the size parameter x = 2 tiR/\ from a single scatterer. Indeed, a 
monochromatic laser can be used to trap and levitate a single large droplet of water, e.g. 
$ = 100/im. Then, the power of the laser can be increased and as a result the droplet 
starts to evaporate. This allows for measuring the variation of scattering cross section 
as a function of x [46, 280]. With the advance of laser technologies, other methods 
have been used by different authors to measure scattering efficiencies. XifrAPerez et 
al. used Fourier Transform Infra-red Spectrometer coupled to a microscope to measure 
the scattering cross section of Silicon particles in the infra-red [281]. Similarly, other 
authors have combined the use of continuum laser sources with spectral filters to make A 
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scans, while keeping the size of the single strnctnre constant [282-284]. In this chapter, 
as it was pointed ont in section 3.5, I nse cylindrically symmetric beams to measnre 
backward scattering efficiencies for single spheres. As it has been observed in chapter 4, 
cylindrically symmetric beams can be nsed to control the scattering of spheres. They 
can enhance their scattering ripple strnctnre, and in case high enongh order modes 
are nsed, they can also excite WGMs. Fnrthermore, the scattering cross section of 
spheres excited by vortex beams has, np to my knowledge, never been measnred as a 
fnnction of A yet. The only similar experimental realisations have been done at one 
fixed wavelength and only nsing LG modes of order I = —1,0,1 [285, 286]. Here, not 
only A scans will be carried ont, bnt also high AM beams snch as m = 6 will be nsed 
to stndy the interaction of symmetric light with these spheres. Even thongh the stndy 
is not fnlly complete, the resnlts presented in this chapter ratify the fact that high AM 
can nnveil some scattering resonances that are hidden nnder a plane-wave or Ganssian 
excitation. 


8.2 Experimental set-up 

The experimental set-np is schematically displayed in Fignre 8.1. Note that the mea- 
snrements presented in this chapter are done in backscattering. The experimental 
set-np allows for measnring forward scattering too, bnt it has not been nsed. The 
preparation of states is very similar to the one carried ont in the previons chapter. 
That is, a monochromatic Ganssian mode is expanded nsing two lenses (Li and L 2 ) to 
match the chip size of the SLM; then, nsing a linear polariser and a half-wave plate (Pi 
and HWP), its polarization is set to optimize the efficiency of the SLM; afterwards, 
the SLM modifies the phase of the beam to create a vortex beam of order /; the first 
diffraction order is selected using a pair of lenses and an iris (L 3 , L 4 and I); then, its 
polarization is modihed to create a state of well-dehned helicity with a linear polariser 
and a quarter-wave plate (P 2 and QWPi). Later, the light beam hits the back-aperture 
lens of MOi and is focused down to the sample. 
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Figure 8.1: Schematics of the experimental set-up. A laser beam of variable wavelength 
A = [760 — 810]nm is expanded with two lenses (Li, L 2 ) and polarized with a linear polariser 
(Pi) and a half-wave plate (HWP) to match the requirements of the SLM. The SLM gives a 
spiral phase to the beam and with the use of two lenses (L 3 , L 4 ) the first diffraction order is 
selected by an iris (I). Later, the polarisation of the beam is modified with a linear polariser 
(P 2 ) and a quarter-wave plate (QWPi). 


Here, a subtlety needs to be commented. In the experiments with nano-apertures 
shown in chapter 7, the wavelength of the laser is kept constant, therefore once the 
optical path has been aligned, the SLM is not changed. In the current case, the aim of 
the experiments makes this part a little bit more complicated. The laser used for this 
experiment is a Toptica DLpro tunable diode laser Its wavelength range goes from 
760-810 nm, and its linewidth is extremely narrow (2 ■ 10“® nm). The output power is 
over 50mW at all the different wavelengths, peaking at 785nm, where it gives lOOmW. 
Due to the 50nm range of the laser, the hologram displayed by the SLM produces 


^http://www. toptica.com/products/research_grade_diode_lasers/tunable_diode_lasers.html 
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different diffraction angles at the two boundaries of the wavelength range. That is, 
where A9x = X/Ax (see equation (6.17)). This is specially obvious 
at the focal plane of MOi. Because of the different diffraction angles, the beam hits 
the back-aperture of MOi in a slightly different position, which creates a considerable 
displacement of the focus spot in the focal plane^. Now, the SLM can correct for this 
displacement by displaying more or less 27r phase jumps. The formula to give the same 
diffraction angle for all wavelengths is: 

Ax^ = ( 8 . 1 ) 

-^ref 

where Aref is the reference wavelength for which the alignment has been done, and 
is the total number of 27r phase ramps given by the hologram at A = Aref di¬ 
vided by width of the SLM chip. I have chosen Aref to be the wavelength for which 
the provider of the SLM has made the LUT, which in this case is A = 785nm. Then, 
the alignment is done at A = 785nm with a certain Ax^’'®*^ so that the different orders 
separate and the efficiency is still high. 

The interaction between the incoming beam and the sample then takes place. The 
sample, as explained in section 6.4.2, is a coverslip which has been spin-coated with 
a solution of Ti02 particles on top. The slide is placed on a sample holder which 
is attached to a nano-positioning device. The nano-positioning device is also clamped 
onto a micro-positioning device. The micro-positioner is used to look for single spheres. 
Then the single sphere is centred with respect to the beam using the nano-positioning 
device, whose precision is O.Snm (see section 7.3). Note that MOi is a water immersion 
objective, therefore the incoming beam must go through the coverslip hrst, and then 
is scattered off the particle. 

The scattering is collected both in the forward and backward semi-spaces. Both of 
them analyse the scattering in the same fashion. Figure 8.2 shows the analyser in the 
forward direction. It can be observed that MO 2 (lOOx, NA= 0.9) collects the scattered 


^In fact, this is the basic principle of holographic optical tweezers [218] 
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Figure 8.2: Picture of the forward scattering measuring set-up. MOi focus light on the 
sample, which lies on top of a sample-holder. The sample-holder is clamped onto a nano¬ 
positioning device. The forward scattering is collected by MO 2 , which sends the light up to a 
mirror that reflects the light towards a helicity projector. The helicity projector is formed by 
a quarter-wave plate and a linear polariser. The projected state is imaged by a CCD camera. 


light and images it on the CCD camera. Before that, a qnarter-wave plate projects 
the two helicity components onto two different linear polarisations so that a linear po¬ 
lariser can select one of the two and remove the other (see section 7.4 and Figure 7.5). 
The backward scattering acquisition set-up is shown in Figure 8.3 and schematically 
in Figure 8.1. The same MOi used to focus the light on the sample is used to collect 
the scattering. Then, a beam-splitter (BS) is used to separate the scattered light from 
incident beam. The scattered beam is separated into its two orthogonal helicity com¬ 
ponents with a quarter-wave plate and then another linear polariser is used to select 
one of the two components. Afterwards, the beam goes through a lens that images the 
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Figure 8.3: Picture of the backward scattering acqusition set-up. The light scattered 
off the particle in the backward semi-space is collimated by MOi, which sends the beam 
downwards. A 45° mirror steers the light back to a horizontal plane. Then the scattered 
light is separated from the path of the incident light by a beam-splitter (BS). Afterwards, the 
light goes through some neutral dentisty filters and the helicity state is transformed using 
a quarter-wave plate and a linear polariser. Finally, a lens images the beam onto a CCD 
camera. 

sample onto a CCD camera. 


8.3 Measuring methods 

The set-up explained above allows for measuring the direct and crossed helicity compo¬ 
nents both in forward and backward scattering. However, there are a few fundamental 
differences between the acquisition of data in forward or backward scattering. Forward 
scattering encompasses the scattering of the spheres in the forward semi-plane as well 
as all the direct light which does not interact with the sphere. That is, in Mie Theory 
notation, the measured electric held is the total held E*, which is the addition of the 
scattered held and the incident one E^ (see section 3.3). Due to the size of the 
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particles used in this experiment (<l> = 950nm), the contribution of E® to the scattering 
is much smaller than the direct light Eb Thus, an intensity recording of E* with a CCD 
camera does not yield a very good signal-to-noise ratio of the scattered held, since the 
latter is blurred by the incident beam. Two possible ways of overcoming this problem 
are explained next. 

The hrst one uses the fact that the incident beams of light used in this experiment 
have a well-dehned helicity. Due to this fact, all the incident light is vertically polarised 
after it has been collected with MO 2 and transformed with a QWP. Now, the scattered 
held has two helicity components, as the sphere is not dual, i.e. E® = E® + Eb^. The 
direct component E®^ is transformed to a vertical linear polarised state. In contrast, 
the crossed component is horizontal linearly polarised after the QWP. Thus, the linear 
polariser can hlter out the light with helicity p and leave exclusively Eb^. This method, 
even though is elegant and easy to implement, does not always perform as well as one 
would wish. If we dehne as the extinction ratio of a linear polariser, then the inten¬ 
sity of the incident beam P must be smaller than ■ t~^ for this projection to work. 
This method has been tried for the set-up shown in Figures 8.2, 8.3 and P ~ P_p ■ 
therefore the experiment has not been done in this way. 

The second method is based on dark-held microscopy principles [146, 287, 288]. MOs 
manufacturers make most of the objectives so that they create an image at inhnity. 
Certainly, these are the ones used in this experiment. Then, their back-aperture plane 
contains the Fourier decomposition of the image held. Thus, a 4-f hlter system (see 
Appendix B.3) can be used to sort the Fourier content of the image held. The idea 
is that most of the incident light travels at small fcp’s, whereas the sphere scatters 
light for all angles, in particular those kp ~ k. Hence, if the light at the center of 


^In the same way as I did in the previous chapter, I denote ‘direct component’ to the scattering 
component with the same helicity as the incident light. The term ‘crossed component’ is used to refer 
to the scattered light whose helicity is the opposite to the incident field. 

Hr is defined here as the ratio between the unwanted component and the wanted one. Thus, it can 
take values between 0 and 1. 
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Figure 8.4: Schematics of the helicity content in different scattering events, a) A plane 
wave with a well-defined helicity (red) hits on an air-glass interface. The transmitted light 
preserves the helicity, whilst the reflected light flips the helicity value (blue), b) A non-dual 
sphere scatters intensity in both helicity components (red and blue). The division between 
left and write semi-spaces is merely aesthetic: both helicities are scattered in all directions, c) 
Backward crossed helicity (blue) component of the scattering of a sphere sitting on top an air- 
glass interface. The scattering encompasses the crossed helicity component in the backward 
semi-space as well as the reflection of the incident beam on the interface, d) Backward direct 
helicity (red) component of the scattering of a sphere sitting on top an air-glass interface. 
The scattering only encompasses the direct helicity component. 


the Fourier plane is blocked, only the spatial frequencies with large fcp’s will be im¬ 
aged by the CCD. Those Fourier components mostly correspond to the scattering of 
the particle. Even though this method can potentially give good signal to noise ratios, 
we opted for the method presented in the next paragraph because it was more versatile. 

The scattering can also be measured in the backward semi-space. And this is the way 
the measurements presented in this chapter have been carried out. An intuitive model 
behind the measurements in backward scattering is presented in Figure 8.4. Remember 
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that the sample that is going to be used to carry out A scans with different optical 
vortices is the one depicted in Figure 6.11 - single Ti02 particles with <l> = 950nm on 
a glass coverslip. Then, Figure 8.4(a) sketches the interaction between an incoming 
plane wave (grey straight arrow) with a well-defined helicity (red spinning arrow) and 
a glass-air interface. It can be seen that a small fraction of light (generally it is about 
4% of the incident intensity) represented by a dashed arrow is reflected back. The rest 
of the 96% of light goes through and maintains its helicity value^. Now, note that the 
4% of the light that is reflected changes its helicity (blue spinning arrow). The interface 
acts as a mirror. Without loss of generality, the ^ axis can be placed in the direction 
normal to the interface. Then, the ^ component of the incident k is flipped, and 
ky maintained. The same mirror symmetry applies to the polarisation vectors. That 
is, z ——z, and x —)■ x, y —y. That implies that the polarisation in the real space 
(computed with respect to some external reference frame in the lab) is maintained, but 
its projection to the linear momentum direction yields an opposite helicity [165]. The 
same conclusion can be reached using the anti-commutation relations between A and 
Mn (see section 2.2 and equation (7.13)). 

Note that the previous explanation is only valid if the incident field is a plane wave 
with a well-defined helicity propagating in the direction normal to the interface. In the 
experiment, the field that impinges the coverslip is a fully vectorial focused beam with 
many different plane waves components, and therefore these considerations do not fully 
apply. In particular, the reflected light also has some components with the same helicity 
value as the incident beam, as the symmetry arguments why the helicity value must flip 
only apply for plane waves travelling at normal incidence to the interface (see section 
5.2.) Then, the reflected beam is computed applying the reflection Fresnel coefficients 
to each of the s' and p waves of the incident field and then adding them all together [52]. 

Figure 8.4(b) shows the fact that a non-dual sphere scatters light in both helicities. 

®Even though the glass is not dual, the helicity must be preserved in the forward direction due to 
the cylindrical symmetry of the interface (see section 5.2 for a more detailed explanation). 
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The division between left and write semi-spaces is pnrely aesthetic. Both helicities 
are scattered in all directions in space. The red spinning arrows represent the direct 
helicity component, while the bine ones represent the crossed component. It can be 
observed that the polarisation of a plane wave component scattered in the forward di¬ 
rection with helicity p is the same one as the polarisation of a plane wave with helicity 
—p scattered in the backward direction. This property is the key why the acqnisition 
of data has been done in backward scattering. 

Figure 8.4(c) shows the excitation of a sphere with a beam with helicity p (red), its 
reflection on the glass-air interface (dashed grey arrow), and the scattering crossed 
helicity component of the sphere in the backward semi-space. Now, because both the 
reflection from the interface and the scattering crossed component have the same he- 
licity, they all acquire the same horizontal linear polarisation after they are projected 
to a linear basis by the QWP in the backward set-up (see Figure 8.3). Similarly, the 
direct component of the scattering also gets projected onto the linear basis. However, 
its state becomes vertically polarised. Note that in this approximate model, all the 
light reflected off the glass is in the crossed helicity component. Therefore, if the linear 
polariser filters out the crossed component (or horizontal polarisation), all the remain¬ 
ing light will exclusively be from the direct scattering component of the sphere. This 
is depicted in Figure 8.4(d), where it is shown that all the light with helicity p on the 
backward semi-space comes from the scattering of the sphere. This is the way the cross 
section curves as a function of A have been measured. That is, instead of measuring 
all the scattered intensity, only the direct component of scattering in the backward 
semi-space has been recorded. The used measuring protocol is the following one: 


1 . Alignment at A = 785nm . As mentioned previously, the system is aligned at 
A = 785nm. That is, the SLM is aligned with the two MOs and the camera so 
that all the different orders that will be used in the experiment (Z = —5, ..,5) are 
properly focused and collimated. 
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Figure 8.5: Forward imaging of the sample. Left: snapshot of a particle which is not 
isolated enough to perform single scattering experiments with it. Right: snapshot of a single 
sphere. The experiments have been carried out with this particle. There are no particles 
within 20 /rm. All particles are made of Ti02 and their = 950nm. 


2. Finding a suitable single sphere . The sample is moved with the nano¬ 
positioning device to single out a Ti02 sphere on the glass substrate. In order 
to see that in screen of the computer, the focusing objective (MOi) is moved out 
of focus so that the illumination is constant in the held of view of MO 2 . Then, 
the position of the particles can be tracked with the CCD camera in the forward 
semi-space. It is important that the selected sphere is isolated enough, so that 
no multi-scattering phenomena take place. Figure 8.5 shows two CCD snapshots 
of single particles. Clearly, the single sphere on the left is not far enough from 
other spheres, therefore it has not been used. In contrast, the single sphere on 
the right snapshot is isolated enough. This is the particle that has been used for 
the experiments presented in this section. 

3. Bagkground measurement at every a . Once a sphere has been singled 
out, MOi is positioned to focus light on the particle. At the same time, MO 2 is 
moved mm’s away from the focus, to avoid reflections from its lenses. Then, the 
particle is laterally moved tens of microns away from the position where the beam 
is focused. At this point, a snapshot of the beam is taken. This snapshot will 
be referred to as the background. The idea is that, this reflected light cannot be 
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removed from the scattering measurements®. Thus, it will be used to normalise 
the scattering of the particle. The background measurement depends on the 
incident mode, consequently a background snapshot will be taken for each mode 
Ep”, with p = —1,1 and I = —5,.., 5. The process is repeated at every different 

A. 

4. Scattering measurement at every A . The sphere is placed back at the 
centre of the beam. Now, centring the sphere with respect to the beam is not an 
easy task to do. In fact, it gets more and more cumbersome as the topological 
charge of the incident beam gets larger. The main reason why high optical 
vortices are hard to centre when they are focused was already mentioned in 
section 7.4. An optical vortex of charge I breaks and it gives rise to |/| vortices 
of charge signal) [274-277]. Then, |/| singular points can be found in the optical 
plane, which slightly break the cylindrical symmetry. In order to minimise this 
symmetry-breaking for every A, E *^"_3 (for p = 1 helicity modes) and 3 (for 
p = —1 modes) have been used to centre the particle with respect to the beam. It 
can be observed (see Appendix D) that a multipolar held in the backward 
semi-space is identical to the multipolar held in the forward semi-space. 

That is, 

(J < 0) = A'S. (z > 0) (8.2) 

Now, when is collimated and made paraxial, it yields an optical vortex of 
charge I' = m + p, where m = p + l for beams of the kind E),", giving I' = 2p + l. 
Then, in the direct component of backscattering, E 3 T 3 and E’A'^ 3 yield vortices 
of charge /' = 1 and I' = —1 respectively. Hence, if the CCD image obtained in 
backward scattering is a cylindrically symmetric beam with |/| = 1^, it means that 
the incident beam (E'^'As or E‘A\ 3 ) is properly centred with respect to the spherical 
particle. Then, once the particle has been centred, ten snapshots are taken of 
the scattering of the sphere for each combination oi p = —1,1, I = —5, ..,5. 

®As explained in Figure 8.4(d), the background would be null if the incident field was a plane wave 
with a well-defined helicity propagating in the direction normal to the glass interface. 

^See Figures 7.6(1 = —l,p = 1) or 7.7(1 = l,p = —1) to see an example of a scattered vortex beam 
with |1| = 1. 
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Figure 8 .6: as a function of A. The incoming beams are of the kind where 

I = —5, ..,5. Note the resonant behaviour of the modes with I = ±1 and I = 2 at 785nm. 
/norm jg defined as the back-scattered intensity when the sphere is centred with respect to 
the incident beam over the intensity of the background. Both intensities are recorded with a 
CCD camera. 


5. CCD NOISE RECORDING . To compensate for the inherent noise of the CCD 
camera due to electronic noise and others, thirty snapshots without light on the 
camera have been taken for each of the exposures used to record the background 
and scattering of the particle. In this way, an image of background noise is 
obtained, and it will be subtracted to all the rest of snapshots (background and 
particle). 

8.4 Wavelength scans 

In the previous section, I have described the protocol to measure the backscattering of 
single spheres projected into their direct helicity component, given an incoming beam 
Ep" and a wavelength A. In this section, the collected data are shown and interpreted. 
First, A scans for E)," with p = 1 and I = —5, ..,5 are presented in Figure 8.6. The L 
values in the legend correspond to the I value of the incident mode used to record that 
plot. The horizontal axis shows the excitation wavelength A, and the vertical axis a 
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Figure 8.7; Resonant behaviour of J"™™ as a function of A. The beams depicted in the 
figure are where I = —2,—1,0,1. The coloured lines highlight the resonant behaviour of 
the particle around A = 785nm. The colors are the same ones as the ones used in Figure 8.6. 
The rest of points, in a black dashed line, are left in the plot for completeness. It is clearly 
observed that the resonant behaviour of the particle is not observed with a Gaussian (/ = 0) 
excitation. The use of vortex beams is crucial to unveil the resonance. 


normalised intensity /"°™. /norm jg computed as the scattered intensity divided over 
the background. The scattered intensity is computed when the single sphere that we 
want to study is centred with respect to the incident beam. In contrast, the back¬ 
ground is computed to normalise the scattered intensity, i.e. it is the reflection of the 
beam on the glass when the sphere is tens of microns away from the beam (see previous 
section). Consequently, the value of /“°™ is always /°°''™ > 1 . Because 10 pictures 
were taken for each of the scattering events, an average picture is computed®. Then, 
a region of interest in the average picture is dehned. Finally, the scattered intensity is 
computed summing all the individual values of the pixels in this region of interest of 
the average picture. The same region of interest is dehned for the background picture 

®As mentioned before, the noise of the CCD camera is subtracted at every single snapshot (back¬ 
ground and scattering). 
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Figure 8 .8: Standard deviation crAr(%) of the recorded scattering as a function of the 
mode and the excitation A. The incoming beams have helicity p = —1. 


and the values of the pixels in that region are summed up. That gives rise to the 
value of the background intensity. Then, as explained before, is computed as the 
ratio between the scattered intensity over the background. Hence, every single point 
in Figure 8.6 is a dimensionless number computed as the ratio of two intensities, i.e. 
jnorm jg efficiency factor (see section 3 . 5 ). 

These 10 different pictures are also used to compute statistical errors on the measure¬ 
ments. These are shown in Figure 8 . 8 . Except for one value giving a standard deviation 
of crAr(%) 80, all the rest of errors are of the order of aAr(%) ~ 10. In order not to 
make the intensity plots even more complicated, these errors have not been displayed 
as error bars in Figure 8 . 6 . 
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Figure 8.9: as a function of A. The incoming beams are labelled as a function of 

their AM content, m = I + p, where p = 1 in this case. 


Figures 8.6 and 8.7 show very interesting results. Obviously, the A scans with the 
different modes of light yield different curves. Then, it can be observed that the sphere 
behaves very resonantly around 785nm. 

In fact, the resonance is hidden under the Gaussian excitation (/ = 0), and the use 
of vortex beams is crucial to unveil it, as predicted by the theory in chapter 4 (see 
Figure 4.5). Also, as expected, the Gaussian mode (/ = 0) is the one that yields a 
larger intensity. However, the mode with / = — 1 also gives a comparable amount of 
scattering. It is important to note the huge asymmetry between the modes with I = 1 
and / = — 1. The mode with / = — 1, scatters much more light than the one with / = 1. 
This is due to their different AM content, m = I + p. Now, because a beam of = m 
can only excite Mie resonances of order j > m, it is convenient to re-label the incident 
beams in terms of their AM content. This is displayed in Figure 8.9. 

The A scans obtained for E)," with p = —1 and I = —5, ..,5 are presented in Fig- 
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Figure 8.10: as a function of A. The incoming beams are of the kind where 

/ = -5,..,5. 


ure 8.10. Again, their error bars are not plotted along with the average measnrements. 
Instead, they are depicted in Fignre 8.11, yielding standard deviations smaller than 
14%. A close look at Figure 8.10 shows that the results are very similar to the ones 
presented in Figure 8.6. In fact, it is seen that there is an underlying symmetry relating 
the results for EJ," and The underlying symmetry relates the scattering of 

vortex beams with an topological charge I and a helicity p with the scattering of modes 
with topological charge —I and a helicity —p. This is the same symmetry relation 
that yielded CDg = —CD_g in section 7.3: mirror symmetry. In fact, the sample used 
in this experiment (glass substrate with a single sphere on top) has exactly the same 
symmetries as the sample used in chapter 7 (circular nano-aperture in a metallic him). 
Those symmetries are mirror symmetry with respect to any plane containing the z axis. 
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Figure 8.11: Standard deviation cTAr(%) of the recorded scattering as a function of the 
mode and the excitation A. The incoming beams have helicity p = —1. 


rotations along the ^ axis, and time translations (see section 7.2 for a detailed explana¬ 
tion of the symmetries). The mirror symmetric scattering of the particle is highlighted 
in Figures 8.12, 8.13. In Figure 8.12, the scattering produced by modes with m < 0 
for beams with p = 1, and m > 0 with p = — 1 is compared. The resemblance of 
and /“““^(A) is apparent. Similarly, Figure 8.13, portraying the comparative 
scattering of modes with m > 0 for p = 1 and modes with m < 0 for p = —1, also 
shows a great match between the results using the two different families of modes. 
Actually, one of the few signihcant differences is that in the plots for p = 1, there is 
a drastic decrease of intensity for the modes with I = 0, —2, —3 at A = 800nm. The 
fact that this decrease was not observed in its mirror symmetric scattering counterpart 
leads to believe that it is an experimental human error. 


It is clear that the AM of light plays a crucial role in GLMT. As it was predicted 
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Figure 8.12: Comparative scattering of mirror symemtric beams. Left: the incident 
beams have p = 1 and m = I +p < 0. Right: incident beams have p = —1 and m = I +p > 0. 


in [107, 121], the AM of light can be used to unveil some hidden resonances in the 
Gaussian (or plane wave) excitation. Furthermore, it is clear that symmetry consider¬ 
ations allows for qualitative predictions of the results. In this case, the mirror symmetry 
of the sample allows for a prediction of the scattering intensity of beams _i once 
the scattering of E)," modes has been found. As future work, I would like to apply the 
same technique to structures with spherical symmetry, instead of spheres in contact 
with an interface. In order to do that, the method to embed particles in a polymer 
needs to be improved so that the faces of the sample are rather parallel and the scat¬ 
tering from the polymer is null. Another option would be to trap the particles with 
an optical tweezers set-up, which would be easy to implement in my current set-up. 
In a similar direction, I would like to carry out similar experiments with particles of 
different sizes. In particular, I would like to excite WGMs with direct light using the 
technique described in section 4.4. For that matter, particles of large sizes as well as 
higher order optical vortices need to be used. 

8.5 Characterization of crossed helicity modes 

To conclude with this chapter, some of the GGD images that have been used to compute 
the intensity plots in Figures 8.6-8.13 are shown. Because the experiments have been 
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Experiments with single spherical particles 



Figure 8.13: Comparative scattering of mirror symemtric beams. Left: the incident 
beams have p = 1 and m = I +p > 0. Right: incident beams have p = —1 and m = l+p < 0. 


done with the particles centred with respect to the incident beam, the images used to 
compute the background normalisation and the ones recording the scattering of the 
particles share the same features. That is a consequence of both system having the 
same symmetries. The basic difference between the images is the magnitude of the 
intensity. Next, the background CCD images are shown. As explained before, each of 
the snapshots has been taken using an incident beam of the kind Therefore, the 
images are classihed in terms of the charge of the optical vortex created by the SLM (/) 
and the helicity of the incident beam p. Note that, because of equation 8.2, the images 
show the spatial shape of the crossed helicity component in the forward direction. That 
is, given an incident beam of Jz = m = I + p, the captured image has \l + 2p\ optical 
vortices of charge sign{l + 2p). This fact is easily observable from the pictures when 
\l + 2p\ < 3, yet it is difficult to appreciate when the number of singularities exceeds 
3. Finally, note the mirror symmetry between the image prohles. That is, the images 
obtained for EJ," are practically identical to the ones obtained with E^^ To the best 
of my knowledge, the experimental characterisation of crossed components of vortex 
beams with a well-dehned helicity has not been published anywhere. 
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Figure 8.14: Background scattering intensity given an incident beam with p = ±1 
and I = —5, —4, —3. The sphere is about 20;um away. 
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Experiments with single spherical particles 



Figure 8.15; Background scattering intensity given an incident beam with p = ±1 
and I = —2, —1, 0. The sphere is about 20/im away. 
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Figure 8.16: Background scattering intensity given an incident beam with p = ±1 

P,L 

and I = 1,2, 3. The sphere is about 20pm away. 
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Experiments with single spherical particles 



Figure 8.17; Background scattering intensity given an incident beam with p = ±1 
and Z = 4, 5. The sphere is about 20;um away. 







"Education must provide the opportunities for self-fulfilment; it can 
at best provide a rich and challenging environment for the individual 
to explore, In his own way" 

Noam Chomsky 
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Conclusions 


As technology advances, science is challenged with new problems and vice-versa. In 
particular, the development of dielectric photonic technologies in the 1980’s made it 
clear that photonics could compete with electronics to encode, transmit and retrieve 
information. However, the current drive to shrink down our devices has enabled elec¬ 
tronics to clearly outperform photonics. Nanophotonics is a young field, and there is 
still a lot of room for improvement. In this sense, the work presented in this thesis 
is done from a different approach to the mainstream research in nanophotonics. Most 
of efforts in overcoming the diffraction limit of light have been done in the field of 
plasmonics, where the materials and the geometry the nano-structures are made of 
are crucial to control light-matter interactions at the nano-scale. Here, light-matter 
interactions are described in terms of their symmetries. The fact that light can be 
experimentally modified so that it is symmetric under many different transformations 
is crncial in the process. Within this framework, the main findings of this work are the 
following ones: 

Systematic study of EM modes in terms of their symmetries. In chapter 
2, a systematic stndy of six different basis of transverse solutions of Maxwell eqnations 
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is done. Moreover, the relations between the different basis are also given. Then, in 
chapters 4,5,7,8, different problems are faced. Depending on the symmetries of the 
problem, one of the basis is chosen over the others to describe the light-matter inter¬ 
action. These studies have been partially published as part of [54, 76]. 

Effect of the aplanatic lens on the multipolar content of a beam. The apla- 
natic transformation of a lens, which is described in section 2.6, is used to modify the 
beam shape coefficients of a beam. That is, the NA of an aplanatic lens is used to 
control the multipolar content of a beam. This effect is part of the findings presented 


in [107]. 


Enhancement of ripple structure. In Mie Theory, the scattering efficiency is com¬ 
puted as an infinite summation of Mie coefficients: 



(9.1) 


Due to the infinite summation, the resonances of the high order Mie coefficients are 
only perceived as little ripples in the cross section. These ripples have a very low Q 


factor (5 parameter in section 4.4). In chapter 4, it is proven that the use of LG beams 
(and in general any cylindrically symmetric beam) allows for the enhancement of the 


ripples. In fact, their Q factors can be increased by two orders of magnitude even using 
relatively low-order LG modes such as LGg^o- The downside of it, there is much less 
scattering, as all the contributions of the low orders modes are suppressed from the 
scattering. These results were published in [121]. 

Excitation of single multipolar modes. In chapter 4, a technique to excite single 
multipolar resonances regardless of their size and index of refraction is explained. In 
particular, the technique can be used to excite WGMs. The excitation of single modes 
is achieved in two steps. First, a cylindrically symmetric beam with a large value of 
Jz = l + p is paraxially created, with p being the helicity of the beam and I the order of 
the phase singularity that carries. To increase the cross section, the paraxial beam is 
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focused onto the sphere with a MO with a large NA. Then, the wavelength of the laser 
A and the radius of the sphere R are chosen so that the following expression holds: 
R ~ \{l + p)/(nr)/(27r), where f{nr) is a function of the relative refractive index that 
can be computationally checked. For example, when = 1.5, /(n^) ~ 0.8. This 
technique was explained in detail in [107]. 

Generalization of Kerker Conditions. The Kerker conditions were demonstrated 
by Kerker and co-workers in 1983. The first Kerker condition (Kl) states that a sphere 
with e = /i has zero backward scattering. Then, K2 states that a small particle with 
/i 7^ 1 can have zero forward scattering if oi = —hi. In section 5.2, it is demon¬ 
strated that Kerker conditions can be generalized to cylindrical objects of arbitrary 
size. In fact, it is proven that Kl is a specific case of a system symmetric under duality 
transformations and rotations along the z axis. Then, in a very similar fashion, K2 
is generalized to anti-dual cylindrically symmetric systems. That is, if a scatterer is 
anti-dual and symmetric under rotations along the ^ axis, then its forward scattering 
must be 0. A detailed proof of this generalization was published in [150]. 

Inducing dual or anti-dual behaviours in dielectric spheres. It is very challeng¬ 
ing to find dual materials in nature. However, some materials can behave as such under 
certain conditions. In chapter 5, an analytical method to induce duality in dielectric 
spheres is shown. The method can be applied to any geometry and material, but it 
becomes particularly efficient for spheres, as it can be computed analytically. In fact, 
using the same formalism, the anti-dual behaviour of spheres can also be probed. Then, 
the influence of the AM of light on this process is studied. It is demonstrated that the 
AM of light adds another degree of freedom to the previously discussed method. Thus, 
in principle, arbitrarily large spheres can behave as dual or anti-dual when high order 
AM modes are used. This method as well as some practical applications were shown 
in [132]. 
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CD with non-chiral sample. CD is a widely used technique in biology and chem¬ 
istry. Recently, it has gained a lot of interest among the metamaterials and plasmonics 
community. In chapter 7, a symmetry-based method to induce CD in a non-chiral sam¬ 
ple is shown. It is observed that a circular nano-aperture can induce a giant CD when 
it is excited with vortex beams. The key reason why this phenomenon arises is that 
the two vortex beams used in the measurement are not mirror symmetric with respect 
to the other. Furthermore, it is explained that CD not only compares the differential 
absorption of left and right circular polarisation, but also the differential absorption of 
AM momenta states differing in two units. These results are presented in [289]. 

Transmission of optical vortices throngh sub-wavelength nano-apertures. 

In section 7.4, the transmission of optical vortices through a sub-wavelength circular 
nano-aperture is demonstrated. The decomposition of the transmitted light into its 
two helicity components shows that optical singularities arise in scattering processes 
due to the fact that duality symmetry is broken. Similar experimental effects have 
been presented in [169, 249]. 

Experimental excitation of scattering resonances with optical vortices. Some 
of the predicted effects in chapter 4 have been experimentally verihed in chapter 8. In 
particular, it has been seen that vortex beams can unveil scattering resonances for 
dielectric spheres. That is, it has been observed that a relatively flat scattering cross 
section (as a function of the size parameter x) can be turned into a resonant one when 
the appropriate vortex beam is used to excite the single sphere. The effect has been 
corroborated with two sets of beams with different helicities. As a result, the same 
mirror symmetric properties observed in chapter 7 have been also measured. 

As it has been mentioned before, the hndings of this thesis advance the control of 
light-matter interactions at the nano-scale using a new symmetry-based formalism. 
For example, thanks to the method put forward in section 4.4, Q factors of the order of 
Q = 1- 10® are accessible with spheres of radius R = 1.9/rm. Also, as shown in chapter 
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7, a new wealth of information from nano-samples can be extracted carrying out a 
CD measurement with vortex beams. Figure 7.4 in section 7.3 depicts samples whose 
diameter sizes differ in only 100 nm retrieve CD values that differ in more than 70%. 
Another effect with direct consequences in nano-science is the induced duality described 
in chapter 5. The design of new metamaterials requires a great control of the building 
cells. In that sense, Figure 5.2 can be used as a design map for new highly directional 
metamaterials, as dual spheres have zero-backscattering. The applications described so 
far are direct implications of the results presented in the previous chapters. Neverthe¬ 
less, the findings in this thesis open up new possibilities in other fields. For instance, 
there is a growing interest in controlling the magnetic component of light. Different 
groups have recently fabricated structures that can enhance the magnetic component 
of light in certain regions. The semi-analytical techniques described in chapter 4 allow 
for scattering control of spheres. In particular, the technique can be applied so that the 
scattering is purely magnetic. Then, the analysis of the dual behaviour of a scatterer 
could result to be very useful to achieve very interesting regimes of optical forces. The 
fact that an anti-dual scatterer does not scatter light in forward and a dual scatterer 
does not scatter light in the backward direction could be linked to pushing and pulling 
forces. Also, the excitation of single high order multipolar modes could be used to 
carry out side-band cooling in quantum optomechanics experiments. The idea is that 
the excitation of WGMs with direct light can be easily added in an optical levitation 
set-up, which is a very optimum way of cooling down the centre of motion of a particle. 
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A 

Rotation of electromagnetic fields 


In order to compute the rotation of EM fields, the following formulae are necessary. 
Given a EM field t), a rotation of magnitude (p along an axis n is computed using 
equation (2.14): 

RM [^(r, t)] = Mfi((p) • ^p{M7\p) . r, t) (A.l) 


where is given by the following expression 


Mfi((p) 


cos 93 
sin ip 



— sin p 0^ 
cosp 0 

0 1) 


(A.2) 


if the cartesian axis are rotated in a way such that n = z and p > 0 corresponds to a 
dextrogyrate rotation. Now, this is not the only way of computing rotation matrices. 
The most usual way is using Euler angles. That is, the reference frame is considered 
clumped in the origin and then any rotation is expressed as the product of three 
rotations [39]: 


M(a,/3,7) = M,(«)My(/?)M,(7) 


(A.3) 
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Rotation of electromagnetic fields 


where Mz(a) and My(/3) are given by: 



cos a 

— sin a 

0^ 


COS (3 

0 

sin (3 

M,(a) = 

sin a 

cos a 

0 

My(/3) = 

0 

1 

0 


\ 0 

0 

ij 


1 — sin f3 

0 

cos f3 y 


(A.4) 


In this thesis, especially in chapter 2, the operator of rotations R(0fc, ^k) has been used. 
This operator rotates EM fields following equation (A.l). Then, the rotation matrix 
can be obteined as M(a,/9,7) = Mz(0fc)My(6'fc)Mz(O), which yields the following two 
rotation matrix: 




^cos 9k cos (pk — sin (j)k sin 9k cos 0^^ 
cos 9k sin (pk cos (j)k sin 9k sin (pk 
^ — sin 9k 0 cos 9 k y 


(A.5) 


^cos 9k cos (pk sin (pk cos 9k — sin 9k'^ 


M-\cPk,9k) = 


■ sm I 


cos (pk 0 (A.6) 

^sin 9k cos (pk sin (pk sin 9k cos 9k y 
As an example, the rotation of a plane wave is shown. The equation (2.32) will be 
proven step by step: 


R{(pk,9k)^exp{ikz) = M{(pk,9k)icexp{ikz-M ^{(pk,9k)r) (A.7) 


i/cz ■ M ^{(pk,9k)T = ik {sm9kCos(pkX + sm9ksm(pky + cos9kz) = ik ■ T (A.8) 


M{(pk,9k)x = 


^ cos 9k cos (pk ^ 
cos 9k sin (pk 
y - sin 9k 


= -P 


(A.9) 













B 

Holography 


This Appendix gives supplementary information to fully understand the holography 
techniques used in this thesis. 


B.l Fabrication of phase-only CGH 

Hologram design 

Designing the phase hologram is as simple as computing the desired phase as a function 
of the transverse coordinates x, y. In general, the computed phases will be additions 
of phase singularities of order m and oblique phases. An oblique phase represents 
the tilt that a prism would imprint on the beam, where is related to the apex of the 
prism [45, 268]. The larger kx is, the more the beam will be tilted. In fact, kx = ‘2'xjlS.x^ 
where Ax is the width of a 27r phase ramp (see Figure B.l). Then, the superposition 
of an oblique phase and a phase singularity yields a so-called pitchfork hologram (see 
Figure B.l), which is used to create tilted vortex beams [268, 290]. 
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GA.r 



Twisting phase Oblique phase Pitchfork hologram 


Figure B.l: A pitchfork hologram is obtained adding a twisting and an oblique phase. 
The width of six 27r phase ramps is highlighted. 


Hologram print 


Once the hologram has been digitally created, it has to be printed. In order to have 
a good efficiency, CGHs need to have a smooth gray-scale to go from black (27r phase 
shift) to white (0 phase shift). This fact depends on many factors, and the first of them 
is the printer. The chosen printer needs to provide a good gray-scale. Generally, inkjet 
printers not only have better gray-scales than laser printers, but they are also cheaper 
[291]. Thus, all the holograms made during my thesis were done with an inkjet printer 
Canon iP4850 Ink. Also, printers can print both in black and white or color mode. The 
difference between the two modes is that when printing in black and white, the gray 
tones are done by printing black dots. If the gray tone is darker, then the black dots 
are more packed; otherwise, they are more separated. Therefore, the fact that this is 
seen as gray is an optic effect, but it is not truly gray. Colour printing goes through all 
the gray tones by mixing colors in a RGB (red green blue) basis. I printed holograms 
in both conhgurations to carry out the characterization presented in section 6.2.3, but 
their efficiency was not seen to clearly depend on the printing technique. 
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20 |jm 

200 jjm 

'k 

Figure B.2: Sketch of a transverse section of a photographic film. The emulsion is 
sandwitched between the supercoat and adhesive layer, which separates the emulsion from 
the transparent base. The thickness of the emulsion is about 20 ^m, and the thickness of the 
base is about 200 ^m. 

Photography of the hologram 

Once the hologram has been printed out, a picture needs to be taken. At this stage, 
several factors can play a very determinant role to obtain a good spiral thickness d{(j)). 
The hrst one is the photographic him. A sketch of a photographic him can be found 
in Figure B.2. From top to bottom, it is formed by a supercoat, which protects the 
emulsion; an emulsion made of a photo-reactive material, which is generally composed 
of silver halide (AgBr) grains dispersed in a gelatine substrate; an adhesive layer so 
that the emulsion does not how oh the him; and a transparent base, which supports the 
rest of the him. Now, the key parameter for photographic hlms is the lines pairs per 
millimetre (Ip/mm). One line pair is a pair of black and white lines next to each other. 
The larger the number of Ip/mm is, the hner details the picture can resolve. There is 
a myriad of photographic hlms in the market, but the vast majority of them have a 
resolution of 100 Ip/mm. In contrast, a German company called Gigabithlm^ produces 
high resolution hlms with up to 900 Ip/mm (see Figure B.3). These are the ones that 



^ http : / / WWW. gigabitfilm. com / 
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Figure B.3; Schematics of a photographic emulsion for a) Generic photographic film with 
100 Ip/mm b) Gigabitfilm with 900 Ip/mm. The generic one contains randomly located silver 
halide grains of different sizes, whilst the Gigabitfilm keeps both the size and the position of 
the grains in a much ordered manner 


Table B.l; Gommon lens apertures. 


larger aperture / more light 

1.7 

2.8 

4 

5.6 

8 

11 

16 

smaller aperture / less light 


have been used, granting a higher efficiency for the holograms. Another parameter 
that plays a very signihcant role is the camera. In order to make holograms, a Minolta 
manual reflex SR-T 101 camera was used. It is important that the camera is reflex, 
so that what is seen by the viewhnder of the camera is seen in the same way by the 
objective (and him). Then, two parameters of the camera system need to be adjusted: 
shutter speed and aperture. The shutter speed, usually given in fractions of seconds, 
is the amount of time that the camera is permitting light from the lens to hit the him 
or digital sensor. A fast shutter speed is good for “freezing” the action of a subject in 
motion. A slow shutter speed is good for purposely blurring the motion of a subject. 
After carrying out some batches of holograms, I have seen that holograms generally 
turn out better for relatively large exposition times (see section 6.2.3). The aperture 
(or f-stop) refers to the size of the opening in the lens of the camera when a picture 
is taken. With aperture, smaller numbers (ex. f/4) mean a larger opening. Large 
numbers (ex. f/16) mean a small opening. Any time the amount of light is increased 
or decreased by a factor of two, the exposure changes by “1 stop”. Table B.l shows 
the most common apertures. The aperture size of the lens is also related to the depth 
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Figure B.4: Photography set-up. The camera is mounted on a tripod at a certain dis¬ 
tance of the printed holograms, which are stuck to a wall (or cupboard in this case). The 
illumination is white and uniform. 

of field. A large aperture causes a small depth of field, and a small one causes a large 
depth of field. The depth of field of an image is related to the distance between the 
nearest and farthest objects in the scene that appear acceptably sharp in the image. 
In section 6.2.3, the relation between f-stop, exposure time and hologram efficiency is 
shown. Also, the picture needs to be taken free of vibrations. In order to achieve that, 
the camera needs to be clumped on a stable tripod. Moreover, in order to even reduce 
more the vibrations, a trigger has been used to shoot the pictures. The last parameter 
that play a role on this step is the illumination of the printed CGH. Even though its 
dependence has not been numerically quantified, a uniform white illumination on the 
hologram is necessary to obtain a good photography. 

Developing process 

The development is the most important step in order to obtain efficient holograms. 
First, the film needs to be removed from the camera and placed in a developing tank, 
where the film is protected from light. This process needs to be done in a dark chamber. 
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or somewhere where light is over 690 nm. Then, the whole developing process can be 
started, which is divided into the following steps: 

1. Development . The developer transforms the latent image^ to metallic silver. 
Lots of different developers are available in the market, but I made my own one, 
as both the dims and the application that they are used for (phase holograms) 
are not conventional. The main chemical needed for a developer is dihydroxyben- 
zene or Catechol. This is mixed with NaOH and pure water using the following 
quantities: 

• 0.6g catechol in 5ml of pure water. 

• 0.32g NaOH in 8ml of pure water 

• 300ml of pure water 

To keep its developing properties, the mixing has to be made before being poured 
into the developing tank. It is important to continuously stir the developer while 
it is in contact with the film. This stage lasts 6min. 

2 . Stop bath . Once the negative has been developed, the tanks needs to be quickly 
emptied and the stop bath needs to be poured. The stop bath is an acid mixture 
that halts the action of the developer by drastically changing the pH of the 
medium surrounding the film. The developer is a basic solution, and the stop 
bath is acid. It only needs to be stirred in for 30s. It is made of: 

• 10ml acetic acid 

• 320 ml pure water 

3. Bleaching . Bleaching is a crucial step in this fabrication process. Together with 
the fixing, it transforms the him from an amplitude to a phase-only hologram. 
There are different kinds of bleaching. Hariharan studied the efficiency of the 
different kinds of bleaching and reached the conclusion that the rehalogenising 

^The latent image is an invisible image produced by the exposure to light of a photographic film. 
The chemical structure of the silver halides forming the emulsion changes, but those changes are not 
visible. 
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bleach allows for the greatest efficiency of phase-holograms [174, 182, 183]. The 
principal difference between the rehalogenising bleaching process and the others is 
the addition of a rehalogenising agent, which in Hariharan experiments as well as 
mine is potassium bromide (KBr). The products used to make the rehalogenising 
bleach are the following ones: 

• 0.3g potassium dichromate 

• 5g potassium bromide 

• 0.63ml sulfuric acid 

• 300ml pure water 

A rehalogenising bleach converts the developed silver image into silver halides 
and modihes the gelatin matrix where the silver image has been formed. The 
new-created silver halides will be removed a posteriori in the hxing process, but 
the modihed gelatin will stay. The sulphuric acid turns the medium acid so that 
the different chemical reactions can happen. The potassium dichromate dissolves 
the metallic silver created in the development and the potassium bromide re¬ 
converts the dissolved silver into silver halides. The solution needs to be stirred 
for 3min. 

4. Fixing . The hxer removes the silver salts, the unmodihed gelatin (not exposed) 
and the supercoat. Instead, it leaves the transparent base and the modihed 
gelatin whose thickness is directly proportional to the exposition - the most ex¬ 
posed areas have less modihed gelatin removed. Usually, before the hxer is poured 
into the tank, the him is rinsed with pure water to remove all the substances left 
by the bleaching process. Then, the hxer is stirred in the developing tank for 
l-2min. Two diherent hxers have been used yielding similar results: 

(a) 30ml Kodak T-Max hxer -|- 300ml pure water 

(b) 45g sodium thiosulfate -|- 325ml pure water 
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Figure B.5; Phase-only CGH fabrication process. Grey (red) indicates that the film 
has (not) been exposed, a) Formation of latent image - silver halides become ionized, b) 
Developer converts ionized silver halides into metallic silver, c) Bleaching modifies the gelatin 
and converts the metallic silver back into silver halides, d) Fixing removes the supercoat and 
the unexposed emulsion. 


After the fixing is done, the him is completely transparent and light-resistant. 
It is important to rinse the him for 2min after the hxer has been applied. Any 
residual hxer could corrode hologram. Then, it only needs to be dried and cut. 


Figure B.5 summarizes the whole developing process explained in this section. In 
Figure B.5(a), the latent image is formed after exposure (Ag^Br’*'). Figure B.5(b) 
shows that the developer converts the ionized silver (Ag“) into solid metallic Ag. Then, 
Figure B.5(c) displays the action of the stopping bath -|- bleaching. The metallic Ag 
is precipitated and there are cross-links in the gelatin molecules. The modihed gelatin 
will remain after applying the hxing bath. Finally, Figure B.5(d) shows the action of 
the hxing bath, i.e. removing the supercoat and the unexposed emulsion, while leaving 
the modihed gelatine and the transparent base. A pattern of thickness given by the 
original picture is obtained as a result. 
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Hologram # 

Figure B.6: Efficiency of a batch of holograms as a function of shutter speed and aperture 
size. The holograms are the oblique phase given at Figure B.l. The dashed line is a guide of 
the eye. 

B.2 Characterization of CGHs 

The performance of holograms has been characterized in terms of the shutter speed and 
aperture size of the camera. All the rest of variables, i. e. printing method, photographic 
him, illumination and chemical process have been previously optimized, so they will 
remain constant in this study. All the holograms characterized in this section have 
been imprinted with the oblique phase shown in Figure B.l has been used. Then, the 
efficiency has been computed for many different holograms. Every single hologram has 
been characterized by the shutter time and the f-stop used to illuminate it. The results 
are shown in Table B.2 and Figure B.6. As it is observed in Figure B.6, for every single 
shutter speed, there is an optimum aperture size that maximizes the efficiency. For very 
slow shutter speeds, there is a very strong efficiency dependence on the aperture of the 
camera, whilst this dependence is much more smooth for short expositions. The results 
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also suggest that slow shutter speed can produce higher efficiencies. Furthermore, it is 
seen that the maximum efficiencies are above 40%, in accordance to the typical results 
in the held [207] . Nonetheless, it is worth commenting that higher efficiencies could be 
achieved if the chemical process was more accurate. Indeed, even though the phase- 
holograms are meant to be transparent, the ones used to carry out this characterization 
have absorption. If the chemical process was improved and the absorption was removed 
(maintaining the same thickness prohle), then the efficiencies could reach values of 
60%. In fact, this is the result that is obtained when the efficiency is computed as a 
ratio between the intensity going to the hrst diffraction order and the intensity that is 
transmitted through the CGH. 

B.3 Holograms - separation of diffraction orders 

Here, the experimental basics to separate and isolate the 1st diffraction order of a holo¬ 
gram are explained^. In section 6.1.1, a method to separate the different contributions 
from a hologram has been explained. Indeed, if an oblique phase ex.]){ikxx) is added to 
the hologram expression, the different diffraction orders of a hologram are separated 
by an angle = X/Ax, where Ax = 21^jk^. Now, this implies that at a distance ^ 
from the hologram the different diffraction orders are separated a transverse distance: 


Xq = ztan(A0x) ~ 


(B.l) 


Nevertheless, the width of the diffraction orders also grows following Gaussian optics. 
That is, if the beam waist at the hologram plane is tco, then the waist of every single 
diffraction order at a distance z from the hologram is 



(B.2) 


where I is the order of the phase singularity given by the hologram. It is clear that the 
condition that assures that two diffraction orders are separated is 



(B.3) 


^Note that all the equations given here are valid for both CGHs and SLMs. 
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Figure B.7: Sketch of a 4-f filtering system. The field at the object plane is Fourier- 
transformed by the hrst lens. At the Fourier plane, the diffraction orders split and a filter 
(either an iris or a pinhole) selects the 1st diffraction order. This order is re-imaged by the 
second lens. The hgures has been copied from [45] with the permission of Prof. M.C. Teich 


Note that equations (B.1)-(B.3) are valid at distances such that a free-space Fourier 
Transform of the held at the plane of the hologram is carried out. In this case, the 
complex amplitude of the plane wave components at the plane of the hologram are 
related to different positions in the Fourier plane. Experimentally, a free-space Fourier 
Transform is carried out when the distance to the hologram is 

2'kw‘^ 


z > 


(B.4) 


A 


2'kw 

Ax 


-4 


Due to lack of space or need to make set-ups compact, it is not always possible to 
propagate beams. However, there is another way of doing a Fourier Transform, and 
that is using a lens. A widely used set-up to separate the diffraction orders of a 
hologram is a 4-f hltering system. A sketch of a 4-f hltering system is shown in Figure 
B.7. On one hand, the diffraction orders are separated a distance xo at the Fourier 
plane: 


To = 


feA/ A/ 


(B.5) 


2tt Ax 

where / is the focal distance of the lens. On the other hand, the waist of a diffraction 
order at the focus of a lens is given by Gaussian optics: 


w = 


fxVi 


TTWa 


(B.6) 
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Figure B.8: Efficiency dependence on the polarization angle. The angle formed by the 
normal of the incidence and the beam axis is 8’. The red markers display the efficiency in 
the Y axis, whereas the blue markers display the efficiency in the X axis. 


Similarly to the free-space case, the condition that needs to be fulhlled is Xq S> tc, 
which yields the exact same condition given by equation (B.4). A rule of thumb to 
separate the diffraction orders is: 

Ax < ^ (B.7) 

B.4 SLM characterization 

Here, the efficiency of the SLM is characterized as a function of some controllable pa¬ 
rameters in the laboratory, such as the polarization of the incoming beam, the angle 
formed between the incident beam and the normal to the SLM, and the number of 
271 phase jumps displayed at the hologram (see equation (6.17)). Contrary to CGHs 
made of photographic hlms, the performance of the SLM is very dependent on all these 
parameters. Figures B.8 and B.9 show the influence of the polarization and the angle 
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Figure B.9: Efficiency dependence on the polarization angle. The angle formed by the 
normal of the incidence and the beam axis is 5 degrees. The red markers display the efficiency 
in the Y axis, whereas the blue markers display the efficiency in the X axis. 



Figure B.IO: Efficiency dependence on the polarization angle. The angle formed by the 
normal of the incidence and the beam axis is —14 degrees. The red markers display the 
efficiency in the Y axis, whereas the blue markers display the efficiency in the X axis. 
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of incidence on the efficiency^. In fact, in each of the Figures, two lines of data points 
are plotted. The red markers correspond the efficiency measured when the phase holo¬ 
gram is y) = exp(i/c|/), while the blue markers are obtained with y) = exp(ifca;). 
That is, the first (red) one measures the efficiency steering the beam in the x axis, and 
the second one in the y axis (the y axis is perpendicular to the plane of the optical 
table). For all the cases, the vertical axis displays the efficiency of the SLM measured 
in the way it has been described above; whilst the horizontal axis displays the polar¬ 
ization angle. This polarization angle is a direct readout of the angular position of 
a rotation-mount where a half-wave plate is hold. Hence, the value is arbitrary, as a 
rotation of the SLM would make it vary. However, the number is given so that it can 
be checked how the maximum is displaced when the SLM is tilted. The first conclusion 
that can be drown from Figures B.8 and B.9 is that the efficiency significantly depends 
on the polarization of the incident beam. Second, it is observed that when the angle 
between the normal to the SLM and the beam axis increases, the efficiency decreases 
more rapidly for a same range of variation of polarization degrees. That is. Figure 
B.8, which corresponds to almost normal incidence, shows that a 10 degree variation 
in polarization is matched with a 5% variation in efficiency. Figure B.9, which depicts 
the efficiency data obtained when the angle between the incident beam and the normal 
to the SLM is 5 degrees, shows that a 10 degree variation in polarization is matched 
with a 12% variation in efficiency. Furthermore the maximum value for the efficiency 
drops a little and changes its position. The position of the maximum varies 5 degrees 
approximately, which is approximately the angle that the SLM has been rotated. In 
Figure B.IO, where the results obtained when the angle between the normal to the 
SLM and the incident beam is —14 degrees are displayed, the same trends commented 
previously are maintained. That is, the efficiency of the maximum decreases; the angle 
for which the maximum of efficiency is moved as much as the angle between the normal 
of the SLM and the incident beam; and the gap between the efficiency in the x and y 
axis increases. The pixels are squares, so one would expect that the efficiencies in both 
axis were the same. However, a tilt of the SLM breaks the symmetry between the two 

■^The hologram used in both figure gives twenty-four 27r phase jumps. 
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axis and accounts for the measured difference. 

Next, the importance of the hologram displayed on the SLM is characterized. In 
Figures B.8-B.10 the same hologram was used. The hologram is the oblique phase in 
Figure B.l for the efficiency in the y axis, and a 71/2 rotation of it for the x axis. This 
hologram contains twenty-four 27r phase jumps. In Figure B.ll the efficiency is plotted 
as a function of the number of 27r phase jumps given by the hologram. The angle of 
incidence is set at 8’ (same as Figure B.8) and the optimum position of the half-wave 
plate angle is set to maximize the efficiency. Obviously, as given by equation (6.17), the 
larger number of 27r phase jumps the hologram gives, the larger is the angle between 
the 0th and the 1st diffraction order. Then, it can be observed that there is an almost 
linear relation between the efficiency of the SLM and the number of 27r phase jumps in 
for this range of values. The larger the number of 27r phase jumps given by the holo¬ 
gram, the lesser pixels per ramp are used, and consequently the efficiency is reduced. 
To hnalize, the Lookup Table (LUT) provided by the producer of the SLM (BNS) was 
tested. All the data points collected to plot Figures B.8-B.11 were done at A = 632.8 
nm, using a He-Ne laser (see Figure 6.7). The phase shift given by an SLM highly 
depends on A. The transmittance function created by the liquid crystals depends on A 
as shown in equation (6.4). This means that the efficiency of the SLM can vary a lot 
depending on A. In order to account for this dependence and some other imperfections 
of the SLM such as non-linear behaviour of the pixels, LUTs change the theoretical 
value of pixels to give them an effective one which increases the performance of the 
hologram as a whole. In [268], a linear method to correct the imperfections of the SLM 
was developed. The idea is to multiply the phase of the phase-hologram y) by a 

contrast constant C. Then, when the new t'{x, y) = Ct{x, y) is such that t'{x^ y) < —tt 
or t'{x,y) > TT, those phases are re-set to —vr or tt. Using this method. Bowman et 
al. managed to greatly increase the efficiency of an SLM. Using a hologram with ten 
27r phase ramps and maximizing the polarization angle at an angle of incidence of 8’, 
I tried to use the same method to improve the efficiency of the SLM. The results are 
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Figure B.ll: Efficiency dependence on the number of 2'k phase jumps displayed on the 
SLM. The angle formed by the normal to the SLM and the beam axis is 8’. The polarization 
is set to optimize the efficiency. The red markers display the efficiency in the y axis, whereas 
the blue markers display the efficiency in the x axis. 


shown in Figure B.12. The y axis shows the efficiency and x axis the contrast param¬ 
eter C. Leaving (7 = 1 means not doing any changes to the LUT. It can be observed 
that the highest efficiency is achieved for C = 1.02. Hence, the LUT had been properly 
optimized at the frequency of the HeNe laser. I believe that whenever a new SLM or 
wavelength are used, this is a sanity check that should be done. In fact, some producers 
do not send LUTs along with the SLM. In those cases, it is highly recommended to 
apply the method described in [268] to improve the performance of the SLM. 
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Contrast 

Figure B.12: Efficiency dependence on the contrast. The hologram in consideration 
displays ten 27r phase ramps. The angle of incidence is 8’, and the polarization is optimized 
to maximize the efficiency. 
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Table B.2: Efficiency as a function of the aperture and the shutter speed. An aperture 
size of the kind 2.8-4 means that the aperture was set at an intermediate position between 
2.8 and 4. 


Hologram 

Shutter speed (s) 

Aperture size 

Efficiency (%) 

1 

30 

1.7 

35.3 

2 

30 

2.8 

12.4 

3 

15 

1.7 

44.9 

4 

15 

2.8 

29.2 

5 

15 

2.8-4 

13.2 

6 

8 

2.8 

44.0 

7 

8 

2.8-4 

40.3 

8 

8 

4 

23.1 

9 

8 

4-5.6 

12.3 

10 

4 

2.8-4 

40.6 

11 

4 

4 

43.2 

12 

4 

4-5.6 

34.8 

13 

4 

5.6 

20.3 

14 

4 

5.6-8 

8.6 

15 

2 

4 

36.1 

16 

2 

4-5.6 

39.4 

17 

2 

5.6 

36.5 

18 

2 

5.6-8 

29.5 

19 

2 

8 

14.8 

20 

2 

8-11 

7.0 

21 

1 

4 

27.7 

22 

1 

4-5.6 

32.0 

23 

1 

5.6 

35.9 

24 

1 

5.6-8 

37.2 

25 

1 

8 

43.1 

26 

1 

8-11 

38.3 




Proof of 



In this appendix, the equation (7.18) is demonstrated. In order to do that, the following 


relations given in chapter 7 are used: 

/ OO POO 

/ \EX„_i_,\^dxdy (C.l) 

■OO j —OO 

K,i = (^-2) 

d;-d_p = 0 (C.3) 

Using all of them, it can be obtained that can be expressed as: 

/ OO POO 

/ \Aly_^i{x,y)\‘^+ \Bly_^i{x,y)\‘^dxdy (C.4) 

-OO J —OO 

f? / T 

Following an identical procedure, the following equation yields for /_/ : 

/ OO POO 

/ \At^/^^^_i{x,y)\‘^+ \Bty+i^_iix,y)\‘^dxdy (C.5) 

-OO J —OO 


Choosing the mirror symmetric operator to be M{z} = then equation (7.11) 

can be used to get a relation between the coefficients in equations (C.4,C.5): 

= Al i{x, y)ap + 5*/x, y)d.p 

= (Alp _;(x, y)a.p + 51^ _,(a;, y)dp) (C.6) 

= -Alp^_;(-a;, y)dp - Btp^_^{-x, y)a_p 
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Proof of 


which implies that Api{x,y) = and Bpi{x,y) = B^p_i{—x,y), due to 

the orthogonality of dp and d_p. Consequently, it follows that: 


tL/R 




|dl^i/_i,/(a;,|/)p + \Bl^/_^^{x,y)\‘^dxdy 

I - At^/^^^i{-x,y)\'^ + I - Bti/+i^i{-x,y)\‘^dxdy 

\Aty+i^_i{x',y)\‘^ + \Bti/+i^_i{x',y)\‘^dx'dy = 


as the integrations limits remain the same under the change x —)■ —x'. 


(C.7) 



D 

Intensity plots of multipolar fields 


In the current appendix, intensity distributions of multipolar fields are displayed. Bear¬ 
ing their spatial shape in mind is of special importance to understand the excitation of 
WGMs (section 4.4), the suppression of backscattering (5.2), and the intensity profiles 
recorded by the CCD camera in the experiments done in chapter 8. Furthermore, it is 
fairly common to forget about the fact that there are four different kinds of multipolar 
fields depending on the radial function used: jj{kr), nj{kr), h^^\kr), and h^^^\kr) . 
Also, as it has been shown across this whole thesis, the multipolar fields can be clas¬ 
sified either in terms of parity or helicity. Hereafter, some frequently used multipolar 
fields are displayed. They are displayed for both parities, helicities, and typically used 
radial functions {jj{kr) and h^p{kr)). The plots are done in the z — x plane, 2 ; be¬ 
ing the horizontal axis and x the vertical one. Then, because they are cylindrically 
symmetric around 2 ;, the expression in the other planes can be deduced by applying 
a 27r revolution symmetry. The axis are normalised to the wavelength, so they are 
dimensionless. All the units are normalized to A. Also, due to the properties of the 
Wigner D-matrices Dl^^p{(j)k,9k) [39], the plots have only been done for > 0, as it 
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can be observed that the following properties hold: 



|A7 , 


2 


I A (y) |2 _ I A ( y ) |2 


(D.l) 

(D.2) 


First, the dipolar orders are displayed. It is apparent from Fignre D.l (and it can be 
re-checked with all the rest of Figures displaying multipoles with well-dehned helicity) 
that equation (8.2) holds. In fact, equation (8.2) is a consequence of the fact that 


= \A-^f (D-3) 

Second, it can also be observed that A^- has the same behaviour for the two z < 0 and 
z > 0 semi-spaces, whilst does not. This property is the basis of the measurements 
carried out in chapter 8. Third, the stationary wave behaviour of the Bessel function is 
only observed for but not for A^^^. Remember that a spherical Bessel function 

jj{kr) is a superposition of two spherical Hankel functions jj{kr) = h^p{kr) -|- h^^\kr), 
where the spherical Hankel functions have a respective outwards and inwards radiative 
behaviour (see equations (2.59, 2.60) in section 2.4). In addition, looking at all the 
captions of the different Figures, it is observed that the plot area does not remain 
constant. This has been done to visualize the spatial shape of the multipolar fields. 
Indeed, high order Hankel multipolar helds are highly evanescent, therefore the plot 
area has been reduced. In contrast, high order Bessel multipoles require a larger 
plot area. Also note that a central region surrounding the origin of coordinates has 
been removed from the plots of Hankel multipoles. The reason for this is that Hankel 
functions are singular at the origin. Finally, Figures D.11-D.16 depict multipolar beams 
with j = 15. Note the similarities between Figure D.ll and Figure 4.9 in section 4.4, 
where almost a single mode was excited. And last but not least. Figure D.16, which 
depicts A^g^, shows an intensity plot which highly resembles the typical plots for 
WGMs held distribution. It is interesting to see that those plots are achieved when 
rriz = 1 and j is large, in contrast to Oraevsky’s dehnition of WGM, which states that 
a WGM is a multipolar held with very large j and ruz = j [86]. That, of course, does 
not ahect the Q factor of the sphere, as the Q factor is only given by the linewidth of 
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Figure D.l; ^ and for p = ±1. The plot area is a square of side s = 127rA. A 
contour of O.SvrA around the origin has been removed to plot A^’^. 


a Mie resonance, usually the first one for WGMs [86]. And in Mie theory the value 
is degenerated, therefore a multipolar held has the same Q-factor as 
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(m) 


A 

"^ 1,1 




Figure D.2: and for (y) = (e), (m). The plot area is a square of side s = 127rA. 


{y),h 


A contour of O.SvrA around the origin has been removed to plot 



Figure D.3; A^q and A^’q for p = ±1. The plot area is a square of side s = 127rA. A 
contour of O.SvrA around the origin has been removed to plot A^’g. 











195 


(e) (///) 



Figure D.4: and for (y) = (e), (m). The plot area is a square of side s = 127rA. 

A contour of O.SvrA around the origin has been removed to plot A^^q’^. 


p = 1 p = —1 



Figure D.5: A 22 and A^’^ for p = ±1. The plot area is a square of side s = ISvrA. A 
contour of O.SvrA around the origin has been removed to plot A^’^- 
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(e) (/») 



Figure D.7: and for p = ±1. The plot area is a square of side s = ISvrA. A 

contour of O.SvrA around the origin has been removed to plot A^’^. 
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"^ 2,0 


p = 1 


11 


p = —1 


11 



Figure D.9; A^q and A^’g for p = ±1. The plot area is a square of side s = ISvrA. A 
contour of O.SvrA around the origin has been removed to plot A^’q. 
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p = 1 p = —1 



Figure D.ll: and for p = ±1. For A^g ^g, the plot area is a square of side 

s = SdvrA. For A^g^^g, the plot area is a square of side s = GvrA, where a contour of l.GvrA 
around the origin has been removed. 
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Figure D.12: and for (y) = (e), (m). For the plot area is a square 

of side s = SGttA. For the plot area is a square of side s = GvrA, where a contour of 

O.SttA around the origin has been removed. 


P = 1 p = —1 



Figure D.13: A^gg and for p = ±1. For A^^g, the plot area is a square of side 

s = SGttA. For A^g^g, the plot area is a square of side s = GttA, where a contour of I.GttA 
around the origin has been removed. 
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(e) (i/O 



Figure D.14: and for (y) = (e), (m). For A^|^g, the plot area is a square of 

side s = SGttA. For A^|^g^, the plot area is a square of side s = GvrA, where a contour of O.SvrA 
around the origin has been removed. 



Figure D.15: A^^ ^ and A^^^^ for p = ±1. For A^g the plot area is a square of side 

s = SGttA. For the plot area is a square of side s = d.SvrA, where a contour of I.GttA 

around the origin has been removed. 
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Figure D.16: and for (y) = (e), (m). For the plot area is a square of 

side s = SGttA. For the plot area is a square of side s = d.S-TrA, where a contour of 

O.SttA around the origin has been removed. 
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List of Symbols 


Fields 


b’(x, t). general EM field 

V{r,t) . electric displacement 

S{r,t) . electric field 

B{r, t) . magnetic induction 

. magnetic field 

V{r,t) . polarization density 

Ai{r,t) . magnetization density 

D(r) . complex-amplitude electric displacement vector 

E(r) . complex-amplitude electric field vector 

B(r) . complex-amplitude magnetic induction vector 

H(r) . complex-amplitude magnetic field vector 

E^°'^. electric field at the focus of a lens 

Eoo(0fc, . angular spectrum of the electric held 

E“. incoming electric held 

(E®, H®) . scattered EM held 
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List of Symbols 


(E\H'). incident EM field 

. interior EM field 

(E*,H*). total exterior EM field 

Ep". electric field of an incident vortex beam with topological charge I and 

helicity p 

Ep;. transmitted electric held when the incident beam is E^" 

E®. direct helicity component of the scattered electric held 

Eip. crossed helicity component of the scattered electric held 

Api{x,y) . complex amplitude of the direct component of the transmitted elec¬ 

tric held when the incident beam is El", 

Bp i{x, y) . complex amplitude of the crossed component of the transmitted elec¬ 

tric held when the incident beam is El", 

PA 

Mathematical functions 


(5(k — k'). Dirac delta 

Sppi . Kronecker delta 

Jn{kpp) . Bessel function of the hrst kind and order n 

Ok) ■ ■ ■ ■ Wigner rotation matrix 

jj{kr) . spherical Bessel function of the hrst kind and order j 

. generalized Legendre polynomial, p = cosO 

hj{kr) . spherical Hankel function of the hrst kind and order j 

H{6k — 0^) .... Heaviside step function 
dln*p{0k) . reduced rotation matrix 
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r(; 2 ). Gamma function 

2 -Fi(a, b,c : z)... hypergeometric function 
Modes 


(p, . 

* ‘ * 

jrriz 1 jiriz . 

. 

\{m) Ah) 

juiz 1 jmz . 

A + A “ 

jm,z 1 jnriz . 

bX. 

Ah) 

. 

Operators 


(p,s) plane-waves 

plane wave with well-dehned helicity (-1-1, —1) 
Hankel multipolar fields (with well-dehned parity) 
Hankel multipolar helds with well-dehned helicity 
Bessel beam with helicity p = ±1 
multipolar helds (with well-dehned parity) 
multipolar helds with well-dehned helicity 
Bessel beam with (p) = TE/TM mirror symmetry 
Whispering Gallery Mode 


H . hamiltonian 

Ta* . temporal translation (or evolution) of magnitude At 

P. linear momentum vector operator 

Pa;, Py, Pz . cartesian linear momentum operators 

Tap . linear translation of magnitude Ar 

J. angular momentum vector operator 

S. spin angular momentum vector operator 

L. orbital angular momentum vector operator 
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List of Symbols 


Jx, Jy, Jz . cartesian angular momentum operators 

Rii{(p) . rotation operator of an angle ip with respect to an axis given by n 

. rotation matrix of an angle p with respect to an axis given by n 

Ok) . composition of three rotation operators. 

R(0fc,6'fc) = Rz{(j)k)Ry{0k)Rz{0) 

. angular momentum squared operator 

A. helicity operator 

D^p . generalized duality transformation of an angle p 

n. parity operator 

. Mirror symmetry operator with respect to a plane with normal vector 

n 

M{n}. Mirror symmetry operator with respect to a plane containing the 

vector n 

Coordinates and vectors 
t . time coordinate 

Y = {x,y,z) .... position vector in cartesian coordinates 

Y = {p,(j), z) .... position vector in cylindrical coordinates 

r = {r,(j),6) . position vector in spherical coordinates 

X, y, z. cartesian unitary vectors 

p,4>,z . cylindrical unitary vectors 

Y, (j),0 . spherical unitary vectors 

k = {kx,ky,kz) . wave-vector in cartesian coordinates 


















207 


k = {kp,(j)k,kz) ■ wave-vector in cylindrical coordinates 
k = {k,(j)k,0k) ■ ■ wave-vector in spherical coordinates 

n. unitary vector normal to a surface 

(T+, d-. left and right circular polarisation unitary vectors 

Beam and material parameters 

A. wavelength 

u . frequency 

c. speed of light 

k . wave-number 

LGi^q . Laguerre-Gaussian beam 

g>0. radial index 

I . topological charge (azimuthal index) 

w{z) . beam waist at a plane of 2 : constant 

Wo. beam waist at z = 0 

p . helicity value 

eo,/io. electric permittivity and magnetic permeability of vacuum 

e, fi,n . electric electric permittivity, magnetic permeability, and refractive 

index of a linear, non-dispersive, homogeneous and isotropic medium 

a . charge density on a surface 

js. surface current density 

S(g, m). scattering matrix as a function of the geometrical and material prop¬ 

erties 
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List of Symbols 


/ . 

GLMT 

ei, /ii, Hi 


Ur = rii/n. 

R . 

2'kR 

T=— 

^ji^j . 

. 

42 > 41 - 
(s^). 

a. 

C'.. 

Cext . 

Qi . 

Qs . 

Q ext . 

'-'jrrizp . 

e. 

Q. 


focal distance of a lens 

Fresnel coefficients of the lens for s' and p waves 

electric electric permittivity, magnetic permeability, and refractive 
index of sphere 

relative refractive index 

radius of sphere 

size parameter 

scattering Mie coefficients of order j 

interior Mie coefficients of order j 

beam shape coefficients 

time-averaged total exterior Poynting vector 

interior cross section 

scattering cross section 

extinction cross section 

interior efficiency factor 

scattering efficiency factor 

extinction efficiency factor 

beam shape coefficients for beams with a well-defined helicity 
root number for a Mie coefficient 
Q factor 
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w. energy density 

W. energy of the EM per unit length 

W®. energy scattered 

Wp. energy scattered in modes with the same helicity as the incident beam 

Wip. energy scattered in modes with opposite helicity to the incident held 

Tp{R,nr) . transfer function 

Tmip{R,nr) . transfer function for a cylindrically symmetric beam with = nil 

Holography 

U{r,t) . complex wave-function 

U{r) . complex amplitude 

/(r). intensity 

t{x,y) . complex amplitude transmittance 

d{x, y) . thickness of an optical element 

Uo . object wave 

Ur . reference wave 

Ax . width of a 2ti phase ramp 

Ax^ . width of a 2ti phase ramp for a given A 

A9x . diffraction angle given by a hologram 

Experimental parameters 


$. diameter 

T . target 
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List of Symbols 


T. target integro-differential operator 

CDi(%). CD given an incident vortex beam of order I 

tr . Extinction ratio of a linear polariser 

7 . ratio of light scattered into crossed helicity component divided over 


the direct one 

iP . transmitted intensity for a RCP incident vortex beam of topological 

charge I 

. transmitted intensity for a LCP incident vortex beam of topological 

charge I 

P . measured intensity of the incident EM held 

Jlp. measured intensity of the crossed helicity component of the scattered 

EM held 


Ip . measured intensity of the direct helicity component of the scattered 

EM held 


. normalised (by the background) measured intensity 

(TAr(%). standard deviation with respect to the average 

. normalised (by the background) measured intensity at A given a 

monochromatic incident beam with = m, A = p 
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